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"^ Abstract 

pg I Starting from hyperbolic dispersion relations, we derive a closed system of Roy-Steiner equations 

for pion-nucleon scattering that respects analyticity, unitarity, and crossing symmetry. We work 

out analytically all kernel functions and unitarity relations required for the lowest partial waves. 

In order to suppress the dependence on the high-energy regime we also consider once- and twice- 

Qji subtracted versions of the equations, where we identify the subtraction constants with subthreshold 

parameters. Assuming Mandelstam analyticity we determine the maximal range of validity of these 

equations. As a first step towards the solution of the full system we cast the equations for the 

pH , TTTT — > NN partial waves into the form of a Muskhelishvili-Omnes problem with finite matching 

point, which we solve numerically in the single-channel approximation. We investigate in detail 

-- 1 I the role of individual contributions to our solutions and discuss some consequences for the spectral 

^ I functions of the nucleon electromagnetic form factors. 
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1 Introduction 

Pion-nucleon scattering is one of the most basic and fundamental processes in strong-interaction 
physics. Even though a large data basis exists and numerous investigations based on a cornucopia 
of methods (dispersion relations, quark models, resonance models, chiral perturbation theory, just to 
name a few) have been performed for many decades, the pion-nucleon (vrA^) scattering amplitude is 
still not known to sufficient precision in the low-energy regionj^ This becomes most obvious in the 
scalar-isoscalar sector, which features the so-called pion-nucleon a term (TttAt, i-e. the scalar form factor 
of the nucleon at zero momentum transfer. Its value is a measure of the light quark contribution to 
the nucleon mass (and it can also be related to its strange quark contribution), see e.g. the classical 
paper |5j. The a term has gained renewed interest as it parameterizes the spin-independent cross 
section for possible dark matter candidates scattering off nuclei [BUT] (for a recent review cf. |8| ) . In 
principle, lattice QCD would be the method of choice to pin down the a term — however, a direct 
computation of the scalar form factor necessarily involves disconnected diagrams, which is not yet 
under sufficient control. Similarly, the indirect extraction of a-j^N from the derivative of the nucleon 
mass is still hampered with systematic uncertainties related to the chiral extrapolations utilized, see 
e.g. [9]. Therefore, in this paper we follow a different path, namely setting up the powerful machinery of 
Roy-Steiner (RS) equations that will ultimately allow for a precise determination of the pion-nucleon 
scattering amplitude at low energies. 

More specifically, RS equations are based on hyperbolic dispersion relations (HDRs), a particular 
kind of dispersion relations along hyperbolae in the Mandelstam plane. Dispersion relations are a 
widely used tool that is built upon very general principles, such as Lorentz invariance, unitarity, 
crossing symmetry, and analyticity. There are multiple uses of dispersion relations — they can be used 
to stabilize extrapolation of experimental data to threshold and allow for a continuation into unphysical 
regions, as it is e.g. required for the extrapolation of the pion-nucleon scattering amplitude to the 
so-called Cheng-Dashen point [TO], which is crucial for the extraction of the a term. We notice that 
unitarity constraints can most conveniently be formulated in terms of partial- wave amplitudes. The 
resulting partial-wave dispersion relations (PWDRs) together with unitarity constraints allow to study 
processes at low energies with high precision. We just mention a few examples. The most prominent 
example is of course pion-pion (vrvr) scattering, which is intimately linked to the spontaneous and 
explicit chiral symmetry breaking in QCD. The Roy equations [llj are the appropriate PWDRs, which 
have been extensively studied in the last years |12H17j . leading to a determination of the fundamental 
TTTT scattering amplitude with unprecedented precision. The pion-pion system, however, is special 
as all channels are identical. This is different for the simplest scattering process in QCD involving 
strange quarks, namely pion-kaon [ttK) scattering, which has been investigated in [181119]. As far 
as crossing symmetry and isospin quantum numbers are concerned, the pion-kaon system is similar 
to the pion-nucleon case considered here. Crossing symmetry relates the s-/ii-channel (vrA^ — t- ttN) 
and the t-channel (tttt — )• NN) amplitudes, with the s-channel amplitudes relevant e.g. for cr-term 
physics, while the t-channel amplitudes feature prominently in the dispersive analysis of the nucleon 
form factors. The final aim of solving the full (subtracted) RS system for ttN scattering is a precise 
determination of the lowest partial- wave amplitudes in the low-energy (physical and unphysical) region 
as well as the pertinent low-energy parameters, such as the vrA^ coupling constant and the so-called 
subthreshold parameters, and to provide reliable theoretical errors for the fundamental pion-nucleon 
scattering amplitude for the first time. 

In the low-energy region, the pion-nucleon amplitude is well represented by its S- and P-wave 
projections. Due to the spin of the nucleon, one has in total six partial waves in the s- and u-channel, 



^The exceptions are the S-wave scattering lengths, which can be extracted with high precision from the beautiful 
data on pionic hydrogen and pionic deuterium, see [IHlj- 



commonly denoted as /q_,_, f^_^_, /^_, where the superscript / = it refers to the isospin, / G {0, 1} in the 
subscript to the orbital angular momentum, and the it to the total angular momentum j = / it 1/2. 
Similarly, there are three t-channel S- and P- waves, called f^, /j_, where the superscript refers to 
total angular momentum J and the +/— to parallel/antiparallel antinucleon-nucleon helicities, such 
that there is one wave with even and two with odd isospin (due to Bose symmetry). It was pointed 
out in [20] how to generalize the Roy equations for vrvr scattering to the vrA^ system based on fixed-t 
dispersion relations. These amount to coupled integral equations for the nine partial waves, where the 
effect of the higher partial waves is encoded in the respective kernels of these integral equations. Here, 
we follow a somewhat different path by utilizing hyperbolic dispersion relations as pioneered by Hite 
and Steiner a long time ago [21] . The main advantage of HDRs is that they combine the s- and the 
t-channel (i.e. all three) physical regions, which is obviously not true for e.g. usual fixed-t dispersion 
relations. It is known that a reliable continuation to the subthreshold region in dispersion theory 
can only be made by using input information also from the t-channel, cf. e.g. |22H25j . Furthermore, 
the knowledge of the absorptive parts in the dispersion relations is needed only in regions where the 
corresponding partial-wave expansions converge, and HDRs are considered the best choice fulfilling 
these requirements that yields still manageable angular kernels [21j . In addition, the underlying 
hyperbolic relation {s — a){u — a) = b (with a, b real- valued parameters) also respects s -H- u crossing 
symmetry of the ttN amplitude. Due to the tunable parameters a, b, better convergence properties 
can be achieved with HDRs and they are found to be especially powerful for determining the a 
term ^22j. The derivation of the RS equations for the vrA^ system is given by a series of steps: first, 
one expands the s-/t-channel absorptive parts of the HDRs in s-/t-channel partial waves, respectively. 
Second, one projects the full, partial-wave-expanded HDRs onto both s- and t-channel partial waves, 
resulting in what we will refer to as the s- and t-channel part of the RS system in the following. The 
resulting system of equations exhibits the following general structure: it features the nucleon-pole- 
term contributions, integrals over the imaginary parts of the s-(and ti-)channel as well as integrals over 
t-channel absorptive parts, both from the corresponding threshold to infinity. The generic properties 
of the equations are then determined by the integral kernels. In the equation for each partial wave, the 
corresponding kernels consist of the self-coupling, singular Cauchy kernel and an analytic remainder 
that in addition involves the coupling to all other partial waves. In particular, these kernel functions 
automatically incorporate the analytic properties expected for a given partial wave: the Cauchy kernel 
corresponds to the right-hand cut, while the remainder contains all left-hand-cut contributions. 

Another important issue is the possibility to subtract dispersion relations. This can be advanta- 
geous for various reasons: first, in some cases the asymptotic behavior of the integrand is such that 
subtractions have to be performed to ensure convergence of the dispersive integral. Similarly, if the 
high-energy behavior is not known, it can be subsumed in subtraction constants, which are a priori 
unknown. In some cases, these subtraction constants can be related to phenomenology or the parame- 
ters of a low-energy effective field theory like e.g. chiral perturbation theory (ChPT). Second, one can 
even introduce subtractions that are not necessarily required by the asymptotic behavior in order to 
lessen the dependence on high-energy input, however, at the expense of introducing the corresponding 
subtraction polynomials. Third, subtracting the dispersion relations is especially useful in the ttN 
case, since subtracting at the so-called subthreshold point allows for a relation to the subthreshold 
expansion and is convenient for the continuation to the Cheng-Dashen point. In addition, such sub- 
tractions are well suited for the t-channel problem to be discussed later. In what follows, we will 
consider unsubtracted as well as subtracted versions of the RS equations. 

Next, we will outline the strategy to solve the RS equations, as depicted in Fig. [1) first, one 
solves the t-channel part of the RS system, which takes the form of a Muskhelishvili-Omnes (MO) 
problem [26p27] (using rather well known s-channel partial waves and vrvr phase shifts as input). Then, 
one uses the t-channel MO solutions to solve the s-channel part, and finally the procedure is repeated 
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Figure 1: Flowchart of the solution strategy for the Roy-Steiner system for ttN scattering. 



(iterated) until self-consistency of the partial waves and parameters is reached and the results have 
converged, cf. Fig. [TJ In both the s- and t-channel part of the system one actually solves the equations 
in the low-energy region and for the lowest partial waves, while the amplitudes in the high-energy region 
as well as higher partial waves are needed as input. The separation between both energy regions occurs 
at the so-called matching points Sm and t^i in the s- and f-channel, respectively. Due to the complexity 
of the full problem, we will not yet solve the whole set of RS equations in this article, but concentrate on 
the t-channel part of the system as a first step. The solution of this t-channel subproblem is interesting 
by itself, as it features in the dispersive analysis of the nucleon electromagnetic form factors as well as 
the scalar form factor, which is, in turn, essential for the extraction of the a term. At present, in the 
unphysical region only the KH80 solution |28p29j has been used. It is, however, well-known that this 
solution does not include more recent and precise data and that the nN coupling constant used there 
differs significantly from more modern determinations. Furthermore, no analysis of the theoretical 
uncertainties is performed (apart from an iteration uncertainty, cf. Sect. [53^ . which is an absolute 
requirement for any modern theoretical investigation. Therefore, a new t-channel solution is needed 
as a first step for solving the full system. Finally, a consistent set of partial-wave amplitudes for all 
channels is especially important as far as the u-term extraction is concerned, and it has been pointed 
out that the KH80 solution seems to suffer from internal inconsistencies |22 p 25 t [30]. which emphasizes 
the necessity of a full system of PWDRs. 

The original Roy equations for vrvr scattering |llj were solely based on fixed-t dispersion relations. 
This approach fails for processes involving non-identical particles, since crossing symmetry intertwines 



different physical processes. For this reason, a combination of fixed-t and hyperbolic dispersion re- 
lations was used in [I8l[l9] to construct integral equations for -kK scattering for the partial waves of 
both s- and i-channel, which are therefore referred to as Roy-Steiner equations. In this work, we 
solely consider HDRs, a path that has already proven useful in the construction of RS equations for 
77 — )• TTTT [31]. Our solution strategy for the t-channel MO equations follows [19], however, there is a 
major difference between vrvr — >• KK and vrvr — )• NN as far as inelasticities in the unitarity relation 
are concerned, since the pseudophysical region in the ttN case is much larger due to the large nucleon 
mass. In both cases, the first non-negligible contribution besides vrvr intermediate states originate 
from KK, which play an important role for the 5- wave in view of the occurrence of the /o(980) res- 
onance. For vrvr — )• KK the inelasticities can simply be accounted for by using phase-shift solutions 
for the corresponding partial waves, while physical input for tttt — t- NN is only available above the 
two-nucleon threshold. Once the t-channel problem is solved, the remaining equations take the form 
of the conventional vrvr Roy equations, such that known results concerning existence and uniqueness 
of solutions [32l|33] may be transferred to the s-channel RS equations as well. 

This work is organized as follows: in Sect. [2] we specify our conventions and review HDRs for the 
invariant amplitudes of vrA^ scattering. In Sects. [3] and H] we derive a closed system of RS equations as 
well as a once- and twice-subtracted version, and show how the i-channel equations can be cast into 
the form of a MO problem. Sect. [5] is devoted to the explicit solution of the i-channel MO problem: 
we first review the MO problem with a finite matching point and state the explicit solution for the 
ttN t-channel amplitudes. Then we collect the necessary input and discuss the numerical results. 
Finally, we briefiy discuss the application to nucleon form factors before concluding in Sect. [6j The 
explicit derivation of the s- and i-channel RS equations is described in full detail in Appendices O 
and [HI respectively. In Appendix [C] we determine the range of convergence of our equations, while 
Appendix IDJ contains a discussion of the asymptotic regions in the dispersion integrals. 

2 Preliminaries 

2.1 Kinematics 

We take the s-channel reaction of -kN scattering to be 7r(g) -|- N{p) — )■ 7r(g') -|- N{p') and the i-channel 
reaction to be 7r(g) -|- 7r(— g') — )• N{—p) + N{p') with the usual Mandelstam variables 

s = {p + qf, t = ip-p'f, u={p-q'f, (2.1) 

which fulfill 

s + t + u = Irr? -F 2MI = E , (2.2) 

where m and M,r denote the nucleon and pion mass, respectively. We will use the masses of |34| . with 
the isospin limit defined by the charged particles, i.e. M^^ = M^± and m = rrip (later also Mk = Mj^± 
for the kaon mass). Unless stated otherwise, u is always to be understood as a function of s and t 

u{s,t) = T.-s-t . (2.3) 

We define the generic kinematical Kallen function 

Af« = A(x,M|.,M2) = [x-(Mp-MQ)2][x-(Mp + MQ)2] , (2.4) 



and for the equal-mass case 
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Furthermore, we introduce the general definitional 
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with W as the total center-of-mass-system (CMS) energy, as well as the abbreviation 

K = y/ = \{x,m^,Ml) = [x-sJ\[x-s+], s± = Wl = {m±M^)\ (2.7) 

where W- and Wj^ denote the (s-channel) pseudothreshold and threshold energies, respectively. Ad- 
ditional related useful definitions and relations are 

T.± = m^ ±mI , S+ = So , S_ = W+W- , S^ = s+S- , T. = s+ + s^ . (2.8) 

The CMS kinematics of the elastic s-channel reaction vrA^ — )• ttN above threshold (i.e. for s > s+) 
with CMS momentum q = |q|, nucleon energy E, and scattering angle Zg = cos 9s are then given by 

E{±W) = ±Vm^ + q^ = |±|^ = ±EiW) , 

V = s _ E + ^ . (2.9) 

s 

For the t-channel reaction vrvr — )■ NN with CMS momenta qt for the pions and pt for the nucleons 
and scattering angle zt = cosOt, the CMS kinematics above threshold (i.e. for t > Am?') read 



'^(^) = V l^ ' 




s + u-Y. 
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qt{t) = ^/ 1 - M2 = ^^r = +i(Z- , pt{t) = \l^--m^ = ^af = +ip. 



V4 ^'= 2^* ^ 
s-n 2s + t-S 


= +iq- 
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ztis,t) = = — = , (2.10) 

Aptqt Aptqt Ptqt 

where below the corresponding two-particle thresholds t.,^ and ijv one has to use the quantities 



_(t) = JM^-j>0 \/t<t^ = AMl , p_{t) = Jm'^ -j>0 \/t<tN =Am^ , (2.11) 

whose phases are constrained in general to ptqt = —p-q- and fixed here by convention. Relations 
valid in all kinematical ranges can be written down by relying on the quantities 

i{t) = '-^ = -ql{t) , p?(t) = ^^ = -pHt) , (2.12) 

from which roots in the corresponding regimes may be taken |f| 

^For more on ttN kinematics and for ttN conventions in general we refer to [29]. Note that the convention for ly 
therein and which we have adopted here differs from the choice u — s — u of e.g. [2T] . 

•^We use the non-cyclic convention J^xJ^ for a reaction a + b—^c + d in order to stick to the usual ttN conventions of [29| . 
rather than the cyclic convention ^xf that leads to symmetric kinematical relations for the s-, t-, and u-channel and is 
therefore sometimes used in the literature. While the cyclic convention is especially favorable when all four particles are 
identical like e.g. in the case of tttt scattering, it leads to different sign conventions for the CMS scattering angles and also 
to different isospin crossing matrices (cf. Sect. 12.2] ). The non-cyclic convention, however, is well-suited for s -ir^ u crossing 
symmetric situations like e.g. ttN scattering, with t — corresponding to an undeflected pion (i.e. forward scattering) 
in both the s- and it-channel and thus Zs(t — 0) = 1 = Zu{t = 0) rather than Zs(f = 0) = 1 = — Zu(i — 0) for the cyclic 
convention. 
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Figure 2: Physical regions for s-, t-, and n-channel reactions of vrA^ scattering (shaded) and the 
subthreshold triangle (dot-dashed) enclosing the subthreshold lens. 



The physical regions for the s-, t-, and n-channel reactions are restricted to kinematical regions 
where the Kibble function ^ [35j is non-negative. For vrA^ scattering we have 
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such that the boundaries are given by 
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and the corresponding physical regions are shown in Fig. [2j 

ttN scattering in the isospin limit can be described by the four Lorentz-invariant amplitudes 
A [s,t) and B {s,t), as well as the related amplitudes D {s,t) convenient for low-energy theo- 
rems (all to be defined in Sect. 12. 2|) . These amplitudes are real inside the Mandelstam subthreshold 
triangle defined by the lines s = 5+, u = s+, and t = t^, i.e. below the thresholds for the physical s- 
and u-channel reactions and below the vrvr thresholdjj including in particular the small on-shell but 
unphysical lens-shaped low-energy region (subthreshold lens) close to (i^ = 0,i = 0) depicted in Fig. [2j 

The analytic structure of the invariant amplitudes governs the analytic structure of both the s- 
and t-channel partial-wave amplitudes, for details we refer to [29] and references therein. Here, we 
only mention the different analytic structures of the s-channel vrA^ scattering invariant amplitudes (in 
the complex s-plane) 



*Note that t < is necessary for both the s- and u-channel reaction to be physical. 



• right-hand cut (RHC): physical s-channel cut along s > s+, 

• nucleon pole: at s = m^ from the s-channel nucleon-exchange pole term l/(s — m^), 

• crossed cut: along s < s_ as combination of the n-channel cut s < S- and the t-channel cut 
s< -S_, 

• left-hand cut (LHC): collective name for all cuts in the unphysical region, i.e. for Re{s} < s+. 

In addition, the mapping between the complex s- and g^-planes involves a circular cut in the complex 
s-plane at |s| = S_ = S(+)S(_), where S(_(_)(g^) and S(_)((7^) are the two solutions 



S(±)(g2) = 2q^ + S+ ± 2V(g2 + m^)iq^ + M^) (2.15) 

for a given q^ (note the cut for —rm? < q^ < —M^) with S(_|_)(0) = s+ and S(_)(0) = s_. This circular 
cut becomes relevant once amplitudes are considered as functions of q^ rather than s, e.g. for the 
partial waves. The additional analytic structures of the s-channel partial- wave amplitudes due to the 
partial- wave projection are 

• kinematical cuts: for s < from terms depending on W = ^/s in the partial-wave projection 
formula, 

• short nucleon cut|j along T?_/m? < s < m? + 2M^ from evaluating the u-channel nucleon- 
exchange pole term l/{u{s, Zg) — m?) for Zg = ±1, 

• circular-cut contributions: from t-channel exchange of particles with mass mt > 2M,r) i-e. eval- 
uating l/{t{s, Zg) — ml) for Zg = ±1 and ml = tj^, 

• crossed-cut contributions for s < and singularities at s = 0: from partial-wave projection of 
the aforementioned u- and t-channel exchanges. 

Finally, we mention some kinematical points of specific interest (cf. e.g. |291I37| ): the Cheng-Dashen 
point at (s = u = m?^t = 2M^) = (z^ = 0,z^b = 0) is pivotal for vrA'" a-term physics, since the 
Born-term-subtracted amplitude D~^{y = 0, t = 2M^) = A+(i/ = 0, t = 2M^) - g^ /m is related to the 
a term by a low-energy theorem [10l[20l[38HlD] |j The subthreshold point at (s = it = sq, t = 0) = (z/ = 
Q,vb = ~-^w/(2'^)) serves as expansion point for the subthreshold expansion, while the (s-channel) 
threshold point (s = Sj^,t = 0,u = s_) = {u = Mt^,ub = — M^/(2m)) is relevant for the threshold 
expansion/parameters (e.g. scattering lengths). 

2.2 Isospin structure 

The most general Lorentz-invariant and parity-conserving T- matrix element for the process 7r"(g) -|- 
^{p) ~^ ■""((?') + ^{p') with isospin indices a and h is given in terms of Lorentz-invariant amplitudes 
A, B, and D according to 



T'ftis,t) = -{r^rn^+(s,t) + -[r" , t'^]T-{s , t) = d^^TUs^t) + ie,acT^Tj,{s,t) , 



1 1 

T},is,t) = Uf{p')!^A\s,t) + ^lB\s,t)Y,ip) = Uf{p')!^D'is,t) - IMi?^(s, t)|u,(p) , 
D^{s,t) = A^{s,t) + u{s,t)B^{s,t), I £{+,-}, (2.16) 



■'Actually, there are two short nucleon cuts as discussed in the appendix of [36]. The second one, however, is situated 
on an unphysical sheet. 

•^Note that since af ° = Zs(m^, 2M^) = -Af^/(4m^ - M^) ^ --5.56 x 10~^ is close to zero, the amplitudes at the CD 
point are dominated by the (s-channel) S'-wave. 



where we have introduced the isospin index / = +/— for the part that is even/odd under interchange 
of a and b. Furthermore, the ttN scattering ampHtudes A have definite crossing properties under 
interchange of s and u for fixed t, i.e. under change of sign of u, such that one can work with 
amphtudes 

' A{u,t) a A{u,t) = +A{-u,t) , 

^ a A{u,t) = -A{-u,t) , 



A{u'^,t) 



(2.17) 



which are even functions of v and thus free of kineniatical square root branch cuts in the complex 
i-plane originating from pt or qt- Exphcitly, the above amphtudes fulfill 



A^{v, t) = ±A^{-v, t) , B^{v, t) = ^B^{-v, t) . 



(2.18) 



The amplitudes of all ten vrA^ scattering reactions can be written in terms of only two independent 
matrix elements with total s-channel isospin index Is G {1/2,3/2}. In agreement with [29] (i.e. using 
the usual Condon-Shortley phase convention for the Clebsch-Gordan coefficients |34j . but the non- 
cyclic kinematical convention according to Sect. 12. l|) we assign the isospin-doublets of both the nucleons 
and antinucleons according to the fundamental representation of the Lie-algebra of SU{2) 
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(2.19) 



and the isospin-triplet of the pions according to 



TT 



|1,1), |7rO) = |l,0), K-) = |l,-l), (2.20) 

which leads to the following properties under charge conjugation C 

C\p) = \p), C\n) = -\n), C|7r±) = -Ivr^) , C\7t^ 

Thus, the relations between the spherical and the Cartesian components of the pion-multiplet are 

1 



(2.21) 



\w 



±\ 



=F-7|(|7ri)±i|7r2)) 



vr 



Fs; 



(2.22) 



By decomposing the initial and final isospin states of the vrA^ system into linear combinations of 
s-channel isospin eigenstates, e.g. 



|7r+p) 



3 3 
2' 2 



k p) 



3 1 
2'~2 



1 1 
2'~2 



|7r°n) 



1 



+ 



1 1 

2'~2, 
(2.23) 



we can readily obtain the relations between the vrA^ isospin amplitudes 



A+ = A{tt^p -^ TT^p) = A{Tr n 



IT n) 



A-^ -A- 



^3/2 



A- = A{7r-p -^ TT-p) = A{7r+n -^ n+n) = ^+ + ^^ = h2A^^^ + A^^^) , 

o 

^0 = Ain'p -> 7r°n) = ^(vr+n 
^(A^A) =^(vr°n- 



vrV) 


= -V2A- 


= ^(^1/2 ^3/2 
o 


vr^n) 


= A+ 






A+ + 2A- 


= ^1/2 . 



(2.24) 
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From these we can infer the so-called isospin triangle relation 

A+-A- = V2A0 , (2.25) 

and the relations for the isospin even/odd amplitudes with / = +/— and the amplitudes in the 
s-channel isospin basis Is G {1/2,3/2} can be summarized in matrix notation as 

A- J = ^"^ U3/2J > (^^3/2] = Csu i^j^-j , C,s = 3C.. = -i^^ _ J . (2.26) 

The s-channel isospin amplitudes with /^ G {1/2,3/2} and the corresponding u-channel isospin am- 
plitudes with /„ G {1/2 = N,3/2 = A} can be shown to obey the s -f-T- u crossing isospin relations 

U'^v " ^'" U^J ' U^J " '^'" U'^v ' ^- = chs = 3 (^ 2 ij ' (^•^'^) 

and combining this with (12.26P yields 

^+^ _ r M"^^ r -r r -U^ ^\ r -r-'-(^ -^\ (2 28I 

A- j — ^i^u I ^A I 1 '-'I'M — '-'us^su — o I _]^ ]^ I ' '-^"'^ ~ ^i^u — \^ 1 / ' l^'-^oJ 

For the t-channel reactions, the \NN) isospin states are superpositions of the states \It = 1, (-^4)3) 
and I It = 0,0) 

|np) = |l,l), |nn) = -^(|l,0) + |0,0)) , |pp) = -^(|1,0) - |0,0)) , |pn) = |l,-l), 

(2.29) 
from which we can deducqil 

\1,0) = ^{\nn) + \pp)), |0,0)=-^(|nn)-|pp)), (2.30) 

whereas the decomposition of the Ivrvr) isospin states reads 

Ivr+^O) = ^(|2,1) + |1,1» , Ivr+vr-) = i=|2,0) + -^|1,0) + -^|0,0) , 

Ivr-vrO) = ^(|2,-1) - |1,-1)) , kV) = y||2,0) - -^|0,0) . (2.31) 

By strictly using the non-cyclic kinematical convention together with the properties under charge 
conjugation p.2ip we can obtain the t-channel amplitudes from the s-channel ones via crossing 

A± = -A{pp^TT^7r^) , Ao = A{7r+n^7r'^p) = -A{np^TT+TT°)=A{np^'iT^TT+), (2.32) 

which together with the s-channel isospin relations (j2.24p on the one hand and the t-channel isospin 
decompositions above on the other hand yields the following relations for the reactions with a proton 



'^Note that (|2.29|) and (|2.30p are in perfect agreement with the usual Clebsch-Gordan coefficients [31], but differ 
from [25] wherein different conventions are used in these and corresponding equations. In particular, the analog of H2.33|) 
in [29] seems to (exceptionally) follow the cyclic kinematical convention. Nevertheless, all other relations, especially the 
crossing matrix (|2.34p and the important relations (|2.35[) . are identical. 
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as target particle 

A{pp -^ vr+vr- ) = -A+ = -A+ +A- = -^3/2 ^ _ 1^40 ^ 1^1 

A{pp -^ 7r-7r+) = -A- = -A+ -A~ = -^(2A^'^ + A^'^) = -^^° - \a^ , 

A{np^ 7r+7r'') = -Ao = V2A- = ^(^1/2 _ ^3/2) = 1=A\ 

A{pp ^ Tr\^) = 1{A+ + A-) = A+ =l{A^/^ + 2A^/^) = ^A^ . (2.33) 

2 3 ^/Q 

Thereby we can easily deduce the s -f-)- 1 crossing isospin relations 

(2.34) 
and the fact that A'^ and A~ have well-defined quantum number It = and /t = 1, respectively, 





^_y - C,i (^^1 ) , C,t = C.sCst = K^ J . (2.35) 

Since 

G|7r) = -|7r) ^ G|7r^) = Ivr^) , (2.36) 

the antinucleon-nucleon initial state in the reaction NN — )■ vrvr has to be an eigenstate of G-parity 
with eigenvalue +1, i.e. it can only couple to states with an even number of pions. The result for 
charge conjugation of an antifermion-fermion or antiboson-boson pair 

C\ff) = (-1)^+^1/7) , C\bb) = (-1)^166) , (2.37) 

yields 

G\NN) = {-iy+^'\NN) , (2.38) 

from which we can conclude that for reactions with a two-pion final state (i.e. G = +1) only the 
combinations (J even, It = 0) and (J odd. It = 1) are allowed. The same combinations arise from 
the symmetry properties of the symmetric isosinglet for It = and the antisymmetric isotriplet for 
li = 1 due to the fact that the exchange of two pions in an orbital state with total angular momentum 
J = L yields a factor of (— 1) . According to (J2.35P this leads to the following selection rules for the 
partial- wave decomposition of the t-channel amplitudes: the partial- wave expansion of the amplitudes 
A '~ or A^^^'^ contains only partial waves with even/odd J, respectively, and the transition 
between the two sets of amplitudes involves the isospin crossing coefficients cj with 

, ,„ if J is even , 
II J IS odd . 

2.3 Hyperbolic dispersion relations 



In [21] it was shown how to construct HDRs for the vrA^ scattering amplitudes, using hyperbolae in 
the Mandelstam plane of the form 

{s-a){u-a) = b, (2.40) 
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with hyperbola parameter b and asymptotes s = a and u = a. They obey the relation 

b 



and 



t 



s(t; a, 6) = - (S — t + Amu^t; a, 6)) 



s — a 



+ S — s — a , 



(2.41) 



4mi/(t;a,6) = v'(t - S + 2a)2 - 46 



u(t; a, 6) = - (S — t — 4?TT,i^(t; a, 6)) 



i(±)(z^;a,6) = S - 2a ± 2v^(2mi/)2 + 6 . 

In the following b is considered as a function of s and t for a given value of a, 

b{s, t; a) = {s — a)(S — s — t — a) , 



(2.42) 



(2.43) 



and hence for given s and a one considers a family of hyperbolae wherein all members are uniquely 
defined by t. Under the assumption that no subtractions are necessary (cf. Appendix [Dl), the HDRs 
for the vrA^ scattering Lorentz-invariant amplitudes can be written as 



oo 

A+(s,i;a) = - ds' 



1 1 
+ 



1 



A {s,t;a) 



vr 



ds' 



oo 

T . + / / /N 1 /■ , / Im^+(s',t') 

lmA+{s',t') + -J dt' ^,^^' \ (2.44) 



s — s s — u 



lmA-(s',t') + - I dt' ^ ' 



vr 



oo 

B+{s, t; a) = N+{s, t) + - f ds' 



u' t' -t 

oo 

I [ u Im B+{s',t' 



lmB+{s',t') + - I dt' - 

■K J V' 



t'-t 



B~ (s, i; a)=N'{s,t]a) + - / ds' 

vr 



1 1 

+ 






where we have defined the abbreviation 

u'{s',t') = u{s',t') 

and under the integrals one has to use 
b[s,t] a 



t {s , s,t] a) 



2s' + t' - S 
Am 



1 



(2.45) 



s — a 



+ S-s'-a, s'(t',s,t;a) = - S - t' + A/(t' - S + 2af - 46(s, t; a) 



(2.46) 



since the external kinematics {s,t,u) and the internal kinematics {s',t',u') are related by 

{s-a){u-a) = b= {s' -a){u' -a) , s + t + M = S = s' + t' + n'. (2.47) 

Only the amplitudes B^ contain the Born-term contributions N^ due to the nucleon poles given by 
(cf. [29] for N^) 



N+{s,t) = N+{s,t), 



N-{s,t;a) = N-{s,t) 



m^ — a 



N+{s,t)=g^ 
N-{s,t)=g^ 



rn^ — s rn^ — u 



1 1 

+ 



rn^ — s rn^ — u 



a V 
m u^ — u"^ 



(2.48) 
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where the usual pseudoscalar ttN couphng constant g and the pseudovector vrA^ couphng constant / 
are given bjo 

7^ = ^ - 13.7 . (2.49) 

47r M2 ^ ' 

In order to express the integrands in terms of the corresponding CMS scattering angles according to 

X{t\4) = X{s',t') Xe{A^,B^], (2.50) 



q'{s') = q{s') , 
p'tit') =Pt{t') = ip'.it') , q[it') = qtit') = iq'_{t') , (2.51) 



X{s\z',)-- 


--X{s',t') 


t'=t'{s',z'^) ' 


we define 




<{s\t')- 


= ^^(«''0 = 1+2^,2 


z[{s\t')-- 


= zt{s',t') 


mv' 

p'A' 



which leads to the relations 

t'{s',z'^) = -2c^\l-z'^), z'^{s',s,t;a) 

1 



2q 
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b{s,t;a) 

1^ + a-\ 



s'{t',z[) = -i^-t' + Ap[q'tz[), z'tit', s, t- a) = -^Vit' " ^ + 2a)2 - 46(s, t; a) . (2.52) 

Note that b is linearly related to z'^ for the s-channel, but only to z'f^ for the t-channel, which will have 
important consequences in Appendix [Cl where it will be shown that the HDRs ()2.44p incorporate 
contributions from the direct as well as from the crossed channels, but not from double-spectral 
regions, provided the parameters are chosen appropriately. Furthermore, one can check explicitly that 
A~^ and B~ are indeed functions of i''^, while A~ and -B"*" are proportional to z^. Since moreover 
4r7T,z/' = 4:p[q'i-z'f- = s' — u' is always real, one may also write the above HDRs (j2.44p in terms of reduced 
amplitudes A^ /u and B^ /u, respectively. This fact will be used in Sect. [Hand Appendix ID. 3[ 

In contrast, for usual fixed-t dispersion relations external and internal kinematics are related by 

t = t' , s + t + u = T. = s' + t' + u'. (2.53) 

It is remarkable that the HDRs have the simple form of (12.44p and (I2.48p . which by neglecting the 
terms depending on a (or equivalently for \a\ — )• c«) reduce to fixed-t dispersion relations, provided, 
however, that the f-channel integrals are discarded. Moreover, the hyperbolae then reduce to fixed-t 
lines, and thus we will refer to the limit \a\ — )• cxd as "fixed-t limit" in the following |£| 

3 Roy— Steiner system for pion— nucleon scattering 

In this section, we first collect the results for the partial-wave hyperbolic dispersion relations (PWH- 
DRs) that follow from the HDRs p.44p via partial-wave expansion in and projection onto both s- 
and t-channel partial waves as explained in detail in Appendices |A] and [Bl in order to state the closed 
system of RS equations for vrA^ scattering. Then, we elaborate on the corresponding partial-wave 
unitarity relations for the s- and especially the i-channel. Finally, we use the threshold behavior of 
the t-channel partial waves f±{t) in order to cast the t-channel part of the RS system in the form of 
a MO problem, whose solution will be the subject of Sect. [5l 



Note that [29] quotes a value of 14.28 based on 21] ■ For niore information on conventions as well as the current 
value see [3lll|4l|43] . 

^As explained in Appendix lC.2l only the limit a — >■ — oo is compatible with range-of-convergence considerations. 
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3.1 Partial- wave hyperbolic dispersion relations 

The s-channel partial- wave amplitudes are conventionally denoted by fi_^{W) with isospin (i.e. cross- 
ing) index / G {+> ~} and total angular momentum j = / it 1/2 = ^ib, where the orbital angular 
momentum can take the values / > for j = 1+ and I > 1 for j = I—. Using a shorthand notation for 
the Zs-projections of the invariant amplitudes 



1 



Xl{s)= fdzsPi{zs)X'{s,t) , , ,, , for Xe{A,B}, (3.1) 

J t=t{s,Zs) = -2q''{l~Zs) 

-1 

the well-known s-channel partial-wave projection formula reads [44] 

fl±(W) = j^{iE + m)[Ai{s) + {W-m)B{{s)]+{E-m)[-Al^,is) + {W + m)Bl^,{s)]}. 

(3.2) 

By construction, the fi^{W) obey the MacDowell symmetry relation [35] in the complex VF-plane 

f!+m = -fli+iyi-W) V / > , (3.3) 

due to which only half of the complex VF-plane is actually needed. Alternatively, this relation can 
be used the other way around to derive the partial waves with j = I— from the ones with j = 
l-\-. Expanding the absorptive parts of the HDRs (|2.44p into s-channel and t-channel partial waves, 
respectively, and subsequently projecting the full HDRs onto the s-channel partial waves fl_^{W) yields 
the s-channel PWHDRs of [II] 

1 1 
fl^{W) = Nl_,{W) + - / dW Y.[Ki,{W,W')lTafl,^{W') + KUw,-W')\rafl,_^^^_{W')} 

oo 

+ ^/d*' j;{G,j(W^,t')Im/4^(0 + ^;j(W^,0lm/:'(i')} 

f J 

= -/(U)-(-^) V/>0, (3.4) 

which constitutes the s-channel part of the full RS system. Here, N^jjy^^ represent the contributions 
due to the nucleon pole terms in the amplitudes B^ as given in (j2.44p . Each s-channel partial 
wave /;{|_(VF) is coupled to the absorptive parts of all other s-channel partial waves via the kernels 
Kii,(W, W), which contain the usual Cauchy kernel responsible for the physical cut and an analytically 
known remainder (denoted by dots below) containing only left-hand cut contributions 

Kl,{W,W') = ,^^^^P^ + ... V/,r>0, (3.5) 

as well as to the absorptive parts of the t-channel partial waves f±{t) via the kernels Gij{W,t') and 
Hij(W,t'), where the lower index it denotes parallel(-l-) or antiparallel(— ) antinucleon-nucleon helic- 
ities and the total (i-channel) angular momentum J can take the values J > or J > 1, respectively. 
Due to Bose statistics (i.e. crossing symmetry), the summations over J in (j3.4p run over even/odd 
values of J for the crossing even/odd partial waves (upper index / = +/—), respectively, as explained 
in Sect. 12. 2[ For the sake of completeness and convenience, in Appendix [X] the different contribu- 
tions to (13. 4p will be discussed along the lines of |211I461[48] (correcting several typographical errors, 
adjusting the conventions, and partially extending the presentation therein at the same time). 
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For the t-channel partial-wave projection, by virtue of s -^ u crossing symmetry it is possible to 
use only half the interval in the cosine of the t-channel CMS scattering angle and thus the projection 
can be written as HS] 



p1 



Hit) 



iptqtY' 
1 



m 



jZtB\s,t) 



1 v/J(J + l) 1 
47r 2J+1 {ptqtY-^ 



dzt 



s=s{t,zt) (ptqt)-^ 
Pj-i{zt)-Pj+iizt)\B^{s,t 



s=s{t,zt) 



s=s(t,zt) 



V J>0, 



V J > 1 , (3.6) 



where again / = +/— if J is even/odd, such that the integrands are always functions of the squared 
angle zf. These formulae are valid literally only for t >t]\i, but can actually be used for all kinematical 
cases, cf. the discussion following (|A.73p . For a closed system of RS equations we need to derive the 
analog of (j3.4p for the t-channel partial waves /±(t), cf. |1 H I19 1 [3T]. The result takes the form 

1 /■ 

fi{t) = Ni{t) + - / dW'Y,{Gjiit,W')lmf!^iW') + Gjiit,-W')lmf(,^,y{W')] 

oo 

+ - /di' j; {i^}j.(t,Oim/:^'(t') + ^L,(t,t')Im/:^'(t')} V J > , 
^ t J' 

1 /■ 

fl{t) = N:lit) + - / dW'^{Hji{t,W')lmfl^{W') + Hji{t,-W')lmf(^^,-^_{W')] 

oo 

+ - I dt'^k]j,{t,t')\Tnf':^{t') VJ>1, (3.7) 

^ t J' 

where again / = +/— if J is even/odd and the sum over J' runs over even/odd values of J' if J is 
even/odd (cf. Sect. 12.2]) . As for the s-channel case, the kernels for the corresponding t-channel partial 
waves can be split into the Cauchy kernel and well-defined remainders 



K]At,t') 



t'-t 



+ ... V J, J' > , 



K^jj,{t,t') 



t'-t 



+ . . . V J, J' > 1 , 



(3.8) 



but, in contrast to the s-channel case, only higher t-channel partial waves can couple to lower ones, 
since Kj'j! (t, t') = for all J' < J, which will be a key ingredient in reducing the t-channel part l\3.7\i 
of the RS system to a MO problem in Sect. 13.31 The technical details of the derivation of the different 
contributions to (j3.7p are relegated to Appendix [Bl 

There are three aspects of convergence in the RS system of PWHDRs constructed in Appendices El 
and [HI first, the question of convergence of the integrals in the high-energy regime is linked to the num- 
ber of necessary subtractions of the dispersion relations, which will be discussed in Sect. |H Moreover, 
for the full system of RS equations to be valid, the convergence of both the partial-wave expansion of 
the imaginary parts inside the integrals and the s- and t-channel partial- wave projection of the full 
HDR equations needs to be shown. Analyzing these two constraints yields the ranges of convergence 
in s and t for (j3.4p and (13. 7p . respectively. As explained in detail in Appendix [Cj the hyperbola 
parameter a can actually be tuned in order to obtain the largest possible domain of validity. For the 
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s-channel part (j3.4p of the RS system the combined analysis of s- and i-channel constraints leads to 
an optimal value of a and a corresponding range of convergence in s of (cf. Appendix IC.3P 

a = -23.19 M^ ^ sG [s+ = (m + M^)2, 97.30 M2] ^ VF e [1^+ = 1.08 GeV, 1.38 GeV] , 

(3.9) 
where s+ = 59.64 M^, while for the t-channel part (j3.7p we find (cf. Appendix IC.4p 

a = -2.71M^ ^ t e [U = 4M^, 205 A5M^] ^ Vt G [V*^ = 0.28 GeV, 2.00 GeV] . (3.10) 

Note that different choices of a for the s- and i-channel partial- wave projections are perfectly justified, 
as we may start from different sets of HDRs. However, the choice of a is not only crucial for the ranges 
of convergence, but also influences the high-energy behavior of the imaginary parts, whose estimation 
via Regge asymptotics is discussed in Appendix [Dj For this purpose one splits the corresponding 
integration ranges s+ < s' < oo and tT^ < t' < co of the HDRs (I2.44p at some appropriate values 
Sa = W^ and t^, respectively, in order to describe the asymptotic s- and t-channel contributions 
to the invariant amplitudes in terms of Regge amplitudes. The remaining non-asymptotic parts are 
then given by the corresponding integrals over s+ < s' < Sa and t^^ < t' < ta_, respectively, plus the 
nucleon pole terms N^(s,t) for the amplitudes B^{s,t). However, eventually the high-energy region 
is of only little practical relevance, in particular if subtractions are performed in order to suppress the 
dependence on higher energies (cf. Sect. |3|). 

In order to use partial-wave unitarity relations that are diagonal in the s-channel partial waves, we 
have to work in the s-channel isospin basis Ig G {1/2,3/2} rather than in the isospin even/odd basis 
/ = +/— (as will be explained in Sect. 13. 2p . and therefore in analogy to (|2.26p we define 

(^3/2) = Csu (^!) , (^!) = C,s (J/s) , for X G {fi±,Ni±,Ku'} , (3.11) 

and the abbreviation 

^1/2+3/2 ^^^ ^,^ ^ _^i/2 ^^^ ^,^ ^ ^3/2 ^^^ ^,^ ^ 2i^+ (W, W) + K~ {W, W) . (3.12) 

The full closed RS system of PWDRs for both s- and t-channel partial waves in the corresponding 
isospin bases Is G {1/2, 3/2} and It G {0, 1} that follows from rewriting (j3.4p and ()3.7p reada^i 

-If 00 ^ 

f!i\w) = nII\w) + - / dW'Y,^{K'''iW,W')lmf;/_^{W') + 2Kf/,\w,W')lmff/_^{W') 

+ Ky,\w, -W) lu.fll,l^^_{W') + 2K'J\w, -W) lu.flll^^_{W')] 

+ - I'dt'Y, ^^ ~ ^'^^'"^ {Gu{W,t')lmfiit') + HuiW,t')lmf:l{t')] , 
^ / J=o 

+ KfJ,\w, -W) Im flll^^^ {W) + KII^^'''\W, -W) Im /f/J^)_ {W')} 



■"All sums run over both even and odd values, and the formulae for the fn^^^ are given explicitly for convenience. 
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+ - / dt'Y,{-^)'{Gij{W,t')ln^fi{t') + Hu{W,t'n^f-{t')} , 
^ / j=o 

1 /■ oo ^ 

+ kI/\-W, -W) Im /(yj^)_ (1^0 + 2Kf/\-W, -W) Im /f/^,)_ (W)} 
^ / j=o 

■'TT 

1 /■ oo ^ 

-- / ^W'Y,^{KT{-W,W')ln^fl/^{W') + Kll^^'''\-W,W')lmfl'^{W') 

+ i^,f (-ly, -M^O Im /(V' 1) Ji^O + ^//'+'^'(-^, -W^O Im /f/^i) JT^')} 

--/dt'j;(-l)^{Qj(-H^,Olm/:^(t') + ^L/(-I^,Olm/i(t')}, (3-13) 

^ / j=o 

together with 

oo 

■^ 7 n 



«=o 



Im/,^^^(iy') 



oo 

+ - / d*' E ^V {^)j'(i,i')Im/:^'(t')+i^3j'(t,0lm/i'(t')} VJ>0, 



J' = J 



1 /" C)0 I 



Im//f (M^O + i±^LlZT^ ^3/2,,,,,, 



-Im/;_^^(M^' 



/ J'— ] 



'^ , l + 3(-l)\^.3/2 



I- /('/+!)- (H^')]} 



l + (-l) 



-i^3^,(t,t')Im/i(tO VJ>1. 



(3.14) 



Note that in the above i-channel part ()3.14p the sums over J' are hmited to J' > J due to ()B.46p . 
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3.2 Partial-wave unitarity relations 

From the unitarity of the 5-matrix S = 1 + iT one can easily obtain the general unitarity relation by 
taking matrix elements and inserting a complete set of intermediate states 

(/|r|i)-(/|rt|z) = z^|dn(f;(/|rt|i)(i|r|i), (3.15) 



{j} 



r(i) 



where dUn' denotes the nj-particle Lorentz-invariant phase space (LIPS) for intermediate state j, 
which in the case of rij identical intermediate particles implicitly includes an additional symmetry 
factor 1/Sn- = l/rijl in order to avoid multiple counting in the phase space integral. Imposing overall 
4-momentum conservation 6^'^'{T,pf — T,pi) and using time-reversal invariance of the strong interactions 
immediately yields the generalized optical theorem for the dimensionless invariant amplitudes Tfi 

ImTf, = IY1 /dn(^)(2vr)45W(Sp, - Sp,)r;^.r,, . (3.16) 

{J} ■' 

Under the additional assumption of hermitian analyticity of the S-matrix (i.e. the amplitudes Tfi obey 
the Schwarz reflection principle T%{s) = Tfi{s*) and are real on part of the real axis) it follows 

V)iscTfi{s) = lim [Tfi{s + ie) - Tfi{s - ie)] = 2i limlmr/i(s + ie) , (3.17) 

for the physical limit corresponding to the s-channel process, and hence ()3.16p may also be proven in 
the framework of perturbation theory to all orders. By normalizing the 4-momentum states according 
to {p'\p) = 2Ep(27:)^6^^' {p' — p) for both bosons and fermions, for generic two-by-two scattering 
ab -^ cd with one particular intermediate 2-particle state j = jiJ2 (with CMS 3-momentum modulus 
Pj) and after partial integration of the 2-particle LIPS the optical theorem (13.16P takes the form 



_^j_2p^rdn^ _ A|^^ 

5O) levr yfs J ATI ^^ '' ' ^^ V 4s 



I-^/^ = TTTYTT-^ ^^/.^.^ ' ^. = \K- ' (3-18) 



leading to the usual form of the differential cross section (with pf and pi in analogy to pj ) 



dafi _ PfPi dafi _ pj_ 
dQ vr dt Pi 



'^fi 



SvrVi 



(3.19) 



A partial-wave decomposition of the invariant amplitudes Tfi allows for a reduction of the unitarity 
constraint (j3.18p to unitarity relations for each partial wave separately. In the presence of spin the 
T operator for two-by-two scattering can be diagonalized by using the eigenstates of total angular 
momentum J as basis, which can be achieved most easily in the CMS via the helicity formalism ^50j . 
With \p denoting the helicity of the corresponding particle, one can take the T-matrix elements in the 
basis of single particle momenta and helicities and by applying the respective phase space integration 
in the CMS, the corresponding invariant helicity amplitudes Tr?' '*' "' * can be written in terms of 
states of relative motion for both incoming and outgoing particle pairs. Thereby, the differential cross 
section for a reaction with a given set of helicities can be derived in full analogy to (j3.19p . With the 
usual angular conventions of [291I50J and the azimuthal angle 99 set to zero, the partial-wave expansion 
of these helicity amplitudes in the helicity basis then reads 

T}r''-''-'\s,t) = 7;S^16vr^(2J + l)r/^,,^^,^,,^(sX_,^,,^_,^(0) , (3.20) 

J 
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where d^^; (0) are the Wigner d-functionq^ and the sum runs over integer /half-integer values of J for 
an even/odd number of half-integer spins present in the initial or final state. In the case of spinless 
particles with (ioo(^) = Pj {cos 6) and J = I the expansion simplifies to 



Tfi{s,t) = ^/SjSd67T}_^{2J + l)Tf^{s)Pj{cos9) . (3.21) 

j=o 

Note that we have added here explicit symmetry factors Si and Sf to the partial- wave expansion of }50j 
in order to take care of identical particles in the initial and final state in a symmetric fashion. This 
normalization reproduces the standard normalizations for spinless processes as well as for vrA^ — )• vrA^, 
and furthermore ensures that no symmetry factors occur in the elastic unitarity relations for the 
partial waves, since they always cancel with the symmetry factor implicitly included in the LIPS 
(cf. (j3.18p ). We will explicitly demonstrate the effect of this convention for the symmetry factors on 
the extended unitarity relation for tttt — )• NN partial waves by considering tttt — )• tttt with KK and 
NN intermediate states below. 

Due to the invariance of strong interactions under time reversal and parity, the helicity partial 
waves obey the symmetry properties 

^Aa,Ab;A„Ad(^) = 'PXc,Xd;K,Xb(^) = '^-A^ -Ad;-A„ -AftC^) • (3-22) 

If the particles are spinless or if the matrix T'^{s) in helicity space is diagonal in some appropriate basis 
(as it is e.g. for ttN — )• vrA^ in the s-channel isospin basis Is G {1/2,3/2}), the unitarity relation ()3.18p 
for partial waves of generic elastic scattering ab — )• ab (i.e. f = j = i) reads 

Imr/,(s) = ^|r/,(s)p, (3.23) 



which is solved by a parameterization of T'L{s) via the real phase shift S'f^{s) 



L fj\oj \ idi tllb icai piidoc omit ^fj\ 



T/,{s) = ^smSJ,{s)e^'U^\ (3.24) 

where (|3.23p and (|3.24p are valid for each diagonal element T^^{s) of T'^{s). For s above the lowest 
inelastic threshold Sinei these equations have to be modified by introducing real inelasticities < 
Vfii^) ^ 1 according to 

with r]j-{s) < 1 for s > Sjnei due to additional intermediate states contributing in (j3.15p . These partial 
waves are then related to the diagonal elements of the corresponding S-matrix via 

SUs) = vUs)e'''f^^'^ = 1 + i^Tf,{s) . (3.26) 

After these general remarks, we now turn to vrA^ scattering: the reduced s-channel partial-wave am- 
plitudes fi^{W) in the s-channel isospin basis Ig G {1/2,3/2} are conventionally normalized according 
to (cf. (Km and e.g. [201E9]) 

fi. (^) ^ 1 MWLv^ZLiV^ = l r,/|(Ty)e^<^(^)-l h'<h/.„, sin SJ^JW) ^.sfuw) .^_^^. 

^^ q 2i q 2i q ' ' 



^^A comprehensive review on Wigner functions, in particular a comparison of different angular conventions used in 
the literature, is given in |51| . 
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where for the elastic form we have used the fact that the lowest inelastic intermediate state is ttttN 
and thus r]i^{W) = 1 below the inelastic threshold W-mei = W+ + Mj^. The s-channel partial- wave 
unitarity relation corresponding to the normalization given above reads 



im//|(Ty) = q\ft{w)fe{w - w+) + 



i-{vii{w)y 

4q 



0{W-W,^ei) , 



(3.28) 



leading to the branch cut for W > W+. 

For the (necessarily inelastic) t-channel partial-wave unitarity relations one needs the dimensionless 
partial- wave amplitudes tj{t) of elastic mr scattering. They are conventionally defined from the 
dimensionless isospin amplitudes of vrvr — )• vrvr via (with t-channel isospin It S {0, 1,2}, total angular 
momentum J = I, and symmetry factors y^SfSi = 2 for identical pions, cf. (j3.2ip and [T2ll20] ) 



r^'(s, t) = 32tt ^(2 J + i)tj' (t)Pj(cos e^ 

that are normalized according to 



J=o 



da: 



TTTT — ?7r7r 



dn 



T'^{s,t) 



8iry/i 



The corresponding elastic unitarity relation then takes the form 

Im t^J (t) = aj\t^J {t)\'^e{t- U) , 
and hence the partial waves can be parameterized as 



(7t 



2gt 



T 



T 1 \Sj{t)] -1 



1 r]]'{t)e^'^'}^''> - 1 v'l(^=i sinS^Jit) ^gh 



(J, 



2i 



o-; 



2i 



e^^yw 



(7, 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



The reduced t-channel vrA^ partial- wave amplitudes f^{t) are related to vrA^ helicity amplitudes 
F)^^{s,t) and dimensionless partial waves F^{t) via (cf. 



F++{s,t)= F__(s,t) 
F+_(s,i) = -F_+(s,t) 



qt 

47rVt 

qt 



^(2J + l)F_/(t)Pj(cos0i) 



j=o 

oo 



pt Vt 



^ ^_^ Fi(t)sing,P}(cosg, 
=1 V'^V-' + iJ 



J- 



Pt 



and they are normalized according to 

2 



^'^■K-K^NN Pt 



dO 



qt 



E 

A, A 



F-xx{s.t) 



SnVt 



2pi 
9i 



F++(s,t) 



SnVt 



+ 



F+-{s,t) 



Svr^/t 



4Pi df^iVTV^^TT 



9i 



dO 



(3.33) 



(3.34) 



The general formulae (j3.33p and ()3.34p are also valid for isospin even/odd parts F^Js, t) with crossing 
index / = +/— and J even/odd, accordingly. Note that when referring to the t-channel isospin basis 
It e {0, 1} as in the following, the isospin crossing coefficients cj of (|2.39p need to be included. In 
general, the t-channel partial waves may be parameterized as 



flit) = |/i(t)|e*'^'W = Re/i(t) + ilm/i(t) . 



(3.35) 
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By considering only vrvr intermediate states in the region t < [AMj^)'^ (which is elastic with respect to 
TTTT scattering, but unphysical with respect to the ttN t-channel) in the general unitarity relation (J3.15P 
for NN -^ TTTT, the f±{t) can be shown to obey the "elastic" f-channel unitarity relation 



Imflit) = a^{t'j{t)Yfl{t) d{t -U) V t G [U, 16M2) 



(3.36) 



(where the coefficients cj cancel), which leads to the branch cut for i > i^. Since the imaginary part 
Im/^(t) itself must be real, from (j3.36p together with (j3.32p and (j3.35p one can immediately infer 



fi 



it) 






it) 



e^'j 



it) 



V t G [t^, 16M^ 



(3.37) 



i.e. the phases of the t-channel partial waves f±{t) are given by the phases of the vrvr partial waves tj{t) 
modulo IT (by convention we choose the phases to coincide exactly), which is also known as Watson's 
final state interaction theorem [52] . It is common practice to assume that the contributions due to 
47r and other intermediate states can safely be ignored for t < 40M^ w 0.78 GeV^ (see e.g. p9lH0] ). 
However, as demonstrated in [53] in the context of the scalar pion form factor, this is certainly only 
true in the 5-wave below the threshold tx = 4M|^ ~ 0.97 GeV for the production of KK intermediate 
states, while in the P-wave inelasticities effectively start to set in around the ttoj threshold at 0.85 GeV^. 
It is crucial to note that (j3.36p is invariant under rescaling of f±{t) with real factors, whereas elastic 
unitarity relations as (j3.28p (for W < VFinei) and (j3.3ip are always nonlinear in the corresponding 
partial wave. Hence, fixing the normalization of all different partial waves that are needed in extended 
i-channel unitarity relations (i.e. allowing for additional intermediate states) in a consistent manner 
can only be done resorting to the corresponding elastic reactions, as we will now demonstrate for a 
system of coupled-channel equations with vr, K, and N degrees of freedom. Writing Tn = Tt^t^^t^t^ , 
Ti2 = Tf^j^^^^, Ti3 = Tjqjq^^^ etc. for the T-matrix elements and using the invariance of strong 
interactions under time reversal, the general unitarity relation reads in terms of matrix elements 



S}jSji 



5 



fi ' 



5 



fi - "fi 



5fi + iTfi = 6if + iTi 



fi - "if 



'if 



5, 



if 



(3.38) 



In particular, one can read off the extended elastic unitarity relation for vrvr — t- vrvr and the extended 
unitarity relation for NN — )■ vrvr with tttt, KK, and NN intermediate states 

<5ll = 1 = IS*!!! -h 15*121 + 1 5*13 1 , Si3 = = SliSi3 + SI2S23 + Sl^S33 , (3.39) 

and thus, by dropping the NN intermediate states in the second relation (since we are finally interested 
in the extended t-channel unitarity relation of vr A^ scattering in the region below the NN threshold) , 
we obtain 

2Imrn = \Tnf + \Tu\^ + \Tnf , 2lmTn = T^iTis + T^2T23 ■ (3.40) 

Introducing now the reduced t-channel partial waves gj{t) of ttK scattering (with isospin It = 0/1 
corresponding to J = I even/odd due to Bose symmetry in vrvr, symmetry factors x/SjSi = y/2 and 
the partial waves defined from dimensionless isospin amplitudes, cf. (I3.2ip and [19J) 



G'^is,t) = 16^V2^(2J + l)(A;t(7t)^5j (i)^j(cos0r'^) 



h 



j=o 



4 ^ 2 * 



(3.41) 



the first relation of (|3.40p may be decomposed into partial waves, and performing the angular inte- 
grations of the phase space integrals leads to the partial- wave unitarity relation for vrvr scattering with 
vrvr, KK, and NN intermediate states 



Im t^'(t) = a^ e]{t) 9{t-U) + {ktqtraf' g'j{t) e{t-tK) + 



16g2 c2 



N 



Fi{t) + Fi{t) 



^{t-tN) . 
(3.42) 
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For t < tx (or if Ij + J equals an odd number) this reproduces the elastic unitarity relation for vrvr 
scattering (|3.3ip . which corresponds to the relation (cf. (|3.32p ) 



[s'/it] 



TTTT— ^TTTT 



1+A'j{t)9{t-U) 



Comparing (j3.42p with the elastic unitarity relation for the partial waves (cf. (|3.39|) ) 



[s'jm 



TTTT— >-7r7r 



+ 



[s'/m 



TVK^^KK 



+ 2 



[si{t)t 



■KTT^NN 



+ 



[smt 



mr-^NN 



(3.43) 



(3.44) 



for the cases t < t]\f and t > t^ successively then allows to fix the normalization of the partial-wave 
S'-matrix elements of the inelastic channels (both in the natural f-channel isospin basis /( G {0, 1}) 



[snt)] 



Aihqt 



^J+i 



rnr^KK 



Vt 



g'jit)0{t-tK), [siit)Y; 



niT^NN 



cjV2 V It 

(3.45) 



These S-matrix elements indeed reproduce the correctly normalized differential cross sections 



da 



■K-K^fKK 



drj 



h 
It 



G'^is,t) 



SttV^ 



(3.46) 



and (j3.34p . respectively. Furthermore, from the unitarity bound of the t-channel partial- wave S matrix 
of 7tN scattering only (cf. (j3.44p ) together with its explicit form (j3.45p and the relations (j3.33p to the 
corresponding partial waves /^ we can deduce that the partial waves fall off asymptotically at least 
as fast as (cf. [29] ) 

fiit) ~ t-^+5 , f:[{t) ^t-\ for t ^ oo , (3.47) 

i.e. f^{t) —7- for t — t- oo by unitarity at least for all J > ; this asymptotic vanishing is usually 
assumed to hold for the 5- wave as well. By virtue of similar considerations, the normalization of the 
remaining partial waves in the second relation of (I3.40p can be fixed. We may introduce the reduced 
i-channel partial waves /i±(i) of KN scattering in analogy to the t:N case via dimensionless helicity 
amplitudes (cf. (13.33P and [54] ) 

47rVt 



H++{s,t) 



H+-{s,t) 



h 

AlTy/t 



^(2J + l)i7^(t)Pj(cos0/ 
2J+1 



KN\ 



J=0 

oo 



E 



',^/j{j+T: 



Hi{t)sin9rP'Acos9r) 



Hi{t) 
Hi{t) 



^\p,ktY ^hi{t) 



Pt 

h 
Pt 



ipthYhiit) , (3.48) 



where it is important to note that, in contrast to vrA^ scattering, also the combinations It = with 
odd J and It = 1 with even J are allowed due to lack of Bose symmetry in KK. In order not to bloat 
the notation, we refrain from using an additional index for It, and in the following e.g. /i_|_~'^™" 

1 4- u J J- J , (J=even/odd,/t=0/l) j.- i j j. t (J=odd/even,/t=0/l) - , 

IS always to be understood as n^. , respectively, and not n,^ , smce only 

the former can couple to the t-channel process vrvr — t- NN. The normalization is fixed by 



da 



KK-^NN 

An 



Pt_sr^ 
kt4^ 



H-xxis,t) 



8vrVt 




H++is,t) 



SttVI 



+ 



H+-{s,t) 



SvrVt 



4p? da^^^j^ji ^ ^^^^g^ 



kl 



dO 



^^Note that our symmetric normalization of the hehcity partial waves p.20p together with H3.33[) and (|3.34|l leads 
to an additional factor of 1/v^ to [^^(t)] * ^ in comparison with [291149) . where one should read in addition 



[5iwL.- 



Cj 



\smi- 
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so that the dimensionless partial- wave ainphtudes Hj^{t) are related to the diagonal elements of the 
corresponding S'-matrix according to 

HmiR^NN = ;iv # ^i(^) ^(* - *^) ' cr = \^J- (3.50) 

Cj \j l^t ^ 

Plugging the partial-wave matrix elements into the partial-wave projection of the second relation 
of ()3.40p for either parallel or antiparallel antinucleon-nucleon helicities yields the extended i-channel 
unitarity relation for the t-channel partial waves f±{t) (extending (J3.36P for KK intermediate states) 

Im/i(t) = ar(t^'(i))*/±(t)^(i-t.) +2cJ^/2A:2^af (5;^*(t))*/i;i(t)0(t-ti,) . (3.51) 

Finally, we can use p.44p to derive the inelasticities rij (t) of the vrvr scattering amplitude that are 
consistent with (|3.5ip . Below the NN threshold, inserting (|3.32p and (|3.45p into (|3.44p leads to 



r^'j{t) = ^Jl- Aa-^a^{hqtYJ\g]^{t)Y e{t-tK) ■ (3.52) 

3.3 From Roy— Steiner equations to the Muskhelishvili— Omnes problem 

3.3.1 Threshold behavior of the t-channel partial waves 

The asymptotic behavior of f^{t) for p^ — )■ and qt ^ {) (which is equivalent to t — t- t^v = Am? 
and t — )• tjr = 4M^, respectively) can be derived directly from the partial- wave projection ()3.6p . Since 
A\t, zt) and B^{t, zt) have definite symmetry properties under s ■^ u and since s — u = Amu = AptQtZt, 
we can write down the expansions 

A'{t,zt) = Y^iPtqtY' Pj'{zt)ar{i) , B^t.zt) = Y.{ptqtY' Pr{zt)hj,{t) , (3.53) 

J' J' 

where only even/odd values of J' contribute according to the symmetry properties of A^ and B^ (i.e. 
even J' for A~^ , B~ and odd J' for A~ , B'^). Let us first consider the limit pt — )• 0, i.e. the behavior 
of f±{t) at the f-channel threshold t]\f. As far as the leading asymptotic behavior is concerned, 
the functions aji[t) and hji{t) can be evaluated at t = tjy and will thus be considered as constant 
coefficients in the following. Inserting these expansions into (j3.6p (where J even/odd corresponds to 
I = +/-), we find for J = that 

fl{t^tN) = 0{pl) (3.54) 

at the physical threshold, while for J > 1 we obtain 

1 

fi{t ^ tiv) = ^^ j dz,PM)z,Pj^,{z,) + 0(p?) = ^^^ + 0{pI) , (3.55) 

-1 

1 



-1 

such that 

fiit -^ In) = 0{l) , fi{t ^ In) = 0{l) , V J > 1 . (3.56) 

However, the linear combination 



^\t) = ^\jj^f-{t)-fi{t) VJ>1 (3.57) 
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vanishes at threshold (cf. [29 1 I55 1 I56] ) 



V J> 1 . 



(3.58) 



The same reasoning may be apphed to the hmit qt — )• as weh, but as A^ contributes at the same 
order as B^ in the expansion of f:l{t), no relation between the threshold values of different amplitudes 
may be inferred for g^ — )■ 



fi{t^t^) = 0{l) VJ>0, f:l{t^t^) = 0{l) VJ>1. 



(3.59) 



In fact, the properties of f^{t) at the t-channel threshold are crucial to ensure convergence in the 
RS equations. From the partial- wave expansion (]A.73p we can easily derive the leading contributions 
to the invariant amplitudes (given explicitly for J < 2) 



A+{iy,t) 



Att 



— + —m V 

PI 2 



'' ''^ + UfUt) + 



PI 



^-(^'^)-3^.^ + ..., 



B+{iy,t) 15 ^ 



Att 

B-{u,t) 

4tt 



Pt 



(3.60) 



demonstrating how the threshold behavior (j3.54p and (j3.58p ensures that the partial-wave expansion 
does not introduce spurious kinematical poles at pt ^ into the expansion of the invariant amplitudes 
and thereby into the HDRs (I2.44p . To illustrate the consequences of this point, we briefly comment 
on the several places in our RS system (j3.4p and (j3.7p where the threshold behavior of f±{t) features: 



1. Although Gij{W,t') and Hij{W,t') diverge as p[ ^ for t' — > t^ according to (|A.92p . the relation 

Hij{W,t'),t' = tN 



Res 



-m 



J 



J+1 



Res 



Gij{W,t'),t' = tM 



(3.61) 



together with ()3.54p and (13.580 ensures that the corresponding integrals in (13. 4|) are well defined. 
We have checked that the explicit expressions in ()A.89P fulfill this equation. 

2. The pf"^ divergence ([B^B of k]j,{t,t') and Kjj,{t,t') for t' -^ tn cancels in ([321) provided 
that 



Res 



KjJ'i't,'t'),t' = tjy 



-m 



J' 



Res 



^jj'i^^t'),t = In 



J' + l 
This relation can easily be verified for the kernels given in ()B.48p and ()B.50p . cf. ()B.49p . 



(3.62) 



3. Based on the asymptotic forms (IB.SP of the pole-term projections Nil{t), one may check their 
threshold behavior to be analogous to (j3.54p and (|3.58p . Note that in this special case the 
relations hold for g^ — )• as well, since A does not contribute to the pole terms: 



Nliptqt ^ 0) = Oip^^q^) , m 



J 



j^-^Ni {ptqt ^Q)-Ni {ptqt ^ 0) = Oipfqt) V J > 1 . 

(3.63) 



3.3.2 Muskhelishvili Omnes problem for the f-channel partial v^aves 



Using the properties of the kernel functions for i-channel exchange as given in Appendix IB . 31 together 
with the threshold behavior of the partial waves as discussed in Sect. 13. 371] we can rewrite the i-channel 



25 



part (j3.7p of the (unsubtracted) RS system as 

oo oo 

n,^ n.x 1 /" ,Im/9(t') 1 f ,lm/0(t') 






/^">--i'"^^^' '"^^'1-t'"'^'"^ ^^^ -.,, 



C77 Cti 



/_^(t) = A^(t) + i|dt'^^^^ VJ>1, (3.64) 

where we have defined the abbreviations 

^i{t) = Ni{t) + Ki{t), 

A:[(t) = - / dTy'^{Gjz(t,iyOlm/4(VF') + Gjz(t,-H^')Im/{+i)-(W^')} 

+ ^/d*' E ^V {^}j'(i,0lm/^'(t')+^3j'(t,0lm/:''(0} VJ>0, 

/ J'=J+2 

A:{(t) = - / dM^'^{^j;(i,M^')Im/z'+(VF') + ^jKi,-W^')Im/(',+i)_(VF')} 

+ ^/d*' E '^^^y^ — i?5^,(t,0im/_^'(t') vj>i, 



(3.65) 



j' J'=J+2 



for the inhomogeneities A^(f), which besides the t-channel projections N^{t) of the nucleon pole terms 
contain the couphng to all s-channel partial waves as well as to the higher i-channel partial waves. 
Note that A;[(t) only contains the left-hand cut and therefore is real for all t > t.,^ . By virtue of (j3.57p 
and the analogous definition 



J -.J 



^rit) = m^j—jAi{t) - Aiit) , (3.66) 

the equations p.64p can be cast into the form of a MO problem for f^{t), f_{t), and the linear 
combinations r"'(t) 



00 

(+i\ 



;;,).<(,, i^/a.^^ 



{t'-tN){t'-t) ' 



00 



r^(,)=A^(t)+^^/dt' ,■■"";,";' , vj>i 



00 

fl{t) = Aiit) + ^Jdt' ^'^j2^P VJ>1, (3.67) 
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where for /+(i) and T (t) combining the integrals effectively yields one subtraction at the threshold 
tN and the additional roots at t' = tN in the denominators are canceled by the threshold behavior of 
the numerators. The solution for /+(i) can then easily be recovered via (j3.57p . 

How these equations (or their subtracted analogs derived in Sect. 14. 4p can be used to determine 
f+it), /±(i)) and f^{t) with the help of MO techniques will be described in the following sections. Note 
that such an easy rewriting scheme is not possible for the s-channel part ()3.4p of the RS PWHDRs, 
since in the corresponding s-channel integrals also < I' < I contribute. 

4 Subtracted Roy— Steiner system for pion— nucleon scattering 

The Proissart-Martin bound [57^158) limits the number of subtractions necessary for the convergence 
of the integrals in the high-energy regime to 2, since the total cross section does not increase faster 
than log s for s — )• ooj^ The influence of the high-energy contributions to dispersion integrals may be 
reduced by means of suitable subtractions for the trade-off of introducing corresponding subtraction 
polynomials with subtraction constants that are a priori unknown. For the MO integrals in (j3.67p 
subtracting in t at subtraction points below t-,^ with the additional constraint s = u in order to 
preserve crossing symmetry is favorable. A particularly useful choice is the subthreshold expansion, 
which amounts to subtracting in t at zero: first, it is very convenient for extrapolation to the Cheng- 
Dashen point in order to elaborate on the vrA^ a term (cf. Sect. 12. ip : second, subtracting at the 
subthreshold point facilitates matching to chiral perturbation theory, which is expected to work best 
in the subthreshold region^ To this end, we first briefly review the subthreshold expansion of the 
scattering amplitudes and then discuss its application in order to write down both the once- and 
twice-subtracted form of the HDRs ()2.44p . 

4.1 Subthreshold expansion 

The subthreshold expansion refers to the expansion of Born-subtracted amplitudes around the sub- 
threshold point {s = u = so,t = 0) = {v = 0,t = 0) (cf. Sect. 12. ip . where the nucleon pole terms 
are subtracted since they are rapidly varying in this kinematical region. Subtracting the pseudovector 
Born terms (indicated by bars) yields 



A+{s,t) = A+{s,t) - ^ , B+{s,t) = B+is,t)-g^ 
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1 1 



A- (s, t) = A- (s, t) , B- (s, t) = B~{s,t)- g 



'- - s 
1 



+ 



+ ^ > (4-1) 



while for the pseudoscalar Born-subtracted (indicated by tildes) amplitudes A^ and B^ the terms 
—g^/m and ^g"^ /Im? need to be dropped (cf. I2.48p . Due to the crossing symmetry of the ampli- 
tudes (|2.18p (similarly for HDRs (|2.44p ) one can write the subthreshold expansion generically for 



^^While the original Froissart bound assumes validity of the Mandelstam representation for the scattering amplitude, 
the result by Martin is based on somewhat less restrictive assumptions. 

'^^For the application of heavy-baryon ChPT to tvN scattering in the subthreshold region see [59]. Conversely, analyt- 
icity and unitarity are used in [BD] to stabilize the extrapolation of nN partial waves derived from ChPT amplitudes in 
the subthreshold region into the physical region, thus enabling the determination of the chiral parameters by matching 
to experimental information in terms of s-channel phase shifts. 
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crossing-even amplitudes as (cf. p9] ) 

^-^ ^ ' V V V V V V \ 

m,n K. J 

(4.2) 
and thus explicitly for the pseudovector Born-subtracted amplitudes as 

oo oo 

m,n=0 m,n=0 

oo oo 

A-iu,t)= Y ^^n^''^^'t\ B-{u,t)= Y b-r^n^^'^e . (4-3) 

m,n=0 m,,n=0 

where the corresponding subthreshold parameters of the amplitudes D = A + uB are related by 
Note that due to 6^^ „ = in particular 

< = 4n ■ (4-5) 

Prom the expansions 

1 1 Amv ,^ / q \ 1 1 2 t .^ / o 9 o\ 

1 — = r' — ^ + o{u\ut) , - — + - — = — ^ + o{,^\i.H,t') , 

s — s s — u [s — soj s — s s — u s — So [s — soj 

(4.6) 
one then can read off the subthreshold expansions of the Born-unsubtracted amplitudes up to and 
including first order 

(,2 

A+{u,t) = ^ + d+Q + d+^t + 0{iy^,iyh,t'^) , 

A-{u,t) = ua^Q + 0{v'^,vt) , B+{u,t) = g^-- -^ + vh+Q + O {v^^ , vt) , 

' [rri'^ — So) 



B-(v,t)=g^ 



m? — sq {m? — so)^ 



^ + b^, + h-,^t + 0{v\v\f) . (4.7) 



4.2 Sum rules for subthreshold parameters 

Subtracting simultaneously at sq = S/2 < s+ and to = < tjr corresponds to the subthreshold 
expansion around (i/ = 0, t = 0) and thus allows for the determination of sum rules for the subthreshold 
parameters. Matching the expansions (|4.7p to the corresponding expansions of the HDRs (|2.44|) by 
equating the coefficients (where it is crucial to keep track of all implicit dependencies in the expansions) 
together with introducing the abbreviation 

ho{s) = -^ ^ (4.8) 

S — Sq s' — a 

then yields the following sum rules for the lowest subthreshold parameters 

oo oo 

4o = -^ + 11 ds' hois') [lmA+{s',z',)]^^^^^ + i||: [im^+(t',4)] 



(0,0) ' 

■5+ t^r 
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-"00 



^01 



2m? 



m^ — a 



+ lJas'ho{s')[lmB-{s',zi)]^^^^^ + ^l^[lmB-{t',4)] 



(0,0) ' 



«+ 



[lmA+{s',z',)] 



is' - soY 



■^01 






(5' - so)^ 






"00 ^ ^ fdc' ^^'^^ (■'^''^^)](o,o) 1 /"df 
4m~ TT J {s' - soY vr y t' 

4r?T, vr 



Imyl"(t',zP 



■00 _ - / j^/ L v-:-.vj(0,0) _^ 1 f dt 



is' - sof 




(0,0) 



(0,0) 



(4.9) 



The subscript (0, 0) indicates that z'^ and z'^ in the s- and t-channel integrals, respectively, are to be 
evaluated at {v = 0, t = 0), which according to (|2.52p and (|2.43p amounts to using 



Ws. 



'2' 



1 



is' - so)' 



(0,0) - 2g'2(s'-a) 
t'it' - 4(so - a)) 



I (0,0) 



/2^/2 



W/q't 



1 + 



t'Aa — tNt-K 



[dtz't 



_ sq- a 
(0,0) - 2q''^{s' -a) 
m So -a 



(0,0) ^2^. 



.12 ' 



(4.10) 



where again we have used the fact that the t-channel integrands depend on the squared angle z'^' only. 
Note that these sum rules as such are valid independent of the choice of a, but in practice one will incur 
an a-dependence once approximations are made (such as truncation of the partial-wave expansion, 
approximation of the high-energy region by Regge theory, etc.). 

4.3 Subtracted hyperbolic dispersion relations 

A single subtraction at (z^ = 0, t = 0) only affects A^{v,t) and i?^(i/, t), since both A~{v,t) and 
B~^{v,t) are proportional to v. Based on the unsubtracted HDRs (|2.44p . the explicit subthreshold 
expansions (|4.7p . and the corresponding sum rules (|4.9p . we obtain the once-subtracted HDRs 



^+(s,t;a) = — + 4o + - dt'l 
m IT J 



1 r .flmA+(t',z^) [ImA+(t',zO] 



(0,0) 






+ 



t' -t 
1 



t' 



(4.11) 



lva.A^{s,Zg) — ho{s) [lmyl^(s', z^^ 



I (0,0) 
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B {s,t;a) = g^ 



1 1 

+ 



m? — s m? — u 



s' +.-„+! /dJ'-^~(''-» [i-B-(*'.--;)i,„,„, 



2m^ 



^00 



vr 



t'-t 



t' 



+ - I ds'< 

IT 



1 1 

+ 



luiB {s , z'g) — ho{s') [luiB (s',z^)] 



(0,0) 



together with the unaltered equations (12.44p for A~ and B~^. Note that the dependence on a of the 
Born-term contribution N~ is canceled by the sum rule (j4.9p for 6qq, which is why the subtraction 
constants are formally included in the subtracted nucleon pole terms in the following for convenience 
(i.e. preserving the generic form of the HDRs (|2.44p ). 

Similarly, a second subtraction at {u = 0,t = 0) yields the twice-subtracted HDRs 



A+{s,t;a) = — + d++d^,t 
m 






+ - I ds'< 

IT 



1 



Im A+{s',z',)- hois') t[dtImA+ is', z',)]^^^^^ 



hois') 



is' - so)' 



[lmA+is',z',)] 



(0,0) 



B is,t;a) = g' 



+ 



m? — s m? — u 



^ + &00 + bo,t 



+ - fdt'l '^B:('':<) -(Ij + ^) [Imi?-(t',z0],o,o) - ^ [5JmB-(t',zO](o,o) 



IT 



t'-t 



t' i'2 



I (0,0) ^z 



+ - / ds'< 

vr 



+ 



ImB is' , z'g) — hois') t[dtlTaB is' , z'g)~\ 



(0,0) 



A is, t; a) = ttQQU -\ — dt 



hois')- (y_^^)2 ) P"^^ (^''^^)](o,o) I ' 
.j Im A- it', z't) [lmA-it',z[)/i^'] 



IT 



(0,0) 



u'it' - t) 



t' 



+ - / ds'< 

vr 



lmA-is',z',) ^ ^ '^^^°'°^ 



(s' - so)' 



B+is,t;a)=g' 



w? — s m? — u 



+ .V + ^/d<'^-^"<^-^i' [i-B-(*'.4)M 



^00 *■ 



TT 



/' 



(0,0) 



v'it> - t) 



t' 






'imv[lm.B+is',z'y\(n.n\ 

imi?+(.',4) — ^. ^ :,^'^'^ 

(s - •so)^ 



(4.12) 
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where A and B^ can also be written as 



A-is,t;a) ^ Ooo^l /"^^/J lraA-{t',z[) 



Amv 



4m vr 



Ap[q[z[{t'-t) t' 



I J ~ 



Im^ i^s' ,z\ 







(0,0) 



is' - sor 



B+{s,t;a) 
Amv 



{m? — s){m? — u) 4m vr 



+ :^ + i /"dt'i i"i^^(*''^*) 1 



4pWt4(i'-*) *' 



Img+(t^z;) 



(0,0) 



1 7 ,[ Imij+(s^4) [Imi?+(g^4)](oo) 



(4.13) 



These subtractions require a modification of the nucleon-pole-term projections and the kernel func- 
tions for both the s- and t-channel contributions calculated in Appendices [X] and [B] as well as the 
asymptotic contributions given in Appendix [Dl The differences on the right-hand side of the once- 
/twice-subtracted HDRs (|4.1ip / (j4.12p compared to the unsubtr acted HDRs (|2.44p are the sources 
for the necessary modifications which are derived in Appendices IA.51 IB.51 and ID. 31 for the s-channel 
kernels and pole terms, their t-channel analogs, and the asymptotic contributions, respectively. 



4.4 Subtracted t-channel Muskhelishvili— Omnes problem 



Using the subtracted kernels and pole terms as derived in Appendix lB . Sl leads to the subtracted analogs 
of the unsubtracted t-channel MO problem (|3.64p . which we will state explicitly in the following for 
J < 2 (the equations for J > 3 are unaltered for up to two subtractions). For one subtraction we may 
write 



^+^*^" +' ^^ vri ^S'(t'-t^)+vry ^S'(t'-t) 



oo 



1,1-sub,., , t^ 7^^,^Im/i(t')-Im/|(tO , fAmfiit' 



flit) = A\r-\t) + 



vr 



dt' 



OO 



t'{t' - tN) 



t-jr 



1 |i-sub^^^ , t f ^^ilm.ft{t!)_ 

t) 






' vr 7 t' - tAf vr y t' -t 



f_{t) = Al\'-\t) + IJ dt'^-^ 



2/+/^ 



(4.14) 
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while two subtractions yield 



oo oo 

00 
2 r T^ ^1 /'+M 



/i(t)^Air--(.)+^/dt'^ 



t) 



flit) = Air^it) + ^ / dt- V- ^^^^^- -^ ^ ■ ■ + ^ / df 



vr J t'{t' -In) tt J t'{t' -t) ' 



2-sub,,, , % 7_^^-Vl!:-^ (*')-I-^+(*') , t f^.lrnflit 

m = Ai\'--\t) + l]di''^r^'' 



t) 



(4.15) 



It is important to note that S- and D-waves are coupled, as A^ contains contributions from J = 2 
according to (|3.65p . While the integrands containing the Cauchy kernel in (j4.14p and (j4.15p for J = 
and J = 1 clearly show the corresponding number of subtractions at to = 0, for J = 2 there is always 
one subtraction less or no subtraction at all. Note that the integrands containing linear combinations 
of the partial waves are proportional to t^ /t' (if affected by the subtractions at all), which results 
in a suppressed internal high-energy dependence inside the integral due to division by t' without an 
increased external high-energy dependence due to multiplication with t^ rather than t as for a usual 
subtraction at zero. 

The un- p.64p , once- ()4.14p , and twice-subtracted ()4.15p equations are of the original form of the 
(subtracted) MO problem with integrals of the absorptive parts times the Cauchy kernel, if the re- 
maining t-independent integrals (which may, however, come with t-dependent prefactors) are absorbed 
into a redefinition of the inhomogeneities A;{.(t). This problem is well defined due to the threshold 
behavior of the partial waves at t = tisi- However, the price for taking advantage of the convergence 
properties of the integrals this way is that reasonable approximations for the starting values for the 
partial waves are needed as input, since the solutions can only be found iteratively. 

Therefore, we prefer to utilize the threshold behavior of the partial waves and use the linear 
combinations T^ {t) in order to rewrite the equations in analogy to (j3.67p . i.e. to modify the original 
form of the (subtracted) MO problem in a well-defined manner. The general n-times subtracted (with 
n G {0, 1, 2}) versions of the MO equations (j3.67p for all J then read 

00 
■'^^ ' +' ^ ' vr J f^it' -tN){t' -t) 

T^(t) - A^r-^tl + t("-^)^(-^)(t-%) 7 ImT^jt') 

r(t)-Ar| {t)+ - J ^\,in-J)ein-j)^t,_t^^^t,_t^ ^-^^1' 

00 

fit) = A^ r-^""(f) + / dt' , ,!_Z \ ' T V J > 1 , (4.16) 
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where the Heaviside step function is to be understood in its right-continuous form, i.e. 9{0) = 1. Again, 
the equations for /+(t) and F (t) exhibit one additional subtraction at tjy, such that the combined 
number of subtractions for all J > can be given as (n — J + l)8{n — J). For convenience, we explicitly 
show the terms in A5|_(t) that couple the D- to the S- waves 



oo 



A^rV) = -Ai_!^|dt'{[t' + t-(t^ + t.) + 6a]^^^ + ^Im/!(0} + ..., (4.17) 

oo 

in 

oo 

» n i2-sub/ . 5 t — tN f dt' ( 1 r, , . , . ,-, ImF^(t') m ^ ^o , ,^^ 

All (*^ = -I^^r^y ■j^^-[it' + t)tNU-t't{tN+U)]-^r^ + tu^lmfl{t')j+... , 

which converge for t' — )• tjy due to the threshold behavior of F^ and vanish for i — t- tjv due to the 
exceptional behavior of the (n-times subtracted) kernel Kq2 i^^- (|B.52p and ()B.68P ): the respective 
remainder denoted by dots above is then given by 



1 r 
ir''\t) + l / d^'X^{Goz|"-^^^i,H-Olm/i(H^') + Go.r"^'^'(t,-H-')Im/(|,i)Ji^')} 

oo ^ 

+ ^ / di' E '-^¥^ [KL'^'it^ Irr^fi'it') + K^' ^^t, Im fl\t')} . (4.18) 



Note that this D- to 5-wave coupling becomes independent of a by subtracting once or twice, while 
the corresponding F- to P-wave coupling also depends on a in the once-subtracted case (cf. (|B.67P ). 

5 Solving the t-channel Muskhelishvili— Omnes problem 

In this section the solution of the MO problem for the lowest i-channel partial waves f±{t) with 
J G {0,1,2} will be discussed. First, the explicit analytical solutions will be stated. Then, the 
numerical input needed will be collected. Finally, the numerical results will be discussed. 

5.1 Muskhelishvili-Omnes problem with finite matching point 

We assume to know the imaginary part of the t-channel partial waves f±{t) above the finite matching 
point tm as well as the scattering phases Sj{t) of the vrvr partial waves tj{t) for 4M^ = t^^ < t < t^, 
which in the elastic region are also the phases of the /^(t) due to Watson's final state theorem, 
cf. (j3.37p . All inelastic contributions will be neglected. Under these assumptions, we have to solve 
equations of the MO type 



tin OO 

,, ^ A/N 1 f . ,T(t')*f{t') 1 f^,lmf(t') ,, ,, 

for /(t) in the range t.,^ < t < t^ with finite im [E], where the physical values of the integrals are 
obtained in the limit i — )■ i + ie and the discontinuity of f{t) across the right-hand cut is given by 
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unitarity (cf. hermitian analytic! ty ()3.17p and the elastic t-channel unitarity relation (j3.36p ) 

^^^[^ = Im fit) = ntyfit) e{t - t^) , (5.2) 

where the inhomogeneity A(t) contains potential left-hand cut contributions to f{t) (i.e. it is real for 
tn < t) and the elastic amplitude T(t) is given by 

r(t) = sin5(t)e*'^W . (5.3) 

We briefly review the result of |19y27] in the following sections. 

5.1.1 General solution 

To begin with, we consider the homogeneous problem for a function fo{t) with non- vanishing imaginary 
part only for i^ < i < tm- The solution can then be written as 

fo{t) = n{t)^o{t) , (5.4) 

with the Omnes function"^ ^1 

1 1 YT^t \ = '"^*^' """"P {i6{t)e{t - u)e{t^ - 1)} , n{o) = i , 

^/f^| = l^(t)||t.-rW, x(t) = ^, 



\m\ 



'fit-u) 

i'jr 



tjT 



where we have analytically separated the endpoint singularities of the principal value integral. By 
assuming the reasonable asymptotic behavior /o(i) — )• for t — )• oo (cf. (]3.47p ). the only analytic 
structures of So(i) allowed by fo{t) and i^{t) are poles at the endpoints t = t^^ and t = tm- Since Cl{t) 
is regular at tyr due to 5(1,^) = 0, the regularity of fo{t) excludes poles at t = t^r and restricts the order 
of the poles at t = im to 

n = [x\ , X = — ^ (5-6) 

TT 

{[x\ denotes the largest integer < x). In this way, we find 

So(t) = pzl^ , (5.7) 

where Vn-iit) is an arbitrary real polynomial of degree n — 1 that introduces n free parameters to the 
Omnes problem. For n = the homogeneous solution vanishes according to V^iit) = and no free 
parameter enters the problem. 
The general solution reads 

-Pn-ijt) 1 f A{t')sm6{t') 1 / . Im/(f) ] 

t-TT try, 



^^Note that for a finite matching point it is not mandatory to work with a subtracted Omnes function. However, 
subtracting once at i = ensures the usual normalization Q.{Q) = 1. 
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which may also be written in terms of a principal value integral as 



fit) 






A(t)cos5(t) + |^(t)||^^--^ + -^dt|^^^,^l^^,_^^ + -ydt^^^^^ 



e*'W ^ (5.9) 



in accordance with [27J for tm — t- oo. Note that due to Watson's theorem (13.37P the prefactor in square 
brackets can be identified with the modulus |/(t)|, and since the phase 5{t) is known, we only need to 
solve the MO problem for this modulus for tj^ < t < t^. 

5.1.2 Subtractions 

If X > 1, suitable subtractions need to be performed in ()5.9p to ensure integr ability for t' — > tm. Let 
us begin with the case 1 < x < 2, i.e. n = 1. We may write 



oo 



^^, Imfjt') _ 1 j, fdt' lmfjt') , ^ fdt' t^-t'lmf{t') \ 

""^ \n{t')\{t' - 1) - t^ - tV^ J t' Mt>)\ +V t' Mt')\ t'-t f ' ^^-'"^ 

••m ''in I'm 

where the second integral is now convergent. The first integral is still divergent, of course, but it does 
not depend on t any more and can thus be absorbed into a redefinition of the (constant) polynomial 
Vo in (j5.9p due to the common prefactor (t^ — i)~^- For higher values of x this subtraction and 
redefinition prescription needs to be iterated, whereby all n parameters contained in the polynomial 
receive corresponding contributions. Applying this reasoning to both integrals of (|5.9p for general x 
and using the highest number of subtractions allowed by the degree of the polynomial, the result is 
given by 

/ (-in 



|/(t)| = A(t)cos5(t) + J^M_ j,„_,(i) + ^^ 



ttJ f" \n{t')\ t'-t 
r fdt' {t^-tri^f{t') \ ..... 

In order to reduce the influence of the high-energy contributions on the Omnes integrals, subtrac- 
tions may also be introduced already right from the beginning (15. ip . With / such subtractions, the 
analog of ()5.1ip becomes 

ifmi Am .m^J^MLln^ ,^, , t-y dt' it^-tr^{t')sm5it') 

oo -v 

, t" f dt' (V-tO"Im/(O l 
^ n J t'^+i \Q{t')\ t'-t j • ^^ > 

This constitutes the final general result, which for / = 1 and n G {0, 1} reduces to the results quoted 
in pjj. 
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5.1.3 Numerical treatment 

The asymptotic beiiavior of tiie Omnes function \^{t)\ for t ^ t^ requires some care in the numerical 
evaluation of the integrals in (j5.12p . Although by construction the singularities for t' — )• tm are 
integrable, the corresponding cusps generate large contributions to the integral and a fully numerical 
treatment would require a very careful distribution of mesh points in order to catch the effect. In 
the following, we will demonstrate how these endpoint singularities can be separated analytically (cf. 
the appendix of |19j). For the sake of simplicity, we discuss here the case of n = / = 0, which 
already displays all relevant features; the generalization is then straightforward. To this end, we split 
the integrals close to the matching point im and approximate |il(t)| by its asymptotic form in the 
proximity of tm 

io(t«t^)i«ifi(t^)iit^-tr . (5.13) 

For r — >• O"*", we may thus rewrite the integrals above the matching point as 

OO OO im+T 

^^, Imfjt') r ^^, Im/(f) ^ Im/(t^) f dt' 



\n{t')\{t' - 1) J \nit')\it'-t) \n{t^)\ J \t^-tr{t'-t) 



OO 



^^, Im/(f) , Im/(t^) ^ ^^^ 

tm + T 

and similarly below the matching point either for tj^ < t < t^ — t 

2 , Ait') sin 6{t') _ r dt' /A(f) sin (5(f) A{t) sin 6{t)\ A{t)sm6{t) tm - t - t 

J * \n{t')\{t'-t) ~ J ¥^t[ \niF)\ Mt)\ J+ \n{t)\ °^ t-u 

A{trn)sm5{t^) 

+ i«(uifa -» '-'"• '"" 

or for tm — r < t < tm 

7 A(tOsin,5(tO _ f dt' A{t')sm6it') Ajt^) sin 6{t^,) ~ 
J |0(t')|(t'-t) J t'-t |0(t')| ^|f^(t„,)|(tm-t)^ '^^' ^" ^ 

The substitution v{t') = {t' — tni)/(im — t) leads to the integrals (with x G (0, 1)) 



V''{1±V) 1-X J l±v' ^' tra-t' 



r{t) m 

J v^[l — v) 1 — x J 1 — v 



Separating the singularities as shown above, the remaining integrals can be solved by using standard 
integration routines. For sufficiently small r (i.e. if r is of the same order of magnitude as the 
discretization error of the integration routine), the above approximations are well justified and this 
procedure allows for a stable numerical evaluation of the Omnes integrals. 
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5.1.4 Continuity at the matching point 

The continuity of the Omnes solution /(t) at the matching point tm is analyticaUy ensured by the 
asymptotic form of the corresponding integrals of (j5.17p (with < x < 1, cf. the appendix of |19] ) 

oo oo 

/I PA 
; 7 = VrCOSeCTTX , /_(tm) = T 7 T = — TT COt TTX . (5.18) 
u^(l + z;) ' J V'il-v) ^ ' 



Taking equation (15. 9p in the limit t — )• tm from below, plugging in these asymptotic expressions for 
the integrals, and using (j5.6p indeed reduces the square bracket to |/(tm)|- This analytical equality 
may also be used as a check of the numerical evaluation. 

However, the continuity of the first and higher derivatives is not ensured in a similar manner. Since 
the solution must not depend on the value of the matching point, an unphysical cusp or non-smooth 
behavior of the modulus of the solution at the matching point only indicates that the input in terms 
of the absorptive part is not precise enough; moreover, the physical condition of a smooth behavior at 
the matching point ensures the uniqueness of the solution [32j. Physically consistent input given, this 
smoothness constraint may actually be used in order to tune/estimate/fit the subtraction constants 

(cf. m)- 

5.2 Explicit solution of the t-channel Muskhelishvili-Omnes problem 

Here, we will give the explicit solutions for the n-times subtracted t-channel MO problem ()4.16p 
using the general results of Sect. 15.11 The crucial ingredient for the following discussion is Watson's 
theorem (j3.37p . which states that below the onset of inelasticities the phases '-p'^ii) of the i-channel 
partial waves /±(t) are given by the corresponding vrvr scattering phases 5j(t) with It G {0,1}, i.e. 
exphcitly for J E {0, 1, 2} 

^\(t) = S^oit) = Soit) , viit) = 6l{t) = 6i{t) , (^l(t) = 50(t) = 62{t) . (5.19) 

These identities enter the solutions at two places: first, in this kinematical region we can use the same 
Omnes function Qj for both /^ and thus also for the linear combination T . Second, in this range of 
t the linear relation (j3.57p is also valid for the moduli such that after solving for IF'^I we can recover 



|/+(*)| = "^\/tXtI/-(*)I - |r^(*)l • (5-20) 



J + 1 
Using once-subtracted Omnes functions in the convention (cf. ()5.5 

dt' 6j{t' 




nj{t) = exp I - / —Yzi \ = l^-^(*)l ^^P {'^j(.t)Kt - *-)^(*m -t)}, f^j(O) 



1- 


t 


=^j(t) 


t 




tin 




U~ 


- 1 




lo ml 1 ^ ^ 1^^'^*^ J t rdt' 5jit')-5j{t) \ 6j{t) 

|^jW|=1-— r-1 exp<;-/— j^—^ \, xj{t) = —^, (5.21) 
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the general n-times subtracted (n G {0,1,2}) solutions of (|4.16p for t £ [tn,tm\ and [6j{t^)/Tr\ = 
read 



.oin-sub,,, , ^ ,^,t^{t-tN)\ y^. AO+f'-^"^t')sm5o(t') 
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Now, we can use the spectral representations of the inverse of the Omnes functions in the un-, once-, 
and twice-subtracted form 



TT J t'-t ttJ \nj{t')\{t' -t) 

t f ^^, sin 5 j{t') ■ e f^^, smdjjt') 

with the derivative of the Omnes function^^l 



Oj(o) = ^nj{t) 

dt 



i/d.Mp, (,24) 



*=° i 



vr ./ t'2 ' 



in order to explicitly perform the integrals over terms that are either constant or come with appropriate 
factors of t' or p'^, i.e. all terms involving the subthreshold parameters as well as the term proportional 
to 5ji/{m? — a) for the unsubtracted case. For this purpose we define A|[(t) via removing all constant 
or subthreshold-parameter contributions from the inhomogeneities A^(t) (cf. (|3.65p and (]B.6ip ) 

A^f ^t) = AiP'^^t) - AiVif ^t) = Ni{t) + Air-^"'^(t) , (5.25) 



^^Note that for fm — >■ oo (and neglecting inelasticities in the single-channel approximation) this quantity is closely 
related to the pion vector radius for J = 1 : limt„_>oo fii ~ \{'r'^)X ■ 
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and thereby we obtain 
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Note that also in the unsubtracted case the expHcit dependence on a cancelso 



^^Actually, this has to be the case: e.g. the constant term proportional to {m? — a)~^ in the nucleon pole terms (|2.48p . 
which was introduced to the dispersion relations via the hyperbolic kinematical relations and which can be thought of 
as a contribution of the contour integral from the circle with infinite radius, leads to constant pole-term contributions 
to the partial waves (cf. (|B.6|) and HB.62|| '). These (unphysical) contributions do not vanish asymptotically, generate 
an unphysical behavior on a, and thus they must cancel in any (physical) solution. Hence, the dispersion integrals for 
the unsubtracted case both for the Omnes solution and the spectral representation of the Omnes function are strictly 
speaking not correct: there should be contributions from the contour at infinity. However, this problem can be solved 
most easily by removing all "dangerous" parts of the inhomogeneities via (|5.23p . which ensures that all these potential 
contributions from the contour at infinity cancel. 
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Finally, due to Watson's theorem (|3.37p we can separate the unknown moduli from the known vrvr 
phases and solve the MO problem for the moduli directly 
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(5.28) 



On the one hand, the subtraction-independent pole terms iV^ are real for t > t.^- — (M^/m)'^ and 
grow rapidly with J for t in the vicinity of tjr, as discussed in Appendix lB.il On the other hand, 
in the elastic region the phases 6j are given by the corresponding vrvr scattering phases such that 
Sjit-ir) = and thus lTaf^{tn) = 0. Since furthermore phenomenologically the vrvr phases grow slower 
for higher J, we thus expect the partial waves f^ (and thereby also their moduli \f±\) to be increasingly 
dominated by the pole terms for increasing J and t ^^ t-n-. However, we do not solve for ff directly 
but for the linear combination F , for which, in turn, the pole-term contributions N^ cancel at tn, 
cf. Appendix IB. 51 The pole-term domination of |/:|^| enters when calculating these parallel helicity 
moduli from the solutions for the |F'^| and the (pole-term dominated) antiparallel helicity moduli \f£\ 
via (|5.20p . where in addition the relative importance of the latter increases with J due to the factor 
^ J ji^J + 1). Both the pole-term domination and the dependence of |/j|l| on |/_| will be explicitly 
demonstrated in Sect. I 



5.3 Input 



In this section we will discuss all input that is needed to solve the (elastic) t-channel MO problem (j5.28p 
as given in Sect. 15.21 
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Figure 3: Moduli \^j\ of the lowest once-subtracted finite-matching-point Omnes functions for ^/t 
0.98 GeV. 



5.3.1 Pion pion phases and Omnes functions 

We use TTTT scattering phase shifts Sj{t) of [MlllB] for J G (0,1, 2} with It G {0, 1} (It even/odd for J 
even/odd) which are constructed for y/i G [2M7r, 1.15 GeV] cj Schenk-like parameterizations p^[62] 



ts.n6''{t) = aUt 



2J 



{A'j+B!;q^ + c!;qt + ...} 



671- 



Jt 



Jt 



(5.29) 



where the parameter rj denotes the point where the corresponding phase shift passes through 7r/2 
and the Schenk parameters Aj etc. may be related to the coefficients of the threshold expansion, 
ensure both the vanishing at threshold 5J{tT^) = as well as the correct square-root-power behavior 
above threshold. Thus, 6q is linear, 5\ cubic, and 82 quintic in aj. 

In Fig. [3] we show the moduli \^j\ of the resulting once-subtracted finite-matching-point Omnes 
functions according to (j5.2ip for J G {0, 1, 2}, where the choice y^ = 0.98 GeV ensures that xj{t) G 
(0,1) and hence nj = [xj(tm)J = for t G [i7r,im]- Therefore, all functions are normalized to unity 
at t = 0, finite for all i, and vanish at t = im due to the finite-matching-point prefactor \t^ — t\^-'^'^' . 
Furthermore, for J = the Omnes function exhibits a cusp (i.e. a discontinuity of the derivative) at 
the physical vrvr threshold t-,^ and decreases approximately linearly over a wide range in t, for J = 1 it 
is fully dominated by the /3(770) peak, and for J = 2 it is almost flat (equaling one again roughly at 
the end of the KH80 energy range at 0.88 GeV and dropping rapidly above). 

Using instead the parameterization of the vrvr phases as given in [14] for the numerical evaluations 
leads to deviations in these Omnes functions, and thereby the solutions of the MO problem (I5.28p . 



For TTTT scattering the validity of the Roy equations can be shown rigorously for i^ < f < 60 A/^ based on axiomatic 
field theory [11]. Assuming Mandelstam analyticity, this range can be extended to i^ < i < 68M^ [61], which corresponds 
to 2M,r < \/i < 1.15 GeV, by reasoning along the lines of Appendix [Cl for tttt scattering. 
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that, however, are much smaller than the effects of the alterations described in Sect. [5? 



5.3.2 General remarks on existing pion— nucleon partial-wave analyses 

Before summarizing the input from vrA^ partial-wave analyses that will be used in the following, some 
general remarks are in order: first of all, we will use the Karlsruhe-Helsinki dispersive partial-wave 
analysis KH80 |28II29) both as input for s-channel partial waves as well as subthreshold parameters and 
as reference for our MO i-channel partial- wave solutions, since KH80 is still the only consistent analysis 
for all the partial waves and parameters entailed in our RS framework. KH80 is based on ttN — )• 
i:N data only (and isospin invariance) and uses Pietarinen's expansion method [64j in combination 
with conformal mapping techniques, aided in particular by fixed-t analyticity Its solutions for both 
channels are given as tables in [29] o Moreover, an iteration uncertainty of about 3% for the iterative 
KH80 procedure is stated, which, however, cannot replace a thorough analysis of the systematic 
uncertainties. The subsequent Karlsruhe analysis KA84 [65] improves on KH80 especially for higher 
partial waves by using a modified PWDR framework and thereby smoothing KH80, but no consistent 
subthreshold parameters are derived in this frameworko The same holds true for the continuously 
updated VPI/GWU(SAID) s-channel analyses, see e.g. |43 U67f[69] . for which at most the vrA^ coupling 
constant and some of the necessary subthreshold parameters are determined. For the i-channel partial 
waves in the unphysical region t G [t7r,tAr], there also exists an unpublished solution [56] extending 
the KH80 energy range ^/l £ [2M^,^/40M^ = 0.88 GeV] to roughly IGeV. While this solution is 
compatible with KH80 within the aforementioned range, it seems to suffer from internal inconsistencies 
for higher energieso For the f-channel partial waves in the physical region t > t]\f, however, there 
exists a partial- wave analysis [TOj, which at least in principle could be used as input. Finally, a partial 
update of the KH80 analysis including new data and using more computational power was reported 
in [7l], but so far only results for forward vrA^ scattering have been published |72j . 

5.3.3 s-channel partial v^aves 

We use the KH80 solution for the s-channel partial waves from [29J for W^ <W< Wg, = 2.5 GeV. On 
the one hand, this is roughly the same energy range as for the continuously updated GWU "current 
solution" [69] such that we are able to compare between KH80 and GWU solutions as input. However, 
the effect of taking the GWU solution (or the "smoothed" KH80 solution [69]) instead as input for 
the i-channel Omnes problem (i.e. the corresponding inhomogeneities A^) turns out to be much 
smaller than the effects discussed in Sect. 15.41 On the other hand, at W^ = 2.5 GeV a reasonable 
transition from the truncated sum of partial waves below Wg, to the Regge model for the full invariant 
amplitudes above Wg can be achieved as we will demonstrate now. Summing up all partial waves with 
I < 5 would encompass all 4-star resonances of [34], but of both I = 5 4-star resonances, A^(2220) as 
HiQ and A(2420) as H^^n, especially the latter is mostly out of this energy range due to its broad 
width of roughly 700 MeV. Hence we expect the best agreement with the Regge model |73], which is 
based on differential cross section and polarization data for vrA^ backward scattering with W > 3 GeV 



^''As stated in [29], the Karlsruhe-Helsinki dispersive partial- wave analyses KH78 and KH80 (see Sect. 15. 331 for more 
details) use as input the mr phase shifts of [63j, which are based on Roy-equation fits. In principle, the differences 



between these phase shifts and the recent results I14H16] are sources of discrepancies between the KH80 results and the 
solutions of the MO problem. However, this point is of minor importance for the results discussed in Sect. 15.41 

^"^In ,2S] the results for the f-channel partial waves are quoted as KH78 solution, but according to [25] these tables are 
actually calculated from the KH80 s-channel solutions. Thus we will speak of the f-channel partial waves in [S^ as KH80 
solution as well. In general, KH80 is an update of KH78 including more recent data and improved fixed-i analyticity 
constraints. 

^^For a comparison of KH80 with KA84 and also an improvement of the formalism outlined in [46], see [66] . 

■^^There are e.g. rather obvious outliers (corresponding to unphysical jumps) in the phases. 
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Figure 4: Matching of the s-channel absorptive parts between KH80 partial-wave contributions for 
/ < 3, / < 4, and / < 5 and the vrA^ backward scattering Regge model [73] . 



as discussed in Appendix ID .H for / < 4 and a scattering angle of z^ = —1 corresponding to backward 
scattering. Since deviations between summing up contributions for / < 3, / < 4, and I < 5 start to 
show up around 1.5 GeV and we are interested in the matching to the Regge model at the end of the 
GWU range of validity around 2.5 GeV, only this region is shown in Fig. H] (in the spirit of [H]). Note 
that only / < 4 yields the correct sign compared to the Regge contribution in all four cases. Moreover, 
it turns out that for I < 5 the agreement is even worse than for I < 3. Hence, in the following all 
higher partial waves with I > 5 will be neglected below W^.- 



5.3.4 t-channel partial vi^aves 

The assumption of elastic unitarity breaks down in the 5-wave as soon as the KK channel opens, 
which manifests itself in the appearance of the /o(980) resonance. In principle, there are several ways 
how this phenomenon may be accommodated in a single-channel description. 

First, inelastic contributions could be included directly in the solution of the MO equations along 
the lines discussed in |74l[75j provided that the inelasticities are well known. However, in the case of 
/^ this would in particular require knowledge of the KK — )• NN S'-wave, but it is unclear how reliable 
input for this partial wave could be obtained independently from the present approach. 

Second, one could retain a rather low matching point t^, but try to model the energy region above 
im by means of a resonance description in order to establish a more meaningful matching condition. 
This strategy proved quite successful in 77 — )• tttt [31J, where the input above the matching point is 
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Table 1: Subthreshold parameter values as given by KH80/Hohler [2U], Stahov [M|, Oades 
and Fettes (heavy-baryon ChPT) [85]. See main text for details. 

dominated by the /2(1270). However, in the case of the /o(980) this strategy is subject to several diffi- 
culties: its pole position is very close to the two-kaon threshold, such that the subtle interplay between 
the TTvr and KK channels can certainly not be approximated by a simple Breit-Wigner description. To 
circumvent this problem, one would be compelled to further decrease the matching point and include 
the /o(980) dynamics by hand using a Flatte-like parameterization [76], which is a modified relativistic 
version of the Breit-Wigner differential mass distribution. However, while the fomr coupling constant 
has been thoroughly investigated ^77j based on the recent dispersive analysis [Tl] (which yields phases 
that are basically consistent with the phases of |15tll6j). the foNN coupling constant is only very 
poorly known, with different meson-exchange models disagreeing significantly on the strength of the 
coupling and the continuation to the physical pole [78l - [8T] . We conclude that including the /o(980) in 
our approach reliably as well as extending the energy range of our representation for /^ beyond the 
two-kaon threshold will require a full solution of the underlying two-channel Omnes problem |82j . In 
this work we will content ourselves with the single-channel approximation. 

Since therefore we can solve the single-channel MO problem in the elastic region only and further- 
more iteration with the s-channel RS solutions (for which, in turn, accurate MO solutions are needed 
as input) as well as a consistent determination of the vrA^ coupling and the subthreshold parameters 
is necessary to finally arrive at precise quantitative results for the partial waves of both channels, here 
we will only give qualitative results for the t-channel partial waves by comparing with KH80. 

Hence, in the following all t-channel absorptive parts above t^ are set to zero. Consequently, all t- 
channel Regge contributions are omitted (since tm < ^a! ci. the discussion of the i-channel asymptotics 
in Appendix ID. 2p . Note that otherwise one would have to avoid double counting of the asymptotic 
regions of the t-channel partial waves in the MO problem. Finally, also all higher partial waves with 
J > 3 are neglected. 

5.3.5 Subthreshold parameters 

To precisely determine the subthreshold parameters is not an easy task, since there simply is no ex- 
perimental data available to analyze the t-dependence of the amplitudes close to t = and thus means 
of analytic continuation or extrapolation are needed. Accordingly, in the literature there are only few 
determinations of all parameters that enter the subtracted RS system. The KH80 results (cf. [29], 
wherein the error estimates are quoted to be "based on deviations from the internal consistency" and 
the total uncertainty to be "somewhat larger" ) and all more recent dispersion theoretical analyses that 
we are aware of are collected in Table [Tj (cf . [2QJ ) . Note that there are several determinations of only 
some of these parameters, which are therefore not listed in Table [TJ In [23] the subthreshold parame- 
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ters are determined by means of interior dispersion relations together with fixed-t dispersion relations, 
and by using as input the s-channel partial waves of both KA84 and VPI/SP98 [67", 86] as well as 
the i-channel partial waves of KH80 (and those of [56] in the consistent energy range) o In contrast, 
finite-contour dispersion relations are used in [83] to derive subthreshold parameter values — again for 
both KH80 and VPI/SP98 input (amongst others) IB 

The subthreshold parameters are the standard expansion parameters for the Lorentz-invariant 
amplitudes, but neither these amplitudes nor the kinematical variables v and t are natural for heavy- 
baryon ChPT, and hence values obtained in analyses using this framework are not very satisfactory, 
cf. [84]. However, for comparison we also state the corresponding values for a third-order calcula- 
tion [85l[88] as given in [85] (Fit 1 therein corresponding to KA84); note that according to [85] [89] 
some of the parameter values even deteriorate when calculated up to fourth order. 

As can be seen already from the deviations between the different determinations of subthreshold 
parameters in Table [U the errors on the central values are in general unrealistically small (i.e. only 
statistical fit errors for specific input in a given framework, thus neglecting systematic errors). Hence 
we can conclude that there is no precise and consistent determination of the subthreshold parameters 
including realistic errors. Since we want to compare our MO results with the KH80 solutions, for 
consistency we use the KH80 subthreshold parameters as given in Table [1] as well as the outdated 
KH80 vrA^ pseudoscalar coupling value of 14.28 instead of the new value of 13.7 as given in (I2.49p ]^i 

5.4 Results 

The numerical results that will be presented in this section are to be understood as a qualitative "KH80 
consistency check" in order to show that the i-channel RS-MO machinery works, and as a first step 
towards a numerical analysis of the full RS system. In particular, by variation of either the coupling 
or the subthreshold parameters we can alter the results significantly, since these variations produce 
the most sizable effects on the MO solutions compared to the other variation that will be discussed 
in the following. However, it is by no means clear a priori what the parameter values or their errors 
are, and only a self-consistent determination of all parameters and partial waves in a second step will 
allow for reliable quantitative results. Therefore, the necessary first task in this program is to check 
our method and the internal consistency of the KH80 results by using KH80 input as described in 
Sect. 15.31 and comparing our i-channel MO results with those of KH80. Moreover, we will investigate 
different systematic effects on the (subtracted) MO solutions | /^^ | , which should prove valuable for the 
solution of the full system: after discussing exemplarily the importance of the different contributions 
to the MO inhomogeneities A^(t), we will also discuss both the connection to the "fixed-i limit'o 
a — >• — oo and the effect of changing the matching point tm- Except for the a — ?■ — oo results, we will 
always use the optimal hyperbola parameter value of a = —2.71 M^ as obtained in Appendix IC.4[ 



^^These are most probably the "new subthreshold parameters" mentioned in [S7|, where no explicit reference is given. 

^''Note that some of the results of 'WS are corrected in ^HT, where also a modified version of the finite-contour dispersion 
relations together with conformal mapping techniques is applied (it is mentioned therein that the subthreshold parameters 
do not change substantially). Since the applied fitting procedure does not respect the exact analytic equality of the 
parameters dj„ and aj„, however, the (corrected) values agree only within the given errors, but not exactly. 

^^Note that the ttTV coupling and the subthreshold parameters are related, as the difference d^Q — g^ /{2m) is given by 
an integral over a total cross section, cf. |20| . 

^® Accordingly, the s-channel integral of the HDRs reduces to the fixed-i result, cf. Sect. 12.31 However, even in this 
limit the HDRs contain additional information as compared to fixed-f dispersion relations, since those do not provide 
equations for the t-channel partial waves in the first place. 
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Figure 5: Contributions to the unsubtracted S-wave MO inhomogeneity A5_| . See text for details 



5.4.1 Contributions to Muskhelishvili Omnes inhomogeneities 

In Figs. [51 El and [7] we show the different contributions to the MO inhomogeneities A^(i) exemplarily 
for A*]_ for the un-, once-, and twice-subtracted case, respectively, from the tttt threshold t^r up to 
1.15 GeV, also indicating the upper limit of the KH80 solution as well as the KK threshold tx as 
the uppermost limit of approximate elasticity for J = 0. We choose the S-wscve for the following 
reasons: for J = the nucleon pole term is zero at t-,^ and does not dominate all other contributions 
as it does for the higher partial waves; in addition, for the S-wave we can also show the coupling of 
the D-wave as the leading example for the coupling of higher partial waves. The pole term N^ is 
independent of both the number of subtractions and a and thus serves as reference in all three plots 
(double-dashed). The s-channel contributions are shown separately for the sum of all partial waves 
with ^ < 4 in the range W £ [W+ , Wg] (dashed) and the Regge contributions of the full invariant 
amplitudes for W > Wg, (dot-dashed). Even in the unsubtracted case both the s-channel Regge as 
well as the t-channel D-wave contributions (dotted) are very small and almost negligible in comparison 
to the other parts. From this it is also clear that the coupling of higher t-channel partial waves (e.g. 
F-wave contributions to P-waves) can be completely omitted. The solid line denotes the sum of all 
these contributions and we have checked for J S {0, 1, 2} and n G {0, 1, 2} that the expected threshold 
behavior according to p.54p and (|3.58p (as for the corresponding partial waves) is indeed fulfilled. 
While all results are given for the optimal value of a unless stated otherwise, for comparison we also 
show the non- Regge s-channel contributions in the "fixed-t limit" a — ?■ — oo. Since this is a very drastic 
alteration (the RS system is not strictly valid in this case as will be explained below), the difference 
of this contributions for the two a values gives a very ample bound on the dependence on a. While 
the Regge contributions vanish for a — )• — oo as discussed in Appendix ID.ll the D-wa.ve coupling is 
not even well defined for a — t- — oo in this framework as can be seen from the explicit a-dependence in 
the unsubtracted case leading to an infinite contribution, cf. (j4.17p . By comparing the three plots it 
is clearly seen that in the once-subtracted case all contributions except the pole term are suppressed. 
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while in the twice-subtracted case an additional t-dependence is introduced such that they are strongly 
suppressed at tir, but at least the s-channel partial- wave contributions are comparable to the pole term 
around 0.75 GeV. For small t, the differences between the two a values are also suppressed by each 
subtraction as expected. 
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Figure 8: MO solutions for the S-wave. 



5.4.2 Muskhelishvili Omnes solutions: comparison with KH80 

We will compare our un-, once-, and twice subtracted MO solutions for |/^(t)| with J € {0, 1,2} for 
t G [t7r,im] with the KH80 results given as tables in [29]. Note that for J > 2 the un- and once- 
subtracted solutions coincide. The a-dependence (which is fully contained in A;[) can be used as a 
crude measure for the systematic uncertainties due to neglecting t-channel input above im (i-e. "non- 
analytic" input), since the physical result must be independent of a. Thus, for the five lowest i-channel 
partial waves we show our "KH80 consistency MO solution" for the un-, once-, and twice-subtracted 
case, each for both the optimal value of a and a — )• — oo in Figs. [8l [9l and [TOj Here, we have chosen 
to use the same value of ^/t^ = 0.98 GeV for all considered partial waves, which in principle is not 
necessary (the effect of varying t^ will be considered explicitly below). As discussed in Sects. 15.3.11 
and 15.3.41 this choice is mainly motivated by the 5-wave phase, which is just below vr at this energy 
(reaching vr around the KK threshold y/tx = 2Mk = 0.987 GeV) such that no additional subtractions 
are necessary in the MO scheme. In general, neglecting any input above the matching point enforces 
|/^(im)| = on the MO solutions. Nevertheless, even for the 5-wave we expect reasonable agreement 
with KH80 for this choice of tm, since both KH80 and [56] suggest that the modulus |/^| has a 
minimum or even an approximate zero between 0.9 GeV and ^/tx. 

In general, the solutions are fixed on both ends of the solution interval [tn,tm\- on the left due 
to the pole term and on the right due to the input above tm- In our case the solutions are forced to 
go to zero at the matching point since the input above tm is set to zero. With increasing J the pole 
term becomes larger and thus more dominating. Therefore, the differences between the n-subtracted 
solutions and also the different a values decrease close to t-^. As they furthermore agree very well with 
the KH80 solution in the respective pole-term-dominated regions, we only show the remaining regions. 
Since the D-wave coupling for the unsubtracted case depends on a, this contribution is omitted for 
the a — )• — oo limit (and thus the solutions for the two a values do not coincide at ^7^)- Obviously, the 
occurrence of a negative modulus (i.e. the unsubtracted |/? | for optimal a and the once-subtracted |/i| 
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Figure 9: MO solutions for the P-waves. 

for a — >■ — oo) only indicates that too much input information is missing in this particular case in order 
to yield a reasonable solution — a problem that can be cured by subtractions. The general pattern is 
as expected: the effect of varying a is suppressed by both the subtraction procedure and higher J. 
Furthermore, the agreement with the KH80 solution is strongly aided by subtracting. This is clear 
since each subtraction power on the one hand suppresses the lacking input above tm and on the other 
hand introduces additional consistent information via the subthreshold parameters as subtraction 
constants. Hence, the twice-subtracted solution for optimal a is our central "consistency result". The 
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Figure 10: MO solutions for the D-waves. For J > 2 one subtraction has no effect. 

S-wave shows a nice convergence behavior in n, but around 0.8 GeV it starts to deviate from KH80, 
which is not surprising as the /o(980) is expected to have an important impact (cf. Sect. I5.3.4p . As 
far as the P- waves are concerned, the numerical results confirm the analytic expectation that |/^| is 
much less well determined or constrained than |/1|: basically, the MO equations for \f^\ effectively 
contain one low-energy subtraction less. Moreover, in the necessary intermediate step of solving the 
MO problem for |r^| the pole-term contributions N^ cancel at tjr (as discussed in Appendix IB. 5p as 
for the S-wave, and thus the solution for \f}^\ is less pole-term dominated and hence more sensitive to 
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the values of the subthreshold parameters. Furthermore, the uncertainties of |/1| propagate into |/^| 
when calculating the latter from |r^| via (j5.20p . All this leads to a rather slow convergence behavior 
in n for fixed a as well as the loss of the expected convergence pattern in n of the differences between 
the two a values (note especially the crossing of the unsubtracted solutions for different a values and 
the negative once-subtracted modulus). However, our central twice-subtracted result for |/^| agrees 
rather well with KH80 especially in the p{770) peak, even though our result for |/1| (which enters |/+|) 
seems to be systematically smaller than KH80. Since this underestimation might be due to forcing 
the solution to go to zero at tm, we will investigate the effect of using a higher value for tm below. 
Nevertheless, the a — )• — oo variant of the twice-subtracted solution for |/1| agrees well with KH80 
in the p peak (though the agreement with KH80 becomes worse for |/1| in this limit). The D-wave 
results are systematically smaller than KH80. The change from one (or equivalently zero) to two 
subtractions towards KH80 is roughly one third of this discrepancy and furthermore approximately of 
the same absolute size for both partial waves, which is probably due to calculating |/^| by using the 
result for |/^| together with the fact that Xr = for all J 7^ 1. For both |/^| and |/^| the accordance 
with KH80 (which is based on fixed-f dispersion relations) in the "fixed-i limit" a — )• —00 is striking, 
the effect of varying a being much larger than the effect of subtractions. 

5.4.3 Muskhelishvili Omnes solutions: variation of the matching point 

Up to now we have used the S'-wave-motivated value \/t^ = 0.98 GeV for all considered partial waves. 
The effect of changing ^/t^ to e.g. 1.1 GeV is shown in Figs. [11] and [12] for the un-, once-, and twice- 
subtracted solutions for J G {1,2}. Again, for J > 2 the un- and once-subtracted solutions coincide. 
For J = 1 it is generally assumed that 4tt contributions can safely be neglected up to the vrw threshold 
around 0.92 GeV; however, the vrvr scattering P-wave inelasticity is small even above that energy and 
hence the impact of neglecting it (for both values of t^) should be smaller than the effect of changing 
im- For J = 2 we do not expect sizable deviations from elasticity, since the tttt D-wave is essentially 
elastic in this energy range. The P-wave solutions exhibit the expected behavior: the differences 
between the two matching-point values become smaller with each subtraction, but the convergence 
behavior in n is again less good for |/+|, where a higher value of t^ does not lead to a better agreement 
with KH80, while for |/1| already one subtraction in combination with the higher matching point yields 
a description of the KH80 solution that is even better than the twice-subtracted version for a — )• —00 
discussed before. Therefore we conclude that on the one hand the KH80 solution for |/1| can be 
reproduced well with a higher matching point already in the once-subtracted case, but on the other 
hand the KH80 solution for |/^| calls for a second subtraction and is hard to be accommodated in our 
MO scheme for energies above roughly 0.8 GeV. The D-wave solutions, however, are hardly affected 
at all in the KH80 energy range by changing t^. As discussed in Sect. 15.21 they are expected to be 
dominated by the pole terms iV^, which for comparison are also shown in Fig. [12] While for |/^| the 
KH80 solution indeed agrees rather well with the pole term itself throughout the whole KH80 energy 
range, for |/^| there are sizable (with respect to the scale) deviations between KH80 and the pole 
term in this region, which again fits the picture that the partial wave with parallel helicity is both 
analytically and numerically less well constrained. Together with Fig. [10] we can conclude that in the 
limit a — )• —00 for |/^| the net effect of adding the dispersive integrals to the pole term is very small, 
while for | /^ | the corresponding dispersive contributions (which thus are not mainly induced by | /^ | 
in this limit) are crucial for the agreement with KH80. For optimal a (and independent from the 
choice of tm)) though, these contributions deteriorate the agreement with KH80 (with respect to the 
pole term) for |/^|, whereas improving the agreement for |/+|; in this case the corrections to |/^| are 
effectively due to \fl\. 
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Figure 11: MO solutions for the P- waves with V^ =1.1 GeV. 



5.4.4 Application to nucleon form factors 

The t-channel partial waves considered in the previous sections are not only an integral part of any 
closed system of dispersion relations for vrA^ scattering fully consistent with crossing symmetry, but 
also an essential ingredient to dispersive analyses of nucleon form factors. The contributions to the 



52 



100 



80 



2M^ 



•^ 60 

O 






=+-= 



40 



20 



I 

> 

O 



-§ 



=^ 



100 



80 



60 



40 



20 



0.3 

2M^ 



40M^ 0.98 1.1 

0-sub(VC'= 1.1 GeV) '- 

0-sub 

2-sub(v^= 1.1 GeV) 

2-sub 

KH80 X 
pole term 




0.4 0.5 



0.6 0.7 0.8 0.9 1.0 1.1 

Vi/GeV 

VlOAf^ 0.98 1.1 

"^ 0-sub(VC'= 1.1 GeV) '- 

0-sub 

2-sub(\/i;;; = 1.1 GeV) 

2-sub 

KH80 X 
pole term 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 

^/t/GeV 

Figure 12: MO solutions for the D-waves with ^/ti^ = 1.1 GeV. For J > 2 one subtraction has no 
effect. 



isovector spectral functions by two-pion intermediate stateqiO in the case of the electromagnetic Sachs 
form factors read [55] (cf. |29y90) for precise definitions and [91j for a recent application) 



ImG^it) 



Tt 



Qt 



m 



Vt 



{F^{t)yfi{t)e{t-U) , ImGMt) = -%{F^(t)yt{t)e{t-t^) , (5.30) 



^^G-parity dictates that intermediate states of an even (odd) number of pions only contribute to the isovector (isoscalar) 
spectral functions; cf. Sect. 12.21 
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Figure 13: Two-pion-continuum contribution to ImG^(i) and Im G^{t). 

while the imaginary part of the scalar form factor is determined by [92j 

3gt 



lma{t) 



^PtVt 



UF^(t))*fl{t)e{t-t^), 



(5.31) 



with the scalar and vector pion form factor F^(t) and F^ (t), respectively. In the case of the scalar form 
factors the approximation by tttt intermediate states breaks down as soon as the two-kaon threshold 
opens, and effects from KK intermediate states are known to be important for a dispersive description 
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of F^{t) |531[93]. In contrast, the two-pion contribution dominates in the vector channel, where 
inelasticities set in more smoothly. It is thus instructive to investigate the impact of our results for 
|/j_(t)| on the spectral functions of the Sachs form factors. To illustrate the corresponding effects we 
approximate the vector pion form factor by a simple twice-subtracted Omnes representatioro (cf. [95] 
and references therein) 

using (r^)]|f = 0.435 fm . The results for the once- and twice-subtracted versions of |/j_| together with 
the comparison to KH80 are depicted in Fig. [I3l As expected from the discussion in Sects. 15.4.21 
and 15.4.51 the relative uncertainty in ImG^ is much larger than in ImG^, which is a result of the 
effectively lower number of subtractions in the calculation of |/^| and its enhanced subthreshold- 
parameter dependence. However, since ImG^'^ is much larger than ImG^, the absolute deviations 
between the individual curves are actually of comparable size in both cases. We conclude that a 
new determination of the subthreshold parameters from a full solution of our RS system should lead 
to improved central values and associated uncertainties for the two-pion contribution to the spectral 
functions of both form factors. 

6 Conclusion 

In this article we have derived a closed system of Roy-Steiner equations for ttN scattering and ana- 
lytically calculated the kernel functions for the lowest s- and t-channel partial waves. Furthermore, 
we have constructed the corresponding unitarity relations in detail, including inelastic contributions 
from KK intermediate states in the t-channel reaction. To pin down the optimal value of the free 
hyperbola parameter a, we have analyzed the domain of validity of the full system (assuming Man- 
delstam analyticity) and determined a for both the s- and t-channel equations such that the range 
of convergence is maximized. We have introduced subtractions at the subthreshold point in order to 
suppress the dependence on the high-energy region and derived the corresponding once- and twice- 
subtracted versions of our Roy-Steiner system as well as sum rules for the subtraction constants and 
the necessary corrections to the kernel functions. 

Casting the t-channel equations into the form of a Muskhelishvili-Omnes problem with finite 
matching point, we have then solved the equations for the t-channel numerically in the single-channel 
approximation. We have assessed the numerical importance of different input contributions for the 
Muskhelishvili-Omnes problem and its solutions and examined the behavior of the Muskhelishvili- 
Omnes solutions for the lowest t-channel partial waves (J G {0,1,2}) with respect to varying both 
the input and/or the framework in many ways, including their sensitivity to the vrvr phase shifts, the 
number of subtractions (n S {0, 1,2}), variation of the matching point tm, and taking the hyperbola 
parameter a to — oo. In general, we find consistency with the KH80 solutions. However, our analysis 
shows that the structure of the equations renders the t-channel partial waves /^ systematically less 
well determined than their counterparts f:[ both due to an enhanced sensitivity to the subtraction 
constants and an effectively lower number of subtractions. Finally, we have briefly discussed some 



■^^This representation ensures that F^ fulfills Watson's theorem, such that the phases in H5.30p cancel. Strictly speaking, 
using any representation that goes beyond the two-pion approximation would be inconsistent unless the corresponding 
inelasticities are accounted for in the determination of f± and the unitarity relation (|5.30|) as well, as exemplified by 
the breakdown of Watson's theorem and the spectral functions' becoming complex. Moreover, the precise value of {r'^)X 
is immaterial in the present context, since we merely wish to convey how the uncertainties in \f±\ propagate into the 
spectral functions. The present choice ensures a decent description of form-factor data, cf. 94 . 
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consequences for nucleon form factors, in particular our analysis gives a first indication where the 
largest uncertainties in the spectral functions are to be expected. 

The next step in the solution of our system of Roy-Steiner equations will be the incorporation 
of KK intermediate states in a full two-channel Muskhelishvili-Omnes treatment of the i-channel 
S-wave, which will have immediate consequences for the extraction of the vrA^ a term via the scalar 
form factor of the nucleon [82) . Having then solved the t-channel part of the system, the s-channel 
equations are solvable with techniques similar to those employed in the context of vrvr Roy equations, 
and finally the iteration of the full system should determine the lowest partial waves as well as the 
subtraction parameters. We are confident that the framework proposed in this article will allow for a 
reliable extrapolation to the Cheng-Dashen point and thus for an accurate determination of the vrA^ 
a term. 
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A Partial-wave projection for the s-channel amplitudes 

In this Appendix, the different contributions to the s-channel part (j3.4p of the RS system are discussed. 



A.l Nucleon exchange 

The useful general definitions 
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for the cases where the crossing-symmetry constraint applies (i.e. J even/odd for / = +/—). Projecting 
the HDR Born terms N\s,t) of ()2.48p onto s-channel partial waves via (j3.2p leads to MacDowell- 
symmetric nucleon pole contributions 
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which, by defining the abbreviatioro 



NLiW) 



9 



IQttW 



{E + m){W -m) 



iQijy) . 26io 



m^ — s 



+ {E-m)(W + m) 



l Ql±l{y) 2(5;-fci,o 

o ~r r, 

(A.4) 



for later convenience, can also be written in the form 
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and where we have defined 



47r 



2W 



m^ — a 



g^ {E -m){W + m) 5ii 
4^ W/ m?-a 



V/>0 
V/> 1 



y{s) = 1 



s + m 



2q^ 



Zs{s,t{s,u = m )) = Xs{s, s = m ) 



(A.5) 



(A.6) 



{xs{s, s') will be introduced in (jA.2ip ). Qi{z) denotes the Legendre functions of the second kind. 

The Qi{z) obey a recursion relation similar to the one for the usual Legendre polynomials Pi{z) 
(for / > 0) 



(/ + l)Pi+i{z) + IPi-iiz) = {21 + l)zPi{z) , 
{I + l)Qi+i{z) + lQi-i{z) = {21 + l)zQi{z) - 5io , 



(A.7) 



which, together with Qi = Pi = foi I < 0, leads in particular to (cf. (IA.56P for the general formula) 

Qi{z) = Pi{z)Qo{z) - 1 , Q2{z) = P2{z)Qo{z) - \z , Q^{z) = Pi{z)Q^{z) -^- + \■ (A.8) 
From the Neumann integral representation for general complex argument z [96] 
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one can read off the lowest function for general real argument y 
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We also need the analytic continuation for purely imaginary argument z = iy, e.g. for y > 1 



Qo(iV) = - log = - log I— = I { arctan-u 
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(A.ll) 



Functions with I > 1 may then be obtained via either the recursion relation (IA.7P or the reduction 
formula ()A.56p . 



^^Of course, also this form of the nucleon pole terms obeys the MacDowell symmetry relation (|3.3p . since the term 
proportional to (5!+i.o/(5;o vanishes for iV;^ as a consequence of I starting at 0/1, respectively. 
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A. 2 s- and u-channel exchange 

By introducing a convenient matrix notation via 

the crossing properties of the even and odd invariant amphtude combinations (j2.18p read 

A^(z., t) = e'as A\-i^, t) , ^3 = (l ^^ ■ (A.13) 

While the s-channel partial-wave projection ()3.2p can be rewritten as 

1 
fi\W)= [dzsK'iW,Zs)A'is,t) (A.14) 

J t=t(S,Zs) 

-1 

where the projection kernel matrix is given by 

R'(^,-.)=fS'+^ !!H, (A.15) 

\'"'i+i,i ^1+1,1/ 

Rl^{W,Zs) = j^{{E + m)Pk{zs) - {E - m)Pn{zs)} = -RU-W,Zs) , 

Rl^{W,z,) = r^{{E + m)iW-m)Pkizs) + iE-m)iW + m)Pnizs)]=-Rlk{-W,Zs), 

the s-channel partial- wave expansion, i.e. the inversion of ()A.14p . takes the form |97] 

00 
A'(^>0 , ,, =Y.^'{W,ZsW{W) , (A.16) 

with the expansion kernel matrix 

s\w,zs)=(%+''^ :$'+^ 

SUW,Zs) = 4vr{^±^P^(z,) + ^^^p;^{z,)\ = -S^^^^-W,zs) , 
I h + m t, — m J 

SUW,Zs) = 4vr j-^ Plizs) - —^ P^(z.)| = -Sl^{-W,zs) . (A.17) 

l£j + m hi — m J 



In accordance with the matrix form of the MacDowell symmetry relation (j3.3p 

il{W) = -a,il{-W), ^i=(j J), (A.18) 

these kernels obey the symmetry relations 

R'(Ty, z,) = -a{R}{-W, zs) , S\W, z,) = -S\-W, z,)cri . (A.19) 
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With the definitions (j2.5ip the s- and li-channel terms of the HDRs (|2.44p thus can be cast into the 
matrix form 



A'{s,t) 



s+u 



1 



Ir' ^'^ 



t=t{s,Zs)=-2q'^{l-Zs) IT 

where the HDR kernel matrix h^ is given by 

H(.s,s';zs) = hiao-e^h2(T3 , (Tq = I2 
1 1 1 



— / ds' hg[s,s';Zs]linA {s',t 



t'=t' is', z'^) = -2q'Hl-z'^) 



(A.20) 



hi{s,s') = — 

C' 



s' — s 2 s' — a 



h2{s,s;zs) 



1 0\ 
l) ' 

1 1 



Xs{s,s) = 1 

1 1 



s + s' -T. 



2q^ Xg — Zs 2 s' — a 



(A.21) 



and [s,s';2;s] indicates that the whole integrand is to be understood as a function of these variables, 
which can be achieved using 



Zg{s,s';zs) = azs + /3 , a{s,s') 



s — a 



q'"^ s' — a 



P{s.s') = l-a- '\ ;^+;l ^ . (A.22) 



s' -a 2g'2 

By expanding the absorptive part of the s- and u-channel HDR terms given in (IA.20p into s-channel 
partial waves via ()A.16P and projecting out s-channel partial waves again by means of ()A.14p . we 
arrive at the partial-wave dispersion relations 



00 ,„ 

f/(VF)r " = - dW ^K"''^(t^,VF')Imf/(H^'), 



(A.23) 



w+ 

where the s- and ti-channel kernel matrix is defined by 

1 



K"''^(VF,Ty') = 2W' f dzs B}{W,Zs)hi[W,W';zs]S^'iW',z'g) . (A.24) 



-1 



Due to the symmetry relations 



criK 



II' J t 



W, W') = K" '^(VF, W') = K" '^{W, -W')(Ti , 



(A.25) 



which follow from the relations (]A.19p . the s- and u-channel kernel matrix can be written with only 
one kernel function according to 



-Kl,i-W,W') -Kl,i-W,-W')J' 



K[i,{W,W') = 2W' j dzs \B}{W,Zs)\ii[W,W'-Zs]S^'{W',z'g)^^^ , (A.26) 



where the subscript denotes the 1, 1-th element of the matrix in the brackets. The PWDRs (1A.23P 
then take the form already stated in (|3.4p 



fLm 



s-\-u X 



1 f ^^ 

''w+ ''=0 



-/(/+!)- (-^) 



S+M 



(A.27) 
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Defining the structure 

^[akn\h{W,W')\ = ^{b{W,-W')akn + h{W,W')ak,n+i 



w 



+ h{-W, -W')ak+i,n + h{-W, W')ak+i,n+i] , (A.28) 

where afc„(s,s') is to be understood as a function invariant under sign changes in W and W, and 
introducing the kinematical abbreviations 

E'{W') = E{W') , HW,W') = ^;f^[W' + W], g{W, W) = ^^^ "^ [W -W + 2m] , 



1 



>i^{W, W) = - 5{W, W) + e^Q{W, W) 
as well as the angular kernels 



E' + m 

E + m 



E' + m 



E' + m^ 
e+(VF' + m) +e_(VF-m) 



1 ± 



Pl{Zs)P'v{z's) 



Xq Zq 



(A.29) 



(A.30) 



-1 
the general s- and u-channel kernel function can be written as 



Kli,{W,W') = hip[Ui,\d{W,W')] - -cp 



Vw ^ Uw 



2 ' s' — a 



q{W, W) 



^[Uiv\5{W,W')] iip[Viv\Q{W,W')] ^[Uiv\^\W,W') 



s — s 



2^2 



s — a 



(A.31) 



Since c^[afc„|6(Wi PF')] encodes the MacDowell symmetry (j3.3p for both pairs {k,W) and (n, VF'), we 
can decompose it in two ways 

^[akn\h{W,W')\ = ^i[akn\h{W,W')\ - ^i[ak+i,n\h{-W,W')\ , 
= ip2[akn\b{W,W')] -ip2[ak,n+i\b{W^-W')] ' 

ipi[akn\biW,W')] = ^{b{W,-W')akn + b{W,W')ak,n+i} , 

^2 [akn\b{W, W')] = ^{b{W, -W')akn + b{-W, -T^')«fc+i,n} , (A.32) 

and with the definitions 

the kernels exhibit the following interrelations 

Kli,{w, W) = kI^^{w, W) - kI^\ ,{-w, W) = kI^?{w, W) - Klf_^^{W, -W) 



(A.34) 



that may be used to write down explicit expressions of the kernels in a compact form. However, 
for numerical evaluations a different prescription is preferable. The part of ()A.3ip that contains the 
s-channel cut can be decomposed according to 



ip[Uu>\5{W,W')] _^n'{W,W') 



+ 



1 W { E + m^^ E -m 



s — s 



W -W W' + W W \E' -m " E' + m 



U, 



w 



-U, 



i+i,i'+i (■ ) 



,,,. ...,, W ( E + m E-m 

JU'{W,W) = -— \ Ui^v+i - — Ui+i^v 

W \ E' + m E' — m 



(A.35) 
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Using the identity 



1 
ldzPl,{z)[Pi^,{z)-zPi{z)]=^S^-^\^ , (A.36) 

-1 



we can easily calculate its residue at the pole W' = W (where a = 1, /3 = 0, and thus z'g = Zg 



Res 



ln'{W,W')^,^^ 



iw {W, W) = Ui,i,+i {s,s)- Ui+i^u {s,s) = 6w , (A.37) 



w -w 

which together with the decompositions 

Uiv{W, W) = Uu'{W, W) + {W - W)Uiv{W, W) , 

WE'±m ^^ '^' c±^vv,vv ) 2W'{E'±m)s ' ^ ^ ^ 

leads us to the alternative form of the kernels Kii,(W, W') 

E-m^^ \ j ip[Vii'\g{W,W')] ^[Uu>\>c'iW,W')] 

-'Jl+i,l'+i ^ - e TT^ ; , (A.3yj 



E' + m '-^^-'-^^J 2q 

where the first term is the usual Cauchy kernel for the s-channel cut (contributing only for 1 = 1') and 
the kernels Kli,{W,W') contain only the left-hand cut. In order to derive explicit expressions for the 
angular kernels UwiW = ^/s, W = vV) from ()A.30P and subsequently their regular parts Uiii{W, W') 
from ()A.38p . we use the following expansion [96j 

P,iz.) = ±aW, ai = t^^j^, .= '—^, (A.40) 



imi-xy. ' 



and hence 

I'-i 

2 



I'—i 
P['{z's) = -lY.(^' + l)ai,^,x"' , x' = i-^=u; + ax, c.(s,s') = ^—^±^ , (A.41) 



A'=0 

together with the binomial theorem and the Saalschiitz identity [96 

/x + A + 1 ^ ^ in - l)l{fi + I + 1 



E ^^ = (-')' , „_Z.,.n. t"^')- (A-^2) 



A=0 

to arrive at the general expression for the angular kernel Un 



Uii'{s,s') 



A — L U, — ( 



These kernels show the following asymptotic behavior: 

Uw ~ q^^ for q^O, Uw ~ g'-2^'+2 for g' ^ , Uw ~ q~^^ for q ^ oo , (A.44) 
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and, in particular, the lowest kernels are given by (note that Uio = and Un = 6io) 

Uiv = for I' <l, Uij+i = a^ , Ui,i+2 = {21 + 3)(3a^ , 

a' 



Ui,i+3 = y { (2/ + 5) [a^ + (21 + 3)/32] _ (2/ + 3) } . 



From (|A.43p and (|A.30P we can easily deduce for W = W 






^w 



for I' < I OT I' — I even , 
for /' — / odd , 

(-l)'+^'(/' + A' + l)! 



2(A' + 1)(/' - A' - 1)!(A' - /)!(/ + A' + 1)! 
which again yields Ui^i'-^-i{W,W) — Ui+i^i'{W,W) = dw. By defining 

"E2 



-,r2 

s — a 



l^ = l + (W'- W)di , di{W, W) = '' ^'' 



1 + {W'-W)d2, d2{W,W') 



4g/2 

w + w 



1 



ss' 



s — a s — a 

we can rewrite the powers of a according to (note that 6o = 0) 



fc=0 



s — a 



s — a 



+ 4+1 



which together with the definitions 



uj = {W' - W)u) , u}iW,W') 



W' + Ws + s' -i: 



Aq 



12 



uu'iw, w) = ^!^ 'x;V' + 1)4+1 E ft') 



s — a 



(a. -/)!(/. + / + 1)! 



^A'-l-|.^M 



allows us to give the explicit form of the regular part Uu/ of the angular kernel Uw as 

/'-I 



Uu>{W,W') = Y,Uu'by + OUw 



X'=l 



from which we can easily obtain the lowest Uu' (note that Uiq = = Uu) 

till' =0 for r < / , Ui,i+i = bi , Ui^i+2 = -{21 + 3){6/+i - 6/ + 2a)a'} 



i7z,z+3 = (2/ + 5)(2/ + 3)<^(/ + 2 



bi+2 h 



2l + 2> 2/ + 5 



5^+1 



uja\l - a + P)) . 



(A.45) 



(A.46) 



(A.47) 



(A.48) 



(A.49) 



(A.50) 



(A.51) 



The angular kernels Vw can be expressed by the kernels Uw as follows: from the integral representation 
of Uw ()A.30P we can deduce that 



/'-I 



Pv{z's) = Y.{2n + l)UnvPn{zs), 



(A.52) 



n=0 
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and inserting this into the integral representation of Vw ()A.30P yields 



z'-i I ^ 

n=0 [ _^i 



p„(z,)PK^.) 



Xe Ze 



> . 



(A.53) 



By using the identity 



1 /dx ^iMf^ = P^(z)Qi{z) forn < Z , 

2 ./ z — X 



(A.54) 



we can write 



I'-i 

Vw = Ql(,Xs)P[,{x'J - Yl {2n+l)Unl'[Pn{Xs)Ql(.Xs)-Pl{Xs)Qn{Xs)} , 

n=l+l 

s' + S-T. 



Xg{s, s ) = axs + /3 = 1 



2q' 



/2 



Xs{s',s) 



(A.55) 



(note that the sum vanishes for I' < I + 1), and with the aid of Wi^i{z), which is a polynomial of 
degree / — 1 in z defined by [96] 

Qi{z) = Qo{z)Pi{z)-Wi.i{z), W_i = 0, (A.56) 

leading to the integral representation 



Wi-i{z) = - dx 

2 / z — X 



(A.57) 



the angular kernels Vw take the general form 

Vii'{s,s') = Qi{xs)P'i,{x'g) - Viv , 

Vw{s,s')= Y, {'in+l)Ur.v[Pl{Xs)Wn-l{Xs)-Pn{Xs)Wl-i{Xs)] . 

n=l+l 

The Vw only contribute for I' >l + 2 

Vw = G for /' < / + 1 , 

and we can immediately read off 

I'-i 
Vl0 = Q, Vli = Ql{Xs), Vov=Qo{xs)Pv{x,)-Y,{'^i^ + ^)UnvWn-i{xs) 



n=l 



(A.58) 



(A.59) 



(A.60) 



where the second equation can also be seen directly by comparing ()A.30p with ()A.9p . Furthermore, 
one easily obtains the asymptotic behavior 



Vw 



^2/+2 



for q^O , Vw ~ ^'-2«'+2 f^j. q' ^0 , Vw ~ g'"^'"^ for q ^ oo . (A.61) 
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From ()A.57p or from Christoffel's formula for / > 1 [96l 



l-l 



1 - 1 for / > 2 even , , , 

2 - ' A.62 
^ for / > 1 odd , 



both yielding (also in agreement with (jA.SP and ()A.56P ) besides W-i = 



Wo = l, Wi{z) = -z 



W2{z) = ^-Z^ 



5 2 2 

-•z , 

2 3 



(A.63) 



where it is useful to note that Wi{z) like Pi{z) contains only even/odd powers of z for I even/odd, 
respectively, we can immediately deduce the non- vanishing angular kernels Uw, Uw, and Vw for /' < 3 

1 



Ui3 = a^5i2 + 5af36n + -{5[a^ + 3/3^] - 3)5^0 , 



Uii = (5/0 , Ui2 = a6ii + 3/36io , 

Ui2 = bi5ii - 3{fei + 2uj}5io , 

Ui3 = b26i2 + 5{6i - 62 - 20a}6ii - 5{36i - 262 + 3w(l - a + /3)},5zo , (A.64) 

yii = Qi{xs), Vi2 = 3x'^Qi{xs) - 3a5io , 



5 15 



2 2 

that are needed for the kernels K^/iW, W) for all combinations (/ > 0, /' < 2) according to ()A.39p 



kUw,w') 



1 W >i^(W,W') , ^ 

+ C+--F77^ -}^10 



W'-W 

ei W' 

2^W 



W s' -a 
q{W, W')Qi{xs) + q{-W, W')Qi+i{xs) \ , 



1 W f E + m 



a^:^^W5io-^^[-2,aQ{W,W')5iQ+\3x',Q{W,W') + Q{W,-W') Qi{xs 



E' + m 



2g2 W 



+ 



3x',g{-W, W) + g{-W, -W')\ Qi+i{x, 

1 w 



■[a>i\W, W')6ii + \3l3x\W, W) + x\W, -W) + ax\-W, W')] 5io\ , 
s — a W ^ L i ) 



Q 



4^^1 ]^ii + I - 2(8fci - 562) - 15w(l - a + /3) + i(5[a2 + S/S^] - 3)c+ 
l!/' + m / I I 2 



ac_ + 777;— -777-777 I 3-77^ Sottt- I }m 



I TI/' 



W' + W W \ E' -m E' + m 



a 



15{|x, + (3}g{W, W) + 3^(VF, -VF') + ^aQ{-W, W')] Sio 



e W 
- ^a^QiW, W')6a + [p!,{x',) q{W, W) + 3x',0{W, -W) \ Qi{xs) 

p:,{x',)q{-w, W) + 3x'^e{-w, -w')] Qi+i{x,) 
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1 w 

s' -aW 
1 



a^x\W, W')5i2 + a \bPx\W, W) + i<\W, -W) + ax\-W, W) 



{b[a^ + 2,f] - 3)>i^{W, W) + 3f3>€^{W, -W) 



+ 5aPx'{-W, W) + ax\-W, -W) 5io 



6ii 



(A.65) 



From KL(W, W) in the form according to ()A.3ip 



kUw,w') 



1 w 

2W E' + m' 



{E + m) 



iW + H^')4-^ + ^\W -W - 2m) 



Ql{Xs) 



S' — S 



(A.66) 



+ e\E -m){W + W' + 2m) ^'+^j^'^ - {E + m) €+{W' + m) + l_{W - m) 



2(5zo 



s' — a 



we can deduce that the nucleon pole terms (IA.3P are reproduced by 

Nl+{W) = -fKloiW, -W = m) = -Nl,^^^^{-W) V / > 



(A.67) 



The expUcit formulae for the additional non- vanishing angular kernels C/;;/, Uui, and Vui for (/ < 2, 4 < 
I' < 6) needed for calculating the additional higher kernels K^, for (Z < 1, 3 < I' < 5) via ()A.39P 
are displayed in Appendix IA.4I Furthermore, we give the asymptotic behavior of the general kernel 
function K^,{W, W), which can be inferred from the asymptotic behavior of the angular kernels (|A.44p 
and fOT]) . 



for g — )• 
for g' — )• 
for g — )■ oo 



KUW,W') ^ q'-^'-^ 



KU-W,W') ^ q^'+\ 
Ki,{W,-W')^q-^''-^ 



{km) 



in agreement with the MacDowell symmetry relation (j3.3p . From ()A.67p we can then read off the 
asymptotic behavior of the nucleon pole terms 



iV/+(VF)~iV(^+i)_(-Ty)~g^' forg^O. 



(A.69) 



A. 3 t-channel exchange 

With definitions (|2.5ip and relations (j2.52p the t-channel terms of the HDRs (j2.44p can be written as 



A^(s,t) 



* _ 1 

t=t{s,zs) vr 



It I 4.1^ 



dt' hj[s,t ;Zs]ImA^(s',t 



,'=s'{t',z[) 



(A.70) 



where the HDR kernel matrix h^ is given by 



hl{s,t';zs) = 
Xi{s,t';zs) = 




2g2 Xt- Zs V 
f '^\ 2 



^0 a1)' xt{s,t') = l + ^ = zs{s,t') 



(with X = 1 V x) , 



(A.71) 
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and the integrand is to be understood as a function of [s, t'; Zg] by using 
// / N "T-^' / F / /^ ci'^(s — a) 



(t'-^ + 2af-A{s-a){2q^ + T.-s-a) 
^{s,t) = ^-^-^^ . (A.72) 

The t-channel partial- wave expansions of the invariant ampUtudes, i.e. the inversion of (j3.6p . read [49] 



where it is crucial that the sums only run over even J for / = + and odd J for I = — due to Bose 
symmetry. Taken literally, the form ()A.73p of the partial- wave expansions is only valid for t' > tN, 
since below the two-particle thresholds t^v and t-,^ the CMS momenta p[ of the nucleons and q^ of the 
pions become purely imaginary and one has to use p'_ and q'_ instead, respectively (cf. (j2.1ip and [29]). 
In particular, in the unphysical range t' £ [i7r,tAr) that we are interested in as the low-energy part of 
the integration range t' G [^,^00), we have g^ G M but Pf,Zf G iW. However, the squares pf and q'-^ 
are always real (albeit not necessarily positive, cf. (|2.12p ) and since the combination p^q'tz'^ = mv' = 
m{2s' + t' — S) is always real as well, so is z^^. Due to the fact that the Legendre polynomials and their 
derivatives have definite parity Pj{—z) = {—iyPj{z) and Pj{—z) = {—l)'^~^Pj{z), a closer look at 
the expansions ()A.73P shows that in all cases only powers of the real combinations Ptqt^t and additional 
factors of powers of the likewise real squares p'j^ and g^^ appear. Therefore, we can symbolically use 
these formulae for all kinematical ranges and factor out powers of the real squared momenta whenever 
necessary in order to form explicitly real quantities. 

By introducing the t-channel partial-wave amplitudes into the matrix notation vislfj 






(A.74) 



the expansions (|A.73P can be rewritten as 



A^(.',t')[,^^,^^,^,^ = X;T^(t',zDf^(t'), (A.75) 

where the expansion kernel matrix is given by 

TV,4) = o('o' I'^y CAt') = M'^J + l){pWty-\ (A.76) 

uj{t',z',) = -^PM) , vj{t\z[) = "l^ ^z'.P'M) , wj{t',z',) = ^L^ P'M) . 

As the sum only runs over even J for I = + and odd J for I = — and thus the full information on 
the crossing properties is already contained in the index J, we can redefine 

i+(-i)"+J 
Xi{s,t';zs) = xiis,t';zs) = [^) ' (with x° = 1 V x) , (A.77) 



^"in order to accommodate the fact that there is no /° to the matrix notation, we define /° = and in the following 
all corresponding quantities (e.g. integral kernels) are also understood to vanish. 
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and omit the index / in favor of J in the fohowing. If we expand the imaginary part of the t-channel 
HDR terms in (lA.TOp into i-channel partial waves via (|A.75j) and project out s-channel partial waves 
again by use of ()A.14p . we can obtain the following PWDRs 

oo 

¥l{W)\' = - fdt' Y,G^'^(.W,t')lm¥-^{t') , (A.78) 

where the t-channel kernel matrix is defined by 

1 
G''^{W,t') = fdzs Ii\W,z,)hl[s,t';Zs]T-\t',z't) . (A.79) 

-1 

Due to the symmetry relation 

G^'^{-W,t') = -aiG^'^{W,t') , (A.80) 

which follows from ()A.19P and is in accordance with the MacDowell symmetry (jA.lSp . the f-channel 
kernel matrix can be expressed by two kernel functions 

where in accordance with /£ = for the matrix notation we set Hiq = 0, and the PWDRs (|A.78p take 
the form already given in ([37 



CX) 



//+Wf = iy"dt' Y.{Gu{W,t')lmfi{t') + Hu{W,t')lmfl{t')} 



t-TT 



With the definitions 



fLi)-i-W) . (A.82) 



2 7-1-1 fo'o'V 
'iP[akn\d{W)]=d{W)akn + d{-W)ak+i,n, VjiW,t') = ^^ ^ "-^'J^ , (A.83) 



and by introducing the angular kernels 

1 



2 J xt- Zs 2 J xt- Zs 

-1 -1 

1 

Cu{s,t') = \Uzs ^j^ Pi{zs)z[P'M) ^ j^^^ ^ ^^^^_^ ^ ^^_g^^ 

^ J Xt Zs 



-1 
we can write the kernel functions as 



Gu{W,t') = -r^jij[Au\E + m] VJ>0, 

Hij{W,t')= J^^ \^^[Bu\{W-m){E + m)\+mi^[ClJ\E + m]\ VJ>1. (A.85) 

\/JiJ + 1) IQt ) 



VJ(J + 1) [Q't 
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If we use the decomposition 
V 1 



, = — + — ) Mi(s,0 

v' xt - Zs z^ z^ xt - Zs 

we find for the angular kernels for even J 

1 



Mit 



l^2{s,t') 



2s + t' -T. 
Mq't 



Aijis,t') = ^jdz, 



Xt - Zs 



-1 



i i 



Pl{Zs)P'M)/z[ 



Xt - Zs 



(A.86) 



(A.87) 



and for odd J 



Au{s,t') = ^ I dzs Piiz, 
-1 



PAzi) ^[^ [^^^ Pl{Zs)PM)/z[ 



Xt - Zs 



-1 



Xt - Zs 



(Ai 



from which we can infer that only even powers of z[ occur and hence a square-root dependence on 
Zs is avoided. We now can work out the kernel functions explicitly, here given for all combinations 
(/ > 0,J<2) 



Gi^{W,t') = - \ ,^ [{E + m)Qi{xt) -{E- m)Qi+i{xt)] , 



AWq^Pt 



E + m 



Gn{W,t') = ^|(2s + t' - ^)Gio{W,t') + 
HniW,t') = l=!^^Zi{W,t')-mGii{W,t'] 

Hi2{W,t') = ^^|(2s + t'- ^)Zi{W,t') -m[4s{s + t' - S) + {t' - S)^] Gio{W,t') 



^E + m 



m{s — a) 



P? 



+ W — m 



^10 



where we have defined 

Zi{W, t') = ^[{E + m){W - m)Qi{xt) + {E - m){W + m)Qi+i{xt)] . 
From the expansion 



Ql{xt) 



'-fTMmi -i^^^ 



(A.89) 



(A.90) 



(A.91) 
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we can finally deduce the asymptotic behavior of the non- vanishing general kernel functions (jA.SSP 

for q^O Gu{W,t') ~ Hu{W,t') ~ q^' , Gij{-W,t') ~ Hu{-W,t') ~ q^'+^ , 

for q't ^ Gu{W, t') ~ Hu{W, t') ^ 1 , 

for p[ ^ Gu{W, t') ~ Hu{W, t') ~ p'f^ , 

for t'^oo Gij{W, t') ~ Hu{W, t') ~ t'-^-'-^ , (A.92) 



in accordance with the MacDowell symmetry relation p.3p . 



A. 4 Higher kernel functions 



Here, we display the explicit form of the additional angular kernels Uw, Uu', and Vw for (/ < 2,4 < 
I' < 6) that are required for calculating the additional higher kernels K^, for (/ < 1,3 < /' < 5) 
via ()A.39P needed to incorporate higher resonances in the s-channel integrals. From ()A.43P we obtain 



u„ = l 



/3|7a2 + 7/32-3} , f/14 = -aJTa^ + 35/3^ - sj , [734 = 70^/3 , (A.93) 

[/o5 = J|l5-70(a2 + 3/32)+63(a^ + 5/3^ + 10a2/32)| , C/15 = ^a/jjoa^ + 15/3^ - sj , 

f/25 = ^a^{9a^ + 63/32-7} , C/oe = y/jjs - 30(a2 + /J^) + ii(3a4 + 3^4 ^ -^q^2^2^| ^ 

C/16 = -a|35 - 126(a2 + 5/3^) + 33(3a^ + 35/3^ + 42a2/32) | , [/ag = -a^phsa^ + 77/3^ - 2l| , 

and ([X50]l yields 

C7o4 = -5|763 - 1462 + 96i + a;[4 - 14a(l - a) + 7/3(1 - a + /3)] | , 

i7i4 = 2I63 - 3562 + 156i - 35a)a(l - a + /3) , [724 = -7|63 - 62 + 2a)a2| , 

- r 15 ~i 

Uo5 = 7| I864 - 5(963 - 862 + 36i) - — cD(l - Q + /3) [(1 -3af + 3/3^] | , 

[7i5 = -7|l264 - 2763 + 2O62 - 56i + Oa [2(5 - 15q + 12a^) + 15/3(1 - a + /3)] } , 

[725 = 3664 - 6363 + 2862 - 630)0^(1 - a + /3) , 

i7o6 = -2l|22(65 - 364) + 5(1563 - 862 + 26i) + £Z'[2(1 - 9q + 31^^ - 22q3(2 - a)) 

+ j(l - a + /3)(3 - lla(6 - 13a) + 33/3^)] | , 

[7^6 = 66(565 - 1464) + 35(2763 - 1262 + 26i) - —Oa{l - a + /3) [(5 - IIq)^ + 55/3^] , 

[726 = -3/11(565 - 1264) + 7(1563 - 462) + Oa"^ [56 - 22q(7 - 5a) + 77/3(1 - a + /3)] | . (A.94) 

From ([08]) it follows that 

-5r i-35r ^ - 7 

VoA = -a^21{axsf + 63/3(ax,) + 7a^ + 63/3^ - 9j , Vu = —a^^{axs) + 3/3| , V2A = -a^ , 

^05 = ^a{3(ax,)3 + 12/3(ax,)2 + (a^ + 18/3^ - 2) (ax,) + 4^{a^ + 3/3^ - l) | , 

Fig = -a2|l5(ax,)2 + 60;5(axj + 9a2+ 90/32 -loj ^ y^^ ^ ^q,3|(q,^j _^ 4^| ^ 
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t/o6 = —ah65{axs)'^ + 825/3(ax^)^ + 5(110^ + 330/3^ - 30) (axs)'^ 

+ 25/3(11q^ + 66/3^ - 18) (ax^) + 25 - 50a^ + 33a^ + 50(110^ -9)^"^ + 825/3^} , 
Fie = ^a'^hbiaxsf + 275/3(0x^)2 + (33a2 + 550/3^ - 50) (ax^) + 5j3{33a^ + 110/3^ - 30) | , 
V2e = —a^l77{axsf + 385/3(aXs) + 2(33a2 + 385/3^ - 35) | . (A.95) 

We refrain from explicitly spelling out the form of b^{W, W) for higher values of /i, as these functions 
follow directly from their definition ()A.48p . 

A. 5 Subtracted kernel functions 

Finally, we summarize the changes that are necessary if the subtracted versions of the HDRs are used 
for the s-channel projection. 

The modified pole terms are given by 



N, 



I in-sub 



/ in-sub. 



in-sub / 



'l+\" {W) = Nl^W) + AiV/+|'™"(Ty) = -Nli^^^_\'-'''\-W) 



2-sub,-,^,s 5lQ r /^2 



/\N+X {W) 



IGttW 



{E + 171) 



^ m ^ ' m 



-(E-™)f(HV(iv + „4)[ + ^ 



(E + m)—(Ad+ +iW- m)^) 
3 \ m J 



AiV^-f"'"^(VF) 



^m 



IGnW 



(E + m) 



{s-So-q^)^+2iW 



m] 



^ + 6o-o - 2q%, 



2) \ m 



{E - m)^ (^-4{W + m)b^^ ] \ + -^{E + m)^ ( ^^ + 4(PF - m)6, 



-sub SlQ 



8ttW 



{E + m){W -m) 



■'oi 

,2 



IGttW 



3 \ m 



01 



2m" 



+ 6, 



00 ' 



(A.96) 



where for convenience we have defined non-vanishing corrections also for the unsubtracted case ac- 
cording to (cf. ()A.5p l 



AN, 



J I O-sub 



l+\ 



(W) 



^ g^ {E + m){W -m) 5iq 
~47r 2W m?-a 



-^H+D- 



I O-sub , 



-W) . 



(A.97) 



The additional contributions to the s-channel kernels that fulfill the MacDowell symmetry rela- 
tion ()3.3p in both {W, I) and (VF', V) can be written for all {I > 0, /' > 0) in the symmetric form 

AKl,{w, W) = Kk\,{w, W) - Kk[,_^{w, -W) + akUw, W) - Kk[i,_,{w, -W) 
+ ^\^Kk\A-w,w') - AKl,_,i-w, -W) 
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AkU^-''''°{w,w') 



W 



k\W, W')h,{s') + ^-\W, W) y ''I )Pl,^, [ KJ,,,, ]5io 



ys' - So 

W r 5{W, W'){s' + s- 2so) + €^p{W, W'){s' - s) 

1^1 {S' - SqY 



(0,0) 



>i\W,W') 



s' — a 



}^'+i(K 



(0,0) 



^m 



^^-^^i^W,W')ho{s')n^,{[z', 






(0,0) 
''J (0,0) 



^10 



x'{W,W')ho{s')[dtz',],^,.P[!J[zi] 



I (0,0) 



;- 1-sub „ 
0/0 — > ^ , 



(A.98) 



where we have used that e = ite . Note that for /' = the term proportional to (s' — a) ^ 
cancels against the corresponding term in K^qIW, W) of ()A.65p as for the nucleon pole terms (cf. the 
relation (jXBTl) ). 

The additional contributions to Gij and Hu may be written as 



AGuiW, t') = AGij{W, t') - AG,+i, j(-VF, t') V (/ > 0, J > 0) , 
AHu{W,t') = AHuiW,t') - AHi+^j{-W,t') V (/ > 0, J > 1) , 



(A.99) 



where for even J 



J J\J 



AG„r-->(»'.0 = 4±;i(2J+ l)(^^ [P,(.;) 



2W 



t'p't' 



I (0,0) 



^IQ 



V(^ [^i(4)](o,o) + [5*^^(-n](o,o)) ('^'0 - f 



^ ^14^^^-^ + ^^^^ [^^(^*)J(o,o)^'o 



A^zj|'""^(H^,t') 



E + m 2J+1 (PtgD W^-"i 



/:V,/2 



2M^ VJ(J + 1) t'pi- 



291^ 



+ m[4P;(4)] U,o-2<? 



/2 (■S - ■so) 
i 

W — m 



P'AO 



(0,0) 



4?^ 



,/2 



P'jiA) 



(0,0) 



+ m 



7K^j(^*)].oo) + [^*(^*^^(^t))] 



and for odd J 

1 2-sub 



AG, 



ij 



AH 



u 



{W,t') 



f/ L-i.- J v-i,/J(0,0) 

.sub E + m 2 J + 1 (ff^gp-^ 
~^ 2W ^J{J + l) t'pf 

E + m,„^ , ^,{p[qiy-^l \Pj{ziy 



(0,0) 



^10 



Sii 



m 



WtP'A4) 



I (0,0) 



Sio 



(A.lOO) 



21^ 



-(2J + 1) 



/„/2 



t'p't 



(0,0) 



[s- So -q^)6iQ + q 



2^11 I 1-sub 







'^'\J~l 



m 



E + m 2J+1 {p[q[ 
' 2W yjiJ+1) t%' 1 V"* ^" "^^ ' T^-^ ^'^^y L^^v^*;j(o,o) 



pfiW-m) + -{s-so)][P'j{z[)\,^^.5io 
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-2q^p'f{W -m) 



+ t[^^(-*)], 



t' 



7[Pji4)]m) + [9^Pj(<)] 



I (0,0) 

Sio - f ) 



(0,0) 



4 L- j\-"J(o,o) 
i.sub E + m 2J + 1 {p'tQiy-^ 



2W //(J+T) t' 

Note that again only even powers of momenta and z'^ occur. 



{W-m)[P'j{4)].^^^5io 



(A.lOl) 



B Partial-wave projection for the t-channel amplitudes 

In the foUowing, we will discuss the different contributions to the t-channel part ()3.7p of the RS system. 

B.l Nucleon exchange 

In order to carry out the projection integrals (j3.6|) we rewrite s and u as functions of t and Zt via 



s{t, Zt) = -{T,-t + 'iptqtzt) , u{t, zt) = -{T,-t- ^ptqtZt) , 



which allows us to cast the nucleon pole terms of the HDRs (I2.44[) into the form 



1 ± ' 



1±1 



m"^ — s m? — u 2{m? — a) 



[t;zt] 



=F 



] '^Ptqt[y-zt {-y)-ztj 2(m2-a) 



libl 



(B.l) 



(B.2) 



where the upper/lower sign corresponds to even/odd J (i.e. to / = +/— ) and we have defined, in 
analogy to ()A.6p . 

y{t) = — = = zt{s = m^, t) = xtit, s' = m^) (B.3) 

wtqt ptqt 

{xt{t,s') will be defined in (JB.ISP ). By noting that the orthonormality of the Legendre polynomials 
yields 



libl 



1 

J dz Pj{z)Pi=2m{z) 



JJI 



1^1 



2/ + 1 



1 

J dz Pj{z)Pi=2n+i{z) VJ,/(m,nG No), (B.4) 





the nucleon pole terms of the PWDRs (J3.7I) can be written as (in analogy to (IA.3P ) 



^^(,)_5^™/yQ.(y) 



4vr i {ptqt 



sJ 



Sjo 



1 6ji 



3 m^ — a 



J.. _ g^ ^JJ{J + l) r Qj-M - Qj+i{y) 5n 



A^l(t) 



4tt 2J+1 



{ptqt) 



m^ — a 



V J>0, 

V J> 1 , 



(B.5) 



which for later convenience may be expressed as (in analogy to ()A.4p i 



Niit) = Niit) ' "^ ^'' 



Nl{t) = Nl{t) 



Att 3 m? — a 
g^V2 5ji 
Air 3 m^ — a 



V J>0, 



^_^(t) = ^ ^^2(Z±5^Zz4).Z^Z±iM , VJ>1. (B.6) 



47r 2 J + 1 



{PtqtY 
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Note that for t G (t7r,iAr) due to pt G iM. also y £ iM and hence we need the analytic continuations of 
Qi{z) as discussed in Appendix lA.il However, the pole-term projections (IB.5|) . (IB.6P are real for all t 
above the logarithmic branch point singularity at f,r — {M'^/rnf' ~ 3.98M^ of the nucleon cut (which is 
the left-hand cut for yit)"^ < 1 along the real axis due to the Zf-projection of the nucleon pole terms), 
since y/{ptqt) and the squares p^ and y^ are always real, and thus we can rewrite the projections solely 
in terms of real quantities due to the defined parity ()A.9P of the Qj{y)- Finally, we comment on the 
asymptotic behavior for p^q^ — )• 0, particularly including the vicinity of the aforementioned logarithmic 
singularity. The ostensible poles in (jB.SP are canceled by the asymptotics of Qj{y) for y — >■ oo. In 
this limit, we may abort the series representation of Qi{z) valid for \z\ > 1 [96] 



Qi{z) 



2'(Z!) 



n2 



^-{i+D + 



(/ + !)(? + 2) _ 



(1+3) 



+ 



(/ + !)(/ + 2) (/ + 3)(/ + 4) 



(2/4-l)!i 2(2/ + 3) 2(2/ + 3) 

after the first term and obtain the leading contributions 



4(2/ + 5) 



Niit) 



9 



J! 



rm 



47r(2J + l)! 

9^ J! [J+T 

47r(2J + l)!!V J 



t-2M3 



(5 JO 



Sji 



m^ — a 



+ 0{ptqt) 



t-2M2 



Sji 



m^ — a 



2„2\ 



+ 0{ptqt 



z-('+5) + ... 



V J> 0, 

V J> 1 . 



(B.7) 



(B. 



In particular, it follows that the leading contribution to N^{t) vanishes, such that N^{t) even involves 
zeros for ptqt — t- 0. However, higher orders need to be taken into account in the approximations (IB.SP 
in order to obtain precise numerical results in particular for gt — )• 0, since the pole terms vary rapidly 
in the vicinity of t^^. Note that (IB.SP and (jB.Sp reduce to the results given in [29] and [56] if the terms 
containing the hyperbola parameter a (that only contribute for J = 1 anyway) are dropped. 



B.2 s- and it-channel exchange 

We may rewrite the t-channel partial-wave projection (j3.6p in matrix form as 



i\t) = j dztT\t,zt)A.\s,t) 



s=s(t,zt) 



(B.9) 



where the projection kernel is given by 



f-'it,zt) = Cj 



UJ VJ 

wj 



Ut) 



1 



4vr(ptgt) 



J-i 



(B.IO) 



Pj^iizt) - Pj+i{zt] 



uj{t, zt) = -Pjizt), vj{t, zt) = mztPj{zt), wj{t, zt) = \^ j , . — 

Qt IJ + I 

For the following, we need the matrix form of both s- and u-channel HDR terms (|2.44p according to 



A\s,t) 



s+u \ 

s=s{t,zt) vr 



//„/ J./^ 



ds h^[t, s ; Zi] Im A (s ,t 



t'=t'{s',z',) 



(B.ll) 



■5+ 



where the kernel matrix h{ is given in ()A.2ip . and [t, s'; zt] indicates that the whole integrand is to be 
understood as a function of these variables, which can be done by using (jB.ip and thereby 



{hi^h2^ 



[t,s';zt] 



1 1 

± 



1±1 



s' — s s' — u 2{s' — a) 



1 



1 



1 



1±1 



[t s'-zt] '^Ptlt [xt- Zt {-xt) - Zt) 2{s' - a) 

(B.12) 



73 



The upper/lower sign corresponds to even/odd J and we have defined in analogy to ()A.2ip 

,, t + 2s'-S 



Xt{t,S 



4^Ptqt 



zt{s',t) 



(B.13) 



According to ()A.72p . the relation between z'g and zt in (jB.lip is given by 



11,1 \ ^t " 
Zs{t,S ■,zt) = — z — 

7 



~u '\ Q (.s -a) , 

71*, s ) = ^ 2 2 = 7('S , t) , 

d[t, s ) = „ „ = d[s , t) . (B.14) 

IGpfqf 



Expanding the absorptive parts of (jB.lip into s-channel partial waves via ()A.16P and projecting onto 
i-channel partial waves by means of (jB.OP leads us to the PWDRs for the t-channel partial waves 



s+u 1 / , _ ,, , 

{■^(t) =- dW ^G-^\t,W')lmii\W') , 



(B.15) 



with the kernel matrix 



G-^'it^W) = 2W' / dzt T-'{t,zt)hi[t,W';z',]S\W',z',) . 



As a remnant of the MacDowell symmetry, (JA.IOP induces the symmetry property 

G'^\t,-W') = G-^^{t,W')ai , 
such that the parameterization with two kernel functions 

^ ^''""^ \Hji{t,W') Hji{t,-W') 



(B.16) 



(B.17) 



(B.18) 



is justified, where again according to /£ = we set Hqi = for the matrix notation. This reproduces 
the s- and u-channel part of (13. 7p 



Hit) 



Hit) 



-u 1 



hu 1 



1 r 

- / <^W'Y,[Gji{t,W')luifl4W') + Gji{t,-W')linfli^^^_{W')] VJ>0, 

If 

- j dW'Y,{Hjiit,W')lrafl^{W') + Hji{t,-W')lTafli^^^_{W')] VJ>1. (B.19) 



w+ 



If we introduce the abbreviations (cf. (jA.SSP ) 

i;[akn\diW')] = d{W')ak,n+i + di-W')akn , Vjit,W') 
we find for the kernel functions 



2W' 



iPtQt) 



J-l 



(B.20) 



Gji{t,W') = fjj\-^i> 
Qt 



A 



Ji 



W + m 



Hjiit,W') = ^j^^^^ 



E' + m 
Cji 



+ mip 



B 



Ji 



1 



E' + m 



1 



E' + m 



V J>0, 

V J> 1 , 



(B.21) 
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where the angular kernels are given by 



i 
Aji{t,s') = Jdzt Pj{zt)[hi T h2}p!{z',: 



lt,s';zt] 



1 
Bjiit,s') = fdzt Pjizt)zt[hi±h2]Pliz';, 



lt,s';zt] 



Cji{t,s') = J dzt [Pj-iizt) - Pj+i{zt 


Decomposing these kernels according to 
1 



hi±h 



4^/(4: 



[t,s';zt] 



1 



^J-l,l - ^J+l,l 



(B.22) 



Aji{t, s') = -^Pl{zs)Qj{xt) - Ajiit, s') , Bjiit, s') = -^Pl{z,)xtQj{xt) - Bji{t, s') , 
Ptqt ptqt 

Cjiit, s') = —P[{zs) [Qj-i(xi) - Qj+iixt)] - Cjiit, s') , (B.23) 

Ptqt L J 



with the real quantity 



Zs{t,s') 



X? — S 



1 + 



2q' 



12 



Zs{s',t) 



(B.24) 



and polynomial parts defined by 

1 

2 



Aji{t, s) = - I dzt Pj{zt) I ^^ T — : + TTTi^ — —Pi (^J r' ' 



_ -ptqt xt - zt 



2{s' - a) 
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B,,(U) = i /d., PX.,) j^M(!;)^iiiM) ^ J^.,PH, 



-1 



Ptqt Xt - Zt 



2{s' - a) 



l\^s) ( > 



Cji{t,s') = ]^ I dzt 



Pj-iizt) - Pj+iizt] 



1 P/(z,)-P/(4) , iTl _^,^^,^ 



Ptqt Xt - Zt 



+ 



2(s' - a 



the kernels Gji and Hji may be written in a recursive fashion 

Gji{t,W') = Gji{t,W')-Gjj_i{t,-W'), Gj,_i = 0, VJ>0, 
Hjiit, W') = Hji{t, W') - Hj,i-i{t, -W') , Hj^^i = , V J > 1 , 



(B.25) 
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E' + my Ptqt 
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qt 
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J,i+i f ' 



(B.26) 



keeping Gj^i just for convenience. Note that since xt/{ptqt) and the squares pf and x^ are always real, 
Aji is real/imaginary for J even/odd and the other way around for Bji and Gji. Therefore, we can 
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conclude that the functions Gji, Hji and hence the kernels Gji, Hji are real for t > t^^ — (M^/m)'^, 
cf. the discussion following ()B.6p . The kernels for all combinations (J > 0,/ < 2) explicitly read 
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s' -a[ 231 
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(B.27) 



The explicit formulae for the polynomial parts Aji, Bji, and Gj^i for (J < 2,/ < 6) needed for 
calculating these kernels and furthermore the additional kernels Gji and Hji for ( J < 2, 3 < / < 5) 



76 



via ()B.26P are given in Appendix IB.4J^^I As a check of our calculation we can reproduce the nucleon 
pole terms (IRK]) by (cf. fOTl) ) 

Nl{t) = -fGjo{t, -W' = m) V J > , N:[{t) = -fHjoit, -W = m) V J > 1 . (B.28) 

The asymptotic behavior of the general kernel functions (IB.2ip can be deduced to be 

for ptqt ^ Gjiit, W) ~ Hjiit, W) '^ 1 , 

for q'^0 Gjiit, W) ~ Hjiit, W) ~ q'"^' , Gji{t, -W) ~ Hji{t, -W) ~ q'-^''^ , 

for g' ^ oo Gjiit, W) ~ Hji{t, W) ~ q'-^^ . (B.29) 

In particular, these kernels are finite for ptqt — ?• and their precise form in this limit may be worked 
out in close analogy to the discussion of the pole terms in Appendix IB. II based on (IB.26p . Note that 
both ()B.28P and ()B.29P obey the MacDowell symmetry relation (|3.3p . as they should. 

B.3 t-channel exchange 

We need the t-channel HDR terms (|A.70P in the form 



A^(s,t) 



= - / dt'h[[t,f';zi]ImA^(s',t'^ 

s=s{t,zt) vr 



^'=s'it',z',) 



(B.30) 



where the kernel matrix h^ is given in ()A.7ip . and the integrand can be written as a function of the 
variables [t,t']Zt] by noting that 



1 1 



2g2 xt - Zs 



[t,t']Zt] 



t'-t' 



V 



_ Ptqt zt 

[t,t';zt] Pt^t ^t 



(B.31) 



and that z'^ and zt are related by (cf. (|A.22p ) 



PtQt 



z[{t,t';zt) = Jaz^ + J3 , a{t,t') - ^^y. , 
^ PtQt 



/3(*' t') = TTTTTli'^ + t' - 2S + 4a) . 



(B.32) 



Wt'q't' 

If we expand the absorptive part of (|B.30p into t-channel partial waves by using ()A.75p and project 
onto t-channel partial waves again via (jB.Op . we obtain the following PWDRs for the t-channel partial 
waves 



f^(t)|* = - fdt' J^K^-^'(t,t')Imf-^'(t') 
y J' 



(B.33) 



where the summation runs over even/odd values of J' for even/odd values of J, accordingly, and the 
kernel matrix is defined by 



K^-^'(t,t') = / dzt f'{t,Zt)hi[t,t';zt]T-'\t',z',) . 



(B.34) 







^^Note that for \a\ — > oo, of all polynomial parts only Boi does not vanish completely and hence /+ receives polynomial 
contributions from the kernels Goi- These remaining contributions, however, are just those that cancel with the leading 
terms of the S-wave pole terms (|B.5|) . cf. the discussion following (|B.8|) as well as the explicit kernels (|B.27|) . 
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Calculating this kernel matrix shows that it can be written with three kernel functions as 

jyjJ' u ,/^ _ fKJj' (*, t') Kjj, (t, t')\ _ Cjr_ (ujj, (t, t') vjj, {t, t') \ 
^'^ V Kjj,it,t')) t'-t{ wjj,{t,t')) ' 

where we have defined different angular kernels for even J and J' 



(B.35) 



i 

ujj' = -^ / dzt Pj{zt)Pj'i4) , VJJ, 
QtPt J 



m 



^J'{J' + 1) QtM 



tPt'it J ^ ■' 



P'j'{4) 



WJJ' 



2J + 



r\/^^l/^-K''^''-^-<=''} 



Zf 



P'j'iO 



(B.36) 



and for odd J and J' 

1 



Pt 

UjJ' = -p2 
Pt 



/ dzt Pj{zt)zt 



Pj'i4) 



m 



VjJ' 



VJ'iJ' + 1) 



Pt_ 
P? 



i 

j dztPj{zt)ztP'j,{z[) 



WJJ' 



2J + 



rV:^!^^ /''-'' {^^-''^'>-^^-'^''}^^''^'>' 



In this way, we recover the form of the i-channel part given in (j3.7p 

oo 

fiit)^ = ^ fdt'Y,{Kh'it,t')lmfi'{t')+KJj,{t,t')lmfl'{t')] 



t-ji 



m) 



t 1 



TT 



dt'Y,k]j,{t,t')i^f':{t') 



(B.37) 



V J>0, 



V J > 1 , (B.38) 



and according to /£ = we set K^j, = = K'j^ . 

From the projection integrals ()B.36p and (IB.37P together with the definitions (IB.32P and 



Pt _t'-t 

P? ¥? ' 



Qt^t -Qt h 



t' -t 



¥/ 



^2 ^2 .'2 ./2 _ _ ^ 4^2^2 _ _ (i + ^z _ 2S + 4a) 



(B.39) 



one can see that the off-diagonal term vjj' is proportional to t' — t, as it should be. Note also that only 
even powers of z^ and zt occur in the projection integrals. Therefore, the kernel functions Kjj,, Kjj,, 
and Kjj, are always real, since the prefactors contain only even powers of momenta. The integrals 
can be performed with the help of |96] 



i-i 

2 



pli^)=E 



ev 2A 



Pli^)=E 



„od 2A+1 



A=0 



A=0 



for even and odd values of /, respectively, where 
(-1)^+^(2A + /-1)! 



aTi 



2'-i(f-A)!(A + |-l)!(2A)! 



alf 



-1)^+^ 



(2A + /)! 



2«-i(^_A)!(A + ^)!(2A + l)! 



(B.40) 



(B.41) 
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which also follow from reordering the expansion 



2'^A!(/-A)!(/-2A)! 



2 for even Z , 
^ for odd Z . 



In this way, the required non- vanishing integrals may be written for even J and J' as 

1 

A ^ TD^t ^ \TD ( J\ _ Y^ „CV \ ^ I ^ \ Xil^ /3A'-M- 



1 4 y ^ 

A'=4 /.=! 



a^/3^'-^a}; , 



|dz, PM){qlz', - q?z?}^ = Y: 2A' 



f-AV 



A'=max{^,l} 



X S9t 



A'-l 

E 

^=max{4T — 1,0} 



(B.42) 



^'^ ') c.^P^'-'-'az+i - Q? E (f ) «"/5^ -'^«}; f , 



rM/5A'-/i~cy 



± 

J dzt [Pj-iizt] - Pj+i{zt)} 






2 



E 2A'a' 

A'=max{|-,1} 



cv 
A' J' 



A'-l 

E 



/i=inax{ -^ — 1 ,0} 
A'-l 



A'-l 



r,^3A'-l-M-od 



" P "j-1,^+1 






A' = ^+l 



(B.43) 



and for odd J and J' as 

1 



j'-i 

2 



d., p,(z,)..^ = E «v.' E (tV'^'"'«°^:^+i 



dzi Pj{zt)ztPj,{zt) = 
A 



^^ A' 

E (2A' + l)a°A E 0') 

/_ ,7-l .,_ ,/-! \^/ 



~^3A'-^~od 
" P ^J,u+1 



1 

jdzt[pj-M)-Pj+M)}p'Az't) 



^ A' 

E (2A' + l)a°,?,, E (^')«'^^''~^«J-i„ 

, J-1 J-1 Vr^/ 



M=^^ 



^ A' 

E (2A' + l)a-^ E 0') 



a'^/^^'-^a^-^. 



J+1,M ' 



(B.44) 
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with the definitions (for even and odd values of J, respectively jfj 



a 



ev 



y- "jj _nJ (^ + 4)K2m)! _ (2/i)! ( J_ 

^2(// + A) + l (^_^)!(2;, + J + i)! (2/i- J)!!(2^ + l + J)!! V 2 



A=0 

j-i 



a 



^+^)!(2^-l)! _ (2^-1)! / ^ '^-l 



^^^2(/. + A) + l (;,_m)!(2;, + j)! (2m-1-J)!!(2^ + J)!! V 2 

(B.45) 
We can conclude that the following kernels vanish: 

K-^J'{t,t') = \/J'<J, (B.46) 

and by using the identities 

(2J + l)a7 ,a7j = 1 , JaJ ,a°^, j = 1 , for even J , 

2 ' ' 2 2 ' ' 2 

(2 J + l)a°Ai ,<2+i = 1 , Ja°/i .a7_^ j^ = 1 , for odd J , (B.47) 

it follows that the non-vanishing kernels for J' = J take the form 

■^■^^ ^ pft'-t t'-t t'-tN t't'-t t't'-tN ~ 

J m J m 



Kjj{t,t') = \- -—-r = J- — VJ>1, 

^jj(i,0 = ^ VJ>1, (B.48) 

from which one can immediately read off the relation (valid for all J) 



k]j{t,t') = mJj—^[K^jj{t,t')-K]j{t,t')] . (B.49) 



J 

f+ 
This together with 

i^02(i,O = ^%{i + i'-2S + 6a}, i^2^(t,t') = -^%{4g2-3(t + t'-2S + 4a)} 



^2 



^i3(*'O = ^^{i + i'-2S + 10a}, i^2^(t,t') = ^^{8p,V + (i/_,)(i + i/_2S + 5a)}, 

Kf3(t,t') = ^{t + t'-2S + 5a} , (B.50) 

completes the calculation of the t-channel kernels with (J < 3, J' < 3). Finally, from ()B.43p and 
()B.44p we may infer the asymptotic behavior of the non-vanishing kernels 

for pt ^ k]j,{t,t')^pi , Kjj,it,t') ~ kh'it,t') ~ 1 , 

for qt ^ k]j, (t, t') ~ Kjj, {t, t') ~ KJj, (i, t') ~ 1 , 

fort ^00 K]j,it,t')^KJj,{t,t')^t''-' , k^jj,{t,t')^f''-'-\ 

forpJ^O K]j,{t,t')r.KJj,{t,t')r.p'f\ k'jj,{t,t') ^ I , 

for q[ ^ k]j,{t,t') ~ i^3j,(i,i') ~ Kjj,{t,t') ~ 1 , 

for t'^00 k]j,{t,t') ~ t'^'-^-2 , K]j,{t,t') ~ ii-5^,(t,t') ~ t'J'-J~^ . (B.51) 

^^These identities are similar to the Saalschiitz formula (|A.42|) employed in [46]. Note that (—1)!! ~ 0!! = 1. 
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Note that the kernel 1^02 (t,t') exceptionally has better convergence properties 



Ki^it, t') ~ pj for pt ^ , kl^{t, t') ~ 1 for t' ^ oo . 



(B.52) 



B.4 Higher kernel functions 



The explicit form of the polynomial parts Aji, Bji, and Cj^i from (IB.25|) for (J < 2,/ < 6) that are 
needed in order to calculate the kernels Gji and Hji for ( J < 2, / < 5) via ()B.26P explicitly read 
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Aii = 0, 
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B.5 Subtracted kernel functions 

Here, we give the modifications of the nucleon-pole-term projections and the kernel functions for the 
i-channel projection, that are required by the subtractions performed in Sect. 14.31 

To start with, the n-times subtracted nucleon-pole-term projections may be written as 

iV^r"""'(f) = iV^(t) + ANir'-^it) , 
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127r V 2m2 
where in analogy to the s-channel projection we have defined unsubtracted corrections (cf. (jB.6p ) 

ANif-^\t) = -^"^J^ , AiV^r^^^t) = -^^4^ , (B.62) 

' 47r 3 m^ —a 47r 3 m^ — a 

which are constant and non-zero only for J = 1, in order to split off all terms that are either constant 
or contain subthreshold parameters. Note that for both one and two subtractions the full nucleon-pole- 
term projections fulfill the threshold relations (j3.63p for pt — >• 0, but no longer for qt — )• 0. However, 
the subtraction-independent parts of the pole terms NJ. still fulfill the relations (I3.63P for ptqt — ?■ 0. 
The necessary update of the s-channel kernels Gji{t, W') and Hji{t, W') may be achieved by adding 
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respectively, to Aji, Bji, and Cj; at the pertinent places in ()B.26p . leading to corresponding AGji 
and AHji. Note that also in both the once- and twice-subtracted case ACji = AAj^i^i — AAj^i^i is 
still valid (for J > 1, here actually ACji = AAj_ij). 

The additional contributions to the t-channel kernels K {t,t') amount, for even J and J', to 
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while for odd J and J' one finds 



^K'jj' 



Pj'iz't) 



(0,0) 
2 



^Jl 1-sub 



AK3,,|--(t,0 = -^==iPdy'-'y{{l - I + ^) [^M4)](o,o) + t [d,P:Az',)] (0,0) }% 



v^ 



i_sub 2J' + 1 .,/„/^J'-l^ [p/ ('^/^i ^1 

1(0,0) 3 



J j'{J' + 1) ^^*^*^ 7 [^J'(^t)] (0,0) ^- ' 



AKU'-'itA) = - 'rv'^l.^ pdy'^'fii' + ^) [^M4)] (0,0) + ^ [d^P'A4)] (0,0) }% 



a/J'(J' + 1) 

^SUb 2J'+1 , , ,.j'_lV2 r -, 6ji 
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Furthermore, AK = for J > 2 or J' < J, the latter being in agreement with (IB.46p . In all cases 
only even powers of z'^ and the primed momenta occur, and hence the additional kernel terms are 
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always real. Here, we refrain from explicitly expanding the Legendre polynomials using ()B.40p as in 
Appendix IB. 3| but only give one example to demonstrate this point (for even J and J') 



[dMP'Az't))]^,,^) = 2 [5.41 (0,0) E«A>A^ [• 



2 r~/2i^-i 



A=l 



1(0,0) • 



(B.66) 



For later convenience, we explicitly state all those subtracted kernels with < J' < 3 that differ from 
their unsubtracted form 



Kl,t'''\t,t') = ^i^o^o(M') '^ l;KUt.t') , Kl,\'-\t,t') = Ul,{t,t') '^ KlAt,t') , 



f^l |2-sub 



K 



02 I 
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1 + 
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t' At' t' 



So 



^N't-K 
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Kl,{t,t' 



K 
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11 



'{t,t') = ^Ki^{t,t') ^^ fKt,{t,t' 



r^2 |2-sub 
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02 
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?>2 |2-sub. / 
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t \ tArtTT t 

1 + - 1^-^ - -So 
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still obeying the threshold-behavior relation (j3.62p . Note that at the level of two subtractions all these 
kernels are independent of a (which is, however, not true for only one subtraction and J > 3 or without 
subtracting), and that the exceptionally safe behavior of KQ^it^t') at t^ is preserved (cf. ()B.52p ): 



r>2 in-sub 



K, 



02 1 



{t^tN,t') = 0{p^t) Vn>0. 



(B.68) 



C Ranges of convergence 

In this Appendix, we will analyze both the convergence of the partial-wave expansion of the imaginary 
parts inside the integrals and the convergence of the partial-wave projection of the full HDR equations. 
For the rest of this section we may work as if no subtractions were necessary. 

C.l Boundaries of the double spectral regions 



The following analysis is performed in the spirit of |19 |[29l[3T|rn The basic assumption is that the 
T-matrix element (and hence the scattering amplitudes {^(s, i), i?(s, t)} oc T(s, i)/(167r)) fulfills Man- 
delstam analyticity jQ^, i.e. that it can be written in terms of double spectral density functions psu-, 



Note that the authors of [19] corrected their discussion of the boundaries of the double spectral regions in [98] . 
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(Ill) 



(IV) 



Figure 14: Box graphs constraining the boundaries of the double spectral regions. Solid lines denote 
nucleons and dashed lines denote pions. 



Ptu, and Pst according t' 



Tis,t) 



1 



TT' 



ds'du' 



Psu{s',u') 

(s' — s){u' — u) 



+ 



vr^ 



dt'du' 



Ptu{t',u') 

{f -t){u' -u) 



+ 



1 



TT' 



ds'dt'- 



Pst{s',t') 



{s'-s){t'-t) 
(C.l) 
The integration ranges are determined by those regions in the Mandelstam plane where the corre- 
sponding double spectral densities have support. The boundaries of these so-called double spectral 
regions will be the central objects of the following discussion. 

The three double spectral densities can be derived by studying the consequences of unitarity in 
the 2-interniediate-particle approximation. We consider the corresponding lowest-lying intermediate 
states as depicted in Fig. [T3] (as unitarity diagrams, i.e. with on-shell intermediate particles), where 
the inelastic (referring to the intermediate state of the s-channel process) diagram (I) and the elastic 
diagram (II) yield the boundary of the support of pst (from which, due to s -f-T- n crossing symmetry, 
the result for p^t directly follows), while (III) and (IV) are relevant for calculating the boundary of 
the support of psu- This leads to boundary functions (cf. f29]) 



4M2)A(s,m2,4M2) -16M^(s + 3S_) 



6i(s,t) = 
for the boundary of pst and thus bi{u, t) and bu{u, t) for the boundary of put-, as well as 



(C.2) 






16Mi 



16Mi 



m su 



m su 



Y?im' 



S2_fm2 



t{s,u)) 
t{s,u)) 



(C.3) 



for the boundary of psu, where we only need to consider 6iii(s, u) = 6iv(^, s) due to s o n symmetry. 
The whole support of all three double spectral densities is then given by the union of the regions 
allowed by the non-trivial constraints that the corresponding boundary functions be non-negative. 
Furthermore, trivial constraints arise from the lower kinematical bounds of the corresponding physical 
regions that are given by the asymptotes of the boundary functions in question, e.g. for the inelastic 
diagram (I) we find the asymptotes s = {m + 2M,r)^ and t = (2M^)^ = 1.,^, and for the elastic diagram 
(II) we obtain s = (m + MT^)'^ = s_|_ and t = (AMj^)'^. Therefore, by defining the following abbreviations 
for the solutions of the implicit equations 



bi/iiis,t) = 



t = Ti/uis) 



Si/iiit) 



(C.4) 



^^Mandelstam analyticity can at any rate be justified in the framework of perturbation theory [99H101| . While for mv 
scattering the vahdity of the Mandelstam representation can even be shown rigorously in a finite region [10211103] , for 
nN scattering (involving unequal masses and spin) at least the uniqueness of amplitudes satisfying this representation 
is ensured by the MacDowell symmetry [10411105] . 
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Figure 15: Double spectral regions for vrA^ scattering (shaded) and boundaries of pst and put reflected 
in the corresponding scattering angle (dot-dashed). 



the boundary of the support of e.g. pst is described by 



Tst{s) = min{ri(s),rn(s)} 



Tiiis) 



for s+ < s < {m + 2M^) 



min{Ti(s), Tii(s)} for (m + 2M^f < s 



(C.5) 



with the functions 



2\2 



Tiis) 



Ttt(s) 



4M2(s-m2-2M2) 
A(s,m2,4M2) 

16M2(s-S_1 



> 4M2 V s > (m + 2M^f , 



A., 



> 16M^ V s > s+ , 



(C.6) 



again limited by the physical constraints, such that by definition Tst{s) > 4M^ for s > s+. The 
boundaries of all three double spectral regions are shown in Fig. [151 The asymptotes of pst are s = s+ 
and t = tn and hence those of put are it = s+ and t = t-,^, while the symmetric asymptotes of psu are 

s = S-^ and u = s+. 



C.2 Lehmann ellipse constraints 

The boundaries of the double spectral regions limit the range of validity of the HDRs in two ways: 

1. The partial- wave expansions of the imaginary parts inside the HDR integrals (internal/primed 
kinematics) in the unphysical regions for both s- and t-channel partial waves converge only for 
CMS scattering angle cosines z' within the corresponding large Lehmann ellipses |106j. These 



ellipses are the largest ellipses in the complex z' plane, centered at the origin with foci at z' = ±1, 
that do not reach into any double spectral region. 

2. For a given value of the parameter a, the hyperbolae {s — a){u — a) = b with asymptotes s = a 
and u = a must not enter any double spectral region for all values of the parameter b that 
are necessary for the partial- wave projections of the full HDR equations (external/unprimed 
kinematics) in given kinematical ranges. Trivial geometrical constraints on a arise already from 
the asymptotes of the double spectral regions. 

In this section we will show how the (large) Lehmann ellipse constraint can be translated for a given 
a into a constraint on b, each for both the expansions in s- and t-channel partial waves. For any 
allowed fixed a, the allowed values of b are those fulfilling both of the above requirements, and the 
(limited) freedom in the choice of a in the construction of the HDRs can be used in order to optimize the 
convergence properties of the PWHDRs. In the two subsequent sections we will investigate numerically 
how these limits on b (for given a) yield the ranges of convergence of the full RS system via the 
restrictions that are necessary for both the projections onto s- and t-channel partial waves to converge. 
For the partial- wave expansion of the s-channel contributions, the Lehmann ellipse constraint 
states that the expansion converges for angles Zg{s',t') = 1 + 2s't' /Xgi (cf. (j2.5ip ) inside the ellipse 

(Re 4)^ jlmz'.r _ 

—42— + — g^ - 1 , i^-V 

with foci at z'^ = ±1 (corresponding to the physical constraint —1 < z'g < 1), i.e. semimajor and 
semiminor axis As and Bg are related by 

A's-B's=l. (C.8) 

Since for given t' the angle z'g is always real in the integration range s' > s+, the maximal value of z'^ 
for given s' not entering the support of pst follows from the corresponding maximally allowed value of 
t' (according to (jC.SP for the internal (primed) variables) and thus reads 

z^--{s') = l+'^Tst{s') = As ys'>s+. (C.9) 

As' 

From the geometrical condition —Ag<Zg<As then follows 

- -^r"" <z's< zT"" V s' > s+ , (CIO) 

and the lower bound due to this reflection in z'^ is actually stronger than the restrictions imposed by 
Psu as shown by the dot-dashed line in Fig. [151 where the z^-reflected boundary of the support of pst 
for 1/ > is given by 

u(^s,t(^s,-Zs{s,t = Tst{s))^) =^+Tst{s), (C.ll) 

with the asymptote u = 4M^ for s — t- cxd due to Tst{s). Furthermore, due to s -^ u symmetry p^t 
yields exactly the same constraints aspst (including the z^-reflected boundary for v < 0), and hence 
we only need to consider the latterO The possible values of t' for given s' are then restricted by 

(cf. m) 

-^-Tst{s')<t' <Tst(,s') ys'>s+. (C.12) 



^Note that both the s- and ^-channel physical regions fit well in between pst, Put, and their reflected boundaries. 
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Via the linear relation (j2.43p for the internal kinematics this range for t' can be translated into a range 
of allowed values of b{s' , t'; a) for given a according to (cf. (IC.lip ) 

b~{s',a) <b<bt{s',a) \/ s' > s+ > a , 
b~ {s, a) = {s' - a) (S - s' - Tst{s') - a) , 

bt{s',a) = is' - a)(s - s' + ^ + Tst{s') - a) = {s' - a)|^ + T,t(s') - a\ , (C.13) 

where we have used that from the asymptotes s = s+ and n = s+ of the double spectral regions it 
is geometrically clear from Fig. [15] that the allowed values of the hyperbola's asymptotic parameter 
a are trivially limited to a < s+ (independent of b), and hence we have s' > a for all s' > s+. By 
invoking the asymptotes 3 = 1-^ and u = t-^ of the z'-reflected boundaries of pst and put (cf. (jC.lip ) 
we can deduce that the allowed range of a is actually geometrically limited by a < tj^, which is the 
reason why the "fixed-t limit" \a\ — )• oo actually reduces to a — )• — oo. Now, we may define the highest 
lower and the lowest upper bound 

6~(a) = max6~(s',a) , ^^(a) = min bf{s',a) , (C14) 

s'>s+ s'>s+ 

as the maximum/minimum value of bs (s', a) within the integration range s' > s+, which then finally 
determines the allowed values of b for given a by 

K{a)<b<b+{a) Vs'>s+>a, (C.15) 

for the s-channel parts of the HDRs. 

The Lehmann ellipse constraint for the partial- wave expansion of the t-channel contributions limits 
the convergence of the expansion to angles z[[s\t') = mv' /[p[q'^ (cf. (j2.5ip ) inside an ellipse similar 
to (IC.7P centered at the origin with foci sX z[ = ±.1 

e^ + fl!^.l. A-m^.. (C.16) 

The argument for the t-channel contributions is more intricate, since inside the integration range 
t' > t^ the angle z[ becomes purely imaginary for t^ < t' < tjy, and hence no relations similar 
to (JC.IOP are possible. However, as the relation between Zj and b is non-linear anyway (cf. ()2.52p ) 

_^,^_ {t'-^ + 2af-Ab{s',t';a) 

h — 1« ,2 12 ' y^-^') 

16pfq/ 

where all squares are real but not necessarily positive, we are interested in the resulting Lehmann 
ellipse constraint for z^^. By squaring equation ()C.16p for general complex z[ we arrive at 

(R.{.n-l)^ (Imtf})' _ 

which corresponds to an ellipse in the complex z[^ plane shifted to the right by (Af — Bf)/2 = 1/2. 
Hence, it is centered at (1/2,0) with the semimajor and semiminor axes given by 

At = ^±^ = Al-]^, Bt = AtBt = A^A^-l, (C.19) 
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such that the foci are at 1/2 =p \J -^ — Bf = 1/2 =p 1/2 (corresponding to the physical constraint < 

z't^ < 1). Since for t' > t-„- we have z'^^ = Re {z[^}, the geometrical condition 1/2 — At < z'^^ < 1/2 + At 
leads to the analog of (JCIOP 

1 - A? = -B^ < z't^ < A't , (C.20) 

where it is important to note that on the right-hand side the relation between z't and At is not fixed 
due to the squares, while the reflection bound on the left-hand side is again more restrictive than the 
corresponding bound due to psu, and hence we only have to look at the boundaries of the support 
of pst- For the following it turns out to be advantageous to rewrite the boundary functions biji{s,t) 
of ()C.2p in terms of (z^, t), since the quantity v[zt,t) = ptqtZt/m is always real 

bi{u,t) = {t-4M^)ij{t-Amu + 6M^f-16m'^M^\ + 8M^h-4mu-^-6^-\ =0, 

buiiy,t) = {t-16M^)i-{t-4muf -4m'^M^\ + 32M^h-Amu-J:\ = 0. (C.21) 

Solving these implicit quadratic equations for ^{t) yields the physical solutions (cf. (29]) 



(t - 2Ml) it + 4M2) + 8M^VtJ(t - AM^)m'^ + M^ 



it - 8M„)^ + AM^^/tJit - l&M^)m? + IGM^ 
^,(*) = 4^(t"l0Ma >° V<>16M,^=4<,. (C.22) 

again limited by the physical constraints, where each sign of the root is fixed by zt{i', t) = mv/{pt(lt) oc 
+v and hence zf^^ = +mv'^^^ / [ptqt) in the physical t-channel region t > Am? = tjy. Defining the 
(positive) combined upper bound on u according to 

AT (^\ ■ r (^\ u\\ /^i(*) fort^<t<4t^, 

N,t{t) = mm{ui{t),uu{t)} = i . (C.23) 

\niin{ui{t),uu{t)\ for 4t^ < i , 

and resorting to the geometrical constraints of the original t-channel Lehmann ellipse ()C.16P for z't, 
the maximally allowed value of the real angle z't = Re z't for given t' > t]\f not entering the support of 
Pst is given by 

777 

^r"(i') = ^T^^*(*') = At yt'yt^, (C.24) 

PtQt 

and thus (|C.20p in this case leads to 

2 2 

1 - -^Nstit'f < z't^ < -^Nstit'f yt'>tM. (C.25) 

Ptflt PtQt 

In contrast, for t-^ < t' < t^ we have p't = ip'_ with real p'_. Accordingly, for the purely imaginary 
angle z't = i Im z't it follows from ()C.16P that 



mu 



, ... 1 1 1 1 7 'III ^ 'I i"* 

iim Zt[t 

- - p_qt • ■ - p^-q-^ 



p'-q't 



<^Nst{t') = Bt => B^ = -^Nst{t'f yt^<t'<tN, (C.26) 



which plugged into ()C.20p yields 



Nstit'f <z't^<l- -^Nstit'f V i^ < t' < tjv . (C.27) 



^ AT Jf'\-i < z'2 < 1 _ ^ AT ^^/^2 
„/2„/2^'^^nM SZt Si ,2 12 

Pt It Pt It 
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However, from both ()C.25P with pf > for all t' > tj^ and ()C.27p with pf < for all U < t' < t^ we 
arrive at the same constraints on u''^ for given t' > t-j^ (cf. |29j ) 

^ - N.tit'f < u'^ < Nstit'f yt'>U. (C.28) 

By virtue of the linear relation (cf. (|2.42p ) 

IQm'^u'^ = (t' _ S + 2af - 46 , (C.29) 

this range for v'^ can then be translated into a range for b{v'^,t';a) according to 

Kii'i a) <b< ht{t', a) M t' > U > a , 

hj{t\a) = h,t' - S + 2af - Am^Nst{t'f , 

6+(t', a) = i(t' _ s + 2af - 4pf g^^ ^ 4^2^^^^^/^2 ^ (^/ _ ^-^^ + ^2 ^ 5.2 ^ 4^2 ^^^(^/)2 ^ (q_3q) 

where we have included the geometrical constraint on a as discussed below equation ()C.13p . Defining 
again the highest lower and the lowest upper bound 

6^(a) = max6^(t',a) , bfi^) — min6^(t',a) (C.31) 

as the maximum/minimum value of b^ (s',a) within the integration range t' > t^^, we can finally 
give the range of allowed values of b for given a by 

Kia) <b<bt{a) V t' > t^ > a , (C.32) 

for the i-channel parts of the HDRs. 

C.3 s-channel partial-wave projection 

As mentioned before, it turns out that the constraints due to put and psu are equal to or weaker 
than the restrictions due to pst- Therefore, we only need to consider the corresponding constraints 
for the s-channel partial- wave projection of both the s-channel partial- wave expanded and the t- 
channel partial-wave expanded HDR parts. However, the strategy to find the optimal value of a and 
the corresponding range of convergence in s is the same in both cases: from the Lehmann ellipse 
constraint it follows that all allowed values of b must obe\o 



KAa)<b<bt/a), (C.33) 

for all s' > s_(_ and t' > t,^, i.e. within the corresponding integration ranges, respectively. The limits 
— 1 < -Zs ^ 1 of the scattering angle for the physical s-channel reaction translate into 

-4g2 = _^<t<0 ys>s+, (CM) 

s 

and hence for given a < s+ < s the bounds on b due to the s-channel partial- wave projection are given 

by (cf. dmsi)) 

6f ""(s, a)<b< 6r^(s, a) y s > s+ > a , (C.35) 

6f>^(s,a) = (s-a)(S-s-a) , 6r"(5, a) = (« " a)(s - s + — - a) =(s-a)|— -a| . 



"Note that the lower bounds coincide: bs (a) = 6^ (a) for all a < s+. 
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The maximally allowed value of s for given a, s^f^{a), is then the largest value of s such that for 
given a both 6™™(s,a) and 6™^^(s,a) lie within the ranges [b~ ^{a) , bf ^{a)'j , respectively. Equating 
the boundary values of b from both the s- and t-channel partial-wave expansions and the s-channel 
partial- wave projection yields 

6r(s,a) = &-,(a) => s = s-,ia), 

br^{s,a) = bUa) => s = stAa), (C.36) 



where s^ ^ and s^^ denote the corresponding maximal solutions for given a, leads to two equations for 
the two wanted unknowns s™f^ and a* ^ defined by 

~sTr = max <r(«) = sTrC<t) ■ (CS?) 

Explicitly, they follow from equating the maximal solutions 

s^,M) = s^tia) => a = alt , s^tio-lt) = st,tiK,t) = s™r , 
stM) =max|s^J"V),'5^1'^V)} . 



'(±) 



(«) = 77 ± V T7 - « - ^s,t(o) ' 



s 




{a) = —{ a^ + Si - 6+ (a) ±J a^ + J^l- 6+ (o) - Aa^^l } , (C.38) 



2 



where for s > s_|_ > S/2 we have s~^ = s^ t ^^^ ^^ ^ addition e.g. a < we have sf^ = s^l . The 
maximum value of the two other (i.e. minimal) solutions for dg^t then yields the highest lower bound 
on s and thus we can write 

s^r = max [s+, s;i~\dlt), <J+^(a^j} for s > s+ and a < . (C.39) 

For the s-channel parts, solving the equations numerically for all allowed a < s+ leads to the fol- 
lowing optimal value of a and corresponding range of convergence in s > s+ = 59.64 M^ = (1.08 GeV)^ 

af = -128.30M^ , s+ < s < sf"'' = 106.09 M^ , bj {d^) = 26860 M^ , 6+(a^) = 34388M^ , 

(C.40) 
in agreement with the unpublished Appendix E of [21j|£3 
For the i-channel parts, this procedure results in 

af = -23.19 M^ , s+<s< sf^"" = 97.30M2 , b-{dt) = 2202M^ , bf{dt) = 5212M^ . (C.41) 

In conclusion, the s-channel constraints are weaker than the i-channel ones, which can also be 
deduced from Fig. [T6l where the situation for df = — 23.19 M^ is shown: for this a the range of b 
limited by &j~(of) and bf{df) for the t-channel partial- wave expansion also lies within the allowed 
range of b for the s-channel partial-wave expansion, and hence for the interval of s given in ()C.4ip 
this range of b covers the interval [6™^'^(s, af ), 6™'^^(s, af )] that is needed for the s-channel partial- wave 
projection. By construction, the resulting family of hyperbolae does cross neither any double spectral 



^^Appendix E of [21] deals with finding the optimal values for the s-channel partial-wave projection of the s-channel 
partial-wave expanded absorptive parts of the HDRs only and follows a similar scheme. The quoted results are sj"'" > 
105 M^ for a^ «-117M^. 
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Figure 16: Allowed ranges of b for s-channel partial-wave projection with a = df 



-23.19 M2 for 



s-channel (left) and t-channel (right) partial-wave expansion. Horizontal lines correspond to b^ (a) 
2202 M^ (solid) and 64+ (a) = 5212 M^ (dashed). 



region nor their z'-reflected boundaries as depicted in Fig. I18( left). and thus (IC.4ip corresponding to 
"v/I™^ = -v/97.30 Mtj = 1.38 GeV constitutes the result for the s-channel partial-wave projection, in 
agreement with |29j r^l 



C.4 t-channel partial-wave projection 



The relation between the range of b permitted by the Lehmann ellipse constraint ()C.33P and the 
corresponding range of convergence in t for the projection of the HDR equations onto t-channel partial 
waves for given a is most easily established on the basis of the squared t-channel scattering angle z^ , 
which must cover the range 



0<zt{t,a,b) 



(t - S + 2a)2 - 46 _ (t - S + 2a)^ 



46 



,2„2 



IGpfqf 



{t-U){t-tN) 



< 1 



(C.42) 



for both the s-channel and t-channel partial-wave expanded parts, since (as discussed after the t- 
channel partial- wave projection formulae (j3.6p ) the integrands are always functions of the real square 



Zj even between the thresholds tj^ and t]\f. Furthermore, < z/ < 1 is not only a necessary condition 
for < Zf < 1 but also equivalent to —1 < zt < 1, which in turn is already sufficient to perform the 
partial- wave projections in our case (cf. the discussion in Appendix [B|). Therefore, the range (|C.42|) 
of z^ constitutes the necessary and sufficient condition not only for the physical region t > t^, but 
for all kinematical regions. Obviously, for t^r < t < t^v and given a, z^ can only be non-negative for 6 
non-negative and large enough. Translating ()C.42p into ranges for 6 while taking care of the signs of 
pf and q^ in the different kinematical regions yields (cf. (|C.30p ) 



br''{t,a)<b<br^t,a) 

bf'^''{t,a) <b< bf''{t,a) 
1 



y tT^ < t < t]sf , 

V t > % (or i < tjr) 



(C.43) 



-(t-S + 2a)2 



"^PtQt 



(t 



Sla + a^ + S^ 



bf^it, a) = -(t - S + 2af > , 6r"(t, a) 

^* As combined result for both s- and f-channel contributions, the numbers s™^'= = 97 M^ for a" — —23 M^ are quoted 
in |29) without further explanation and giving only a vague reference for these numerical values, which is most probably 
meant to be [107] , However, roughly the same numbers are also given more recently in |25j . 
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Figure 17: Ranges of convergence in t for t-channel partial-wave projection from full coverage (shaded 
area) of the physical band < zf < 1 for a = af = -23.19 M^ (left) and a = a\t = -2.71 Ml (right). 
Vertical lines indicate thresholds t.,^ and t^r. 



where the superscripts min/max refer to both the (at least partially) unphysical kineniatical range 
t > t,r needed in our RS system as well as the corresponding min/max values 0/1 of z|. Solving these 
equations for t yields (cf. t(±)(i^ = 0;a,6) of (|2.42p ) 



tf\a,Vl 



min"> 



2a ±2 



^min 



ti{a,bTn 



a + 



imax 



(C.44) 



and the range of convergence in t for given a is the kinematical range in which all values between 
6™™(t,a) and 6™^^(t,a) are covered by both intervals [b~i{a),b~^^{a)\. Between the thresholds (i.e. for 
t^ <t < tj^) this amounts to the conditions h~^{a) < 6™^'^(t, a) and bf^^^{t, a) < 6^j(a), while below or 
above the thresholds (i.e. for t < t^ oi tN < t) we have b~t{a) < bf"^{t,a) and bf"^{t,a) < btj.{a)^ 
Equivalently, we can demand that for given a the band < zf (t, a,b) < 1 must be fully covered by the 
area between zl(t,a,bj^{a)) and zl{t,a,b'^^{a)) in order to determine the range of validity in t. The 
situation that results from using the set ()C.4ip of optimal parameters for the s-channel partial- wave 
projection derived in the previous section is shown in Fig. [TWleft): the t-channel projection is then 
valid for — 5.63M^ < t < 4:4.92 M^ (denoted by the shaded area of coverage), and the reason for this 
rather low upper bound on t is that the curve for 6^ = 2202 M^ changes sign between the thresholds 
and thus enters the critical band < z| < 1, which is hence no longer fully covered by the allowed area. 
Indeed, the range of convergence can be significantly improved if z|(t, a, 6~j(a)) < (and of course 
also z|(t, a, 6j~j(a)) > 1) for all t between the thresholds. From ()C.42p it is clear that for t G (t7r,iAf) 
we have zl{t,a,b) < if and only if 6 < {t — Tj + 2a)^/4, such that the curves for the lower limits 
6~j(a) of b will be tangent to the zero axis provided that b~^{a) = 0. Solving this numerically yields 



b, (a) = 6j (a) = 



ait 



-2.71 M^ 



(C.45) 



■s,t 



as well as bg^{a) < for 



which is unambiguous since it turns out that 6^4(0) > for a < a; 
al ^ < a < s+ (where we have used the numerical equality of the lower bounds for both s- and t- 
channel partial- wave expansion). Furthermore, the curves for bgf{a) start to enter the critical band 
due to change of sign at t^ for a > 2.58 M^. and a > 9.17 M^, respectively (however, the geometrical 



^^Accordingly, at the thresholds the respective min/max values are identical: 
and 6r"(iiv, a) = fer'"'(iiv, a) = (a-f E_)^ > 0. 
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Figure 18: Double spectral regions and limiting hyperbolae for s- and t-channel partial- wave projection. 
Left: for a = af = -23.19 M^ with 67(a) = 2202 M^ (solid) and If {a) = 5212 M^ (dashed). Right: 



for a 



ait 



-2.71 M^ with b-^ia) = (solid) and b+{a) = 2897 M^ (dashed). 



constraint a < tir is partially tighter anyway). Thus, a* ^ is the smallest value of a such that the 
critical band is fully covered between the thresholds, which is shown in Fig. [TTT right). From this 
figure and equation ()C.44p it is clear that in this case we can deduce the corresponding upper and 

™f^(a) on t by the intercepts tj, (a, 6^((a)) of z^{t, a, bft{o.)) with the zero 

axis below and above the thresholds, respectively. Since moreover both t™f^(a) = tg (0,6^^(0)) are 
strictly decreasing in the allowed ranges of a, the minimal allowed value a = a^j. is also the optimal 

one yielding t 



lower bounds f^l^ia) and t 



^s,t 



min 

s,t 



4-)(a* ,&+(«*)) and ef 



s.t 

tp (a* J, 6^j(a* J). This procedure results in 



bfii 



2897 M^ 
3509 Mt 



-9.84 M^ < t< 205.45 M^ , 
-20.67 M^ < t < 216.28 M^ , 



(C.46) 



where the s-channel Lehmann ellipse constraint proves slightly more restrictive, and thus the final 
result for the t-channel partial-wave projection reads 



-2.71 m: 



U <t<f' 



205.45 Mi 



KA) 



0, bt{~alA = 2897M^ 



(C.47) 



which corresponds to vt™^^ = -v/205.45 M^r = 2.00 GeV. Again, ascertaining that the resulting family 
of hyperbolae does enter neither any double spectral region nor their z'-refiected boundaries, which is 
shown in Fig. fTHT right). completes the derivation of the final result ()C.47jl for the t-channel partial- 
wave projection. It is interesting to note that the domain of validity in t is much bigger as the one 
in s, which is refiected by the possibility to use only the positive half < z^ < 1 of the range of 
the scattering angle due to Bose symmetry in the t-channel; in particular the range of convergence 
connects the physical regions for the s- and n-channel reactions, where t < 0, with the t-channel 
physical region t > tjy. 

The complicated interplay between a, bg^{a), and zl{t,a,b) in the different kinematical regions is 
the reason why it is not possible to treat the t-channel projection in analogy to the s-channel projection 
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in the previous section: equating again the corresponding boundary values of b from both the s- and t- 
channel partial-wave expansions and the t-channel partial-wave projection, and subsequently equating 
the corresponding maximal solutions in order to obtain t™f^ as the maximal upper limit on t for t > t]\f 
leads to entering or even crossing the critical band between the thresholds. 



D Asymptotic regions and Regge theory 



The asymptotic s- and t-channel contributions of the HDRs (j2.44p to the invariant amplitudes are 
defined by splitting the corresponding integration ranges s+ < s' < oo and t-,^ < t' < oo at some 
appropriate values Sa = W^ and ta, respectively, which yields the following asymptotic contributions 



+ (s,t) = 

asym ^ ' ' 


ih' 


1 1 

+ - 

s — s s — u 


1 


s' — a 




oo 


(s,t) = 

asym ^ ' ' 


^/-' 


i i 


ImA-{t 


s' — s s' — u 



^ '^ TT J t'-t 

oo 

'^-^ + ,r/ s'-n' t'-t ' 



B+\ (s,t) = - I ds' 

lasym \ ' ^ ' 



B-\ (s,t) = - I ds' 

lasym \ ' ^ yr 




X ..+ / / /N 1 /" w s-^i lmB+(t',z'.) 
ImB+{s',z',) + - I dt' ^ ' *^ 



vr 



Imi? {s , Zg) 



s' -u' t'-t 



1 [ ^^MB-{t',z'^) 



vr 



t'-t 



(D.l) 



and the remaining non-asymptotic parts are given by the corresponding integrals over s+ < s' < Sa 
and t-jr < t' < ta, respectively, plus the nucleon pole terms N^{s,t) for the amplitudes B\s,t). The 
internal (primed) kinematics are given by (cf. Sect. [231 and especially ()2.42|) 1 



'(t'; a, 6) = - S - t' + ^{f - S + 2a)2 - 46 , t'{s';a,b) 



+ T, — s' — a 



s' — a 



1 



n'(t'; a, 6) = - f S - t' - V(t' - S + 2a)2 - 46 j , 

where the parameter 6 is fixed by the external (unprimed) kinematics as 

(s — a)(S — s — t — a) = b = {s — a){u' — a) , 
such that 



s'{u'; a,b) 



b 



+ a 



u'{s'; a, 6) 



a 



(D.2) 



(D.3) 



(D.4) 



u — a s — a 

Thus (for given a and finite 6), for the s-channel integrals we need the asymptotic behavior in the 
limit 

s' — )• oo =^ t' — )• — oo , u' — )• a , (D.5) 

while in the t-channel integrals the asymptotic behavior is determined by 



t — )• oo 



-oo 



a . 



(D.6) 
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Prom (jD.ip the asymptotic parts of the s- and t-channel partial waves may then be deduced by 
the projection formulae (IA.14P and (IB.9P as 

1 

-1 
1 

/(Vd- L,^^ (W) = I dz,{Rl,,{W, zs) A^\^^^ {W, z,) + Rl,,{W^ z,) B^\^^^ {W, z,)] , 

-1 

1 

/+lasym (*) = 0(0 / dz,{uj{t, Zt) A^l^^^ {t, Zt) + Vj{t, Zt) B^l^^^ {t, Z,)} , 



/- lasyn. (^ = 0(0 / dz, Wj{t, Z,) B\^^^^^ {t, Z,) 



(D.7) 



where for the t-channel partial waves we have again / = +/— for even/odd J. Note that for these 
asymptotic contributions we do not expand the absorptive parts inside the integrals in order to take 
into account the high-energy behavior of the full invariant amplitudes as given by Regge theory |108j . 
Therefore, also for the so-called driving terms (i.e. the sums of all higher partial waves that are not 
taken into account explicitly |12yi9] as well as the asymptotic contributions of the lower partial waves 
treated dynamically) the integration ranges are limited by Sa and t^, in order to avoid double counting 
of the asymptotic regions. This procedure follows |121 I19]. motivated by the observations that, first, 
for higher and higher energies one would be forced to explicitly use higher and higher partial waves as 
well in order to ensure the validity of the partial-wave expansion, and second, no available information 
in the asymptotic regime is lost without need. In Sect. 15.3) we explicitly demonstrate the matching of 
the Regge model to truncated sums of the lowest partial waves with I < lma.x for Zmax £ {3,4,5}. 

In the following, the contributions from the asymptotic regions in both channels will be examined 
in the framework of Regge theory. For a general introduction see e.g. |109] . 

D.l s-channel asymptotics 

First of all, contributions from t-channel Regge trajectories, i.e. the leading Pomeron (It = 0) 



trajectory ap{t') ss a 



(0) 



1 (roughly independent of t' but with exponential residue function 



j5p{t') = (Tpexp^l^, where ap represents the asymptotic total-cross-section value for tttt scatter- 
ing and bp is the width of the diffraction peak, cf. [I2l[l9]) as well as the p {It = 1) and / {It = 0) 
trajectories ap{t') = a), + a), t' (and cxf{t') in analogy) should be negligible, since due to (|D.5p they 
will behave as 



Im^(s',t')~;5p(t')s'"^(*'^ 



~ e 2 s ~ e 



Im ^(s',t') 



J a pit') 



Pp{t')s''^'"^- ' ~ S'"" ^""P ~ S 



/a(°'+t'a 



(1) 



(D.8) 
for s' — >• oo, leading to an exponential suppression. 

Let us briefly review the u-channel-exchange contributions to the s-channel reactions of backward 
ttN scattering as discussed in [73]. The invariant amplitudes can be parameterized according to (cf. 
also [291) 



A{s',u') 



^ (if{u')Ci{u') / s' 
^r{a,{u')-^)\sR 



ai{u')- 



b{sW)-e4^^P^(-"'^'' 

^T{a,{u')-'^)\sR 



(D.9) 
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where both sums run over the four trajectories i E {Na,N^, A5, A^}, and the Regge propagators Qiu') 
are given by 

sm[7r[ai{u') - ^l) 

Besides the scahng factor sr = 1 GeV^, the fohowing Regge residues /3j [u') and Regge trajectories 
ai{u') are employed: 

/3f{u') = Oj + biu' , f3f{u') = Ci + diu' , ai{u') = a\' + a'u , (D.ll) 

i.e. both the residues and the trajectories are Hnearly parameterized, and for the latter an identical 
slope a\ = a' is used for all i. The signature Si = (— 1) *~2 of the trajectory i is positive for Na and 
Ap and negative for A'^^ and A^. Since lraQ(u') = —Si, we may conclude that the imaginary parts of 
the invariant amplitudes in the n-channel isospin basis /„ G {1/2 = A^, 3/2 = A} can be written as 

(D.12) 
with the abbreviations 

^fi^') = - w^^^^tF^ > (D-13) 






and analogously for the B amplitudes. Using now the isospin crossing relations (I2.28p . we finally 
obtain the absorptive parts 



I \ai{u')-^ 9 /„'\"i(«')-^ 



1 / c' \ "'I-" ''^2 9 / c' \ ' 

ImA+(.',n'(.',t')) = +^ E A^K) f +1 E A^(^') f 

where the dependence on {s\t') can be translated into dependencies on {s\z'g) for the s-channel 
integrals and {t' ,z[) for the i-channel integrals via (I2.52p . For convenience we also give the numerical 
values of [73] for the 21 real parameters in Table [2j 

As a byproduct, we can use these relations to infer the high-energy behavior of the HDR s-channel 
integrals: from the trajectory parameters given in Table [5] it follows that the high-energy tail of the 
integrals will be governed by the A^ trajectory. Explicitly, the integrands for A^ and i3~ will behave 
as 



.1 ,„/„_i.v is'-'-^^ for a = -23.19 M2 
-'-^•52 for a = -2.71 M2 , 



,'-l,-A,(a)-i ^ ,.'a-lA-r =W ^°^ « = --^-=^-. > (D.15) 



for v! — >• a, whereas the integrands for A and B^ fall off faster by one power in s' (cf. (ID. ID ). We 
thus conclude that the s-channel part of the (unsubtracted) HDRs ()2.44p converges in principle for 

99 



Na N^ As Af 



a [GeV- 


-1- 


-60.68 


47.22 


-75.15 


1419.99 


b [GeV- 


-3" 


326.52 


-215.84 


-138.75 


3052.84 


c [GeV- 


-2" 


546.40 


-101.11 


64.16 


-192.64 


d [GeV- 


-4- 


307.42 


-128.04 


86.77 


-695.81 



a(0) -0.36 -0.62 0.03 -2.65 



a 



' [GeV-2] 0.908 



Table 2: Regge-model parameter values for backward ttN scattering as given in |73j . 

a < 26.57 M^. Note that in order to investigate the behavior of these asymptotic contributions in 
the "fixed-t limit" a — )• — c« (as discussed in Appendix IC.2P it is important to take the limits in 
the correct order, since u' — )• a only after s' — )• oo. Since ai{u') — 2 < —1 for sufficiently large and 
negative a, the s-channel Regge contributions vanish in the limit s' — )• 00 for such values of a. As 
shown in Sect. 15.41 these asymptotic contributions are numerically small for the optimal value of a 
(and a reasonable choice of Sa), and thus they can be safely neglected for a — ;■ — c«, regardless of the 
pathological behavior of the Regge model due to the Gamma function in this case. 

D.2 f -channel asymptotics 

Similarly to the previous section one could use Regge theory to describe the t-channel asymptotic 
region. However, the significance of these contributions in view of the corresponding low-energy 
region differs strongly from the s-channel: while contributing crucially to the dispersive integrals, 
the pseudophysical region t-^ < t < t]\f cannot be constrained from experiment, but requires an 
analytic continuation. Within our system of RS equations this task naturally takes the form of a MO 
problem, as explained in Sects. [3] and [5l The solution of these equations becomes rather involved once 
intermediate states other than tttt are energetically allowed, which happens around 1 GeV (especially 
KK above 2Mk)- In view of the ensuing uncertainty of the t-channel partial waves even below the 
NN threshold it is clear that the inclusion of phase-shift solutions above t^v [70], and even more so the 
modeling of the high-energy region, will be of little practical relevance. Moreover, as shown explicitly 
in Sect. 15.41 already the s-channel Regge contributions are numerically immaterial, in particular if 
subtractions are performed, which provides evidence that also the high-energy region in the i-channel 
can be safely ignored. For these reasons, we will not consider the t-channel asymptotic region any 
further. 

D.3 Subtracted asymptotics 

Here, we show how to incorporate the effects due to subtractions into the Regge description of the 
asymptotic parts of the corresponding subtracted HDRs (14. lip and (14.12J) . However, according to 
Appendix ID. 2 1 all asymptotic t-channel contributions will be neglected. 

For the high-energy tail s' > s^ of the s-channel integrals, according to Appendix ID. li the absorptive 
parts may generically be written as sums of Regge-trajectory contributions 

lmX^-{s',u'{s',t')) = ^^f{u')(—)^ ' ior X£{A,B}, (D.16) 
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with summands of the generic form (i.e. dropping the indices X and i for the time being) 



hw) 



Sf3{u') 



T{aiu') - i) 



l^ ' 



f3{u') = f3^'^^ + /3(i)n' , a(n') = q(o) + a'u' 



(D.17) 



While the evaluation of the Regge contributions is straightforward in the un- and once-subtracted 
case, for two subtractions one furthermore needs the derivative 



dt Ku') 
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-1 
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SR 
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■«}J 



(D.18) 



(0,0) 



where ^{z) denotes the digamma function defined as the logarithmic derivative of the gamma function 

^(z) = ±logr{z) = ^^. (D.19) 



To this end, one may use u'{s' , t') = S — s' — t' and (cf. (j4.10p ) 

is' -so? 



[f] 



(0,0) 



[9tt'] 



(0,0) 



■So - Q 
s' — a 



(D.20) 



After utilizing the crossing relations in order to rewrite the Regge contributions in the / G {+, — } 
isospin basis and expressing t' as well as the corresponding kernel functions in terms of (s', z^), we can 
perform the partial- wave projections of the s-channel contributions onto both s- and f-channel partial 
waves according to ()D.7p . where again the implicit kinematical dependencies have to be taken into 
account accordingly. 

Finally, we demonstrate the projection onto the lowest t-channel partial waves with J < 2 explicitly. 
The n-times subtracted versions of (ID.7P immediately lead to 







ip-sub 1 /■ 2. 

•' + lasym ^ ^ 47r J 




lasym AptqtZt 
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A \„^„^\l,Zt) „_in-sub/^ s\ 



4^2- lasym 
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A^\:-^it,z,)+4mq^zr '^^^°^ 



lasym 



^Ptqtzt 
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(D.21) 
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again written in terms of quantities that are always real since AptqtZt = ^mv. Here, the asymptotic 
s-channel contributions to the invariant amplitudes for e.g. the twice-subtracted case read (i.e. as 
functions of {t,zt), cf. ()B.14p for z'g{t,s';zt)) 



^+|2-sub 1 f,J[ 2{s'-so)+t 1 



s-asym^ ' '^ TT J \ [(g' - Sq + ^)^ - ^pUht «' " « 

t 



1 2-sub 




\niA^{s, Zg 



) [lmA+(.',4)](g_Q) - hois') t[dtlTnA+{s',z';, 



^ ls-asym(^'^*) _ 1 f ^j j lu. A' {s\ z',) [Ini A' {s' , Z^)](o q) 

AptqtZt TT J \{s'-SQ + ^Y-Aplq^z'l {s' - SqY 



(D.22) 



and analogously for B~{t,zt) and B^{t,Zt)/{AptqtZt)- Note that again only real squares of momenta 
and zt occur and hence these formulae are valid in all kinematical regions. Furthermore, by rewriting 
the general t-channel partial-wave projections ()3.6p for both even and odd J in terms of real quantities 
(i.e. i/-even amplitudes and squares of momenta as well as squares of zt) as above, the partial waves 
exhibit ostensible poles at t-,^ for all J > 2 and in addition at t^ for all J > 3, while from the 
discussion of their threshold behavior in Sect. 13.3.11 we know that these poles are immaterial. The 
reason for this behavior can be understood by first noting that for ptqt — ?• the asymptotic (s-channel) 
contributions ()D.22p no longer depend on zt- The orthogonality of the Legendre polynomials Pj{zt) for 
even J > 2 and odd J > 3 then balances the poles and leads to the expected finite (but non-vanishing) 
values of the partial waves at both the pseudothreshold tj^ and the threshold t^ (cf. the explicit case 
for flit) in dEHI). 
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